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Studies of the ;N"(, a) sB™ reaction indicate the emission of heavy particle groups possibly 
as a result of resonance transmutation by fast neutrons from a Ra— Be source. When a nitrogen 
absorber is placed before the nitrogen detector, the heavy particle groups apparently do not 
appear. Upon interposition of carbon as a retarding material between nitrogen’ absorber and 
detector, the reappearance of some of the groups in question occurs. This behavior would 
appear to indicate that nitrogen has a larger than normal probability to undergo transmutation 
by certain fast neutron groups. Two interpretations suggest themselves: (1) the experimental 
points, while based on more than 100,000 pulses, possess statistical fluctuations sufficient to 
produce the observed apparent resonance phenomena for fast neutrons or (2) nitrogen has a 
fast neutron absorption cross section of the order of 60X10-* cm*, much larger than the 
theoretical upper limit of about 0.7 10-* cm*. The resonance interpretation is supported by 
supplementary evidence. Further development of these techniques should make it possible to 
(1) estimate widths of nuclear levels, (2) measure nuclear energy levels, (3) check the reality 
of observed groups, and (4) determine which isotope is responsible for a group. Additional 
investigation of the subject is indicated. 























I, INTRODUCTION energy distribution of the resulting particles, 


N a previous study by the authors! of the when measured, shows peaks. There is evidence 
transformation ;N"(n, a)sB" by fast neutrons that such emission resonance correlates with an 
of continuous energy, discrete groups were found absorption resonance of the neutrons.’ Accord- 
for the emitted heavy particles. The method ingly, an absorber of resonance neutrons, inter- 
employed was that of Wilhelmy,? in which mediate between source and chamber, should be 
neutrons from a continuous source irradiate the responsible for the absence of heavy particle 
target gas in an ionization chamber which is groups. To test if resonance absorption of fast 
coupled to a linear amplifier and recording oscil- neutrons could be detected, the following series 
lograph. If discrete groups are emitted, the of experiments were performed. The apparatus 
t Work done in 1939-1942, while the authors were at and ar used were the same as that previously 
New York University. described,! and the nitrogen pressure in the 
* Major, Corps of Engineers; Chief, Operation Units. sanivats 
1948) i Zagor and F. x. Valente, Phys. Rev. 67, 133 ionization chamber was 336 cm Hg. 
SE. Wilhelmy, Zeits. f. Physik 107, 769 (1937). * W. Maurer, Zeits. f. Physik 107, 721 (1937). 
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Fic. 1. Upper curve is a plot of number of pulses against 
Ea,s for yN“(n, a)sB" reaction. Dotted lower curve is 
plot of number of pulses against pulse size, with NH,NO; 
as absorber. Curves at left, } scale. 


II. RESULTS 


The upper curve in Fig. 1 shows the result of 
a plot of the sum of the energies of disintegration 
products of nitrogen against the number of 
associated oscillographic pulses. The curve shows 
several principal maxima at the energies 1.68, 
2.15, 2.64, 2.98, 3.82, 4.14, and 4.48 Mev. 
Previous experiments indicate these groups to 
be heavy particle groups arising from the fast 
neutron reaction 7N"“(n, a) 5B". 

The experiment was repeated with NH,NO; 
and a nitrogen thickness of 0.38 g/cm?* as an 


absorber. The lower curve in Fig. 1 shows the _ 


result that the distribution is smooth; all the 
groups have faded out. 

Al(NO;)3-9H.O, having a nitrogen thickness 
of 0.26 g/cm?, was used next as the intermediate 
absorber and the experiment repeated. The 
resulting distribution is shown in Fig. 2, and it 
is seen that the distribution is smooth with only 
evidence for the 2.64-Mev group remaining. 
However the intensity of this group is markedly 
diminished. 

To eliminate the possibility that the fading 
out of the groups was caused by the H;0 or Al 
in the Al(NO;)3-9H,O and not to the nitrogen, 


* The background averaged about one count per minute 
and showed a random distribution. It was thus disregarded 
in the calculation. 
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5 See Figs. 4 and 5, reference 1. 











I. ZAGOR 
the distribution was remeasured first with an 
H:.O absorber, and then with an Al absorber. 
The thickness of the H:O absorber in the latter 
measurement was 1.4 g/cm? as compared with 
the 1.1 g/cm? thickness of H,O in the original 
Al(NOs3)3-9H:0 absorber. 

The result with the H,O absorber is shown in 
Fig. 3. The distribution is not smooth and the 
1.68-, 2.15-, 2.98-Mev groups show up clearly, 
though they are not sharp. The reason for their 
lack of sharpness is that the incident neutrons 
suffer appreciable slowing in passing through the 
H,0. This condition is similar to the case in 
which a paraffin plug is placed between source 
and chamber’ and one would expect a diminution 
sharpness of the peaks. 

The thickness of aluminum used in remeasur- 
ing the distribution was 0.91 g/cm? as compared 
with 0.17 g/cm? aluminum in the original 
Al(NOs3)3-9H:O absorber. The resulting dis- 
tribution in Fig. 4 shows the 1.68-, 2.15-, 2.64-, 
and 2.98-Mev alpha-groups appearing sharply. 

Finally, commercial cyanamid (CaCNNH),) 
with a nitrogen thickness of 0.5 g/cm? was used 
as absorber. The distribution, shown in Fig. 5, 
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Fic. 2. Al(NO3)3-9H:O absorber between source and 


chamber. 
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Fic. 3. H2O absorber between source and chamber. 


is completely smooth, with all the groups being 
washed out. 

It thus appears that the interposition of a 
nitrogen absorber between the source and 
nitrogen chamber causes the groups to disappear, 
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Fic. 4. Al absorber between source and chamber. 


RESONANCE WITH NITROGEN 












































i 
=o 8 2 3 4 














5S MEV 6 


Fic. 5. CaCNNH:z absorber between source and 
detector. 


whereas non-nitrogeneous absorbers do not 
cause the groups to fade out. Since the groups 
detected are heavy particles arising from a fast 
neutron reaction, the evidence would appear 
to indicate that resonance absorption of fast 
neutrons has been detected. 

The work of Hornbostel and Valente*® sug- 
gested that the placing of a slowing or moderating 
material between the nitrogen absorber and the 
nitrogen chamber should change the neutron 
spectrum emerging from the nitrogen absorber 
by shifting the spectrum so as to replace, at 
least partially, the absorbed neutrons. The 
neutrons striking the slowing material will have 
a distribution corresponding to a spectrum with 
absorption bands in it. The effect of the retarding 
material is to shift (by retardation of the neu- 
trons) the spectrum, thereby modifying the dis- 
tribution so that the absorption bands are at 
least partially reduced and the distribution 
assumes more nearly its original condition. If the 
proper thickness of retarder is used, some of the 
groups would be expected to reappear. 

A carbon slab of thickness 1.27 g/cm? was 


6 J. Hornbostel and F. A. Valente, Phys. Rev. 55, 108 
(1939). 


























interposed between the NH,NO; absorber and 
nitrogen chamber. The resulting distribution was 
measured, and is shown in the lower curve in 
Fig. 6. The 1.68-, 2.15-, and 2.64-Mev groups re- 
appeared sharply, whereas the remaining groups 
do not. Since all the groups did not reappear, the 
thickness of carbon slab used must have been 
such as to cause only a partial replacement of the 
absorbed neutrons. The upper curve is the ni- 
trogen distribution curve without carbon re- 
tarder and has been shown previously as the 
upper curve in Fig. 1. 


III. DISCUSSION 


The experimental results suggest that we may 
be observing resonance absorption of fast 
neutrons. Interposition of a nitrogen absorber 
between source and nitrogen detector results in 
the disappearance of the groups, and this effect 
is independent of the kind of nitrate used. When 
subsequently a graphite retarder is placed 
between the nitrogen absorber and nitrogen 
detector, some of the groups reappear. The inter- 
position of absorbers, other than paraffin,’ 
between the nitrogen absorber and nitrogen de- 
tector do not cause the disappearance of the 
detected groups. The interposition of a paraffin 
plug between source and detector results in the 
non-appearance of the groups,® indicating the 
observed phenomena to be produced by fast 
neutrons. That instrumentation may not be the 
cause of the groups is shown by the relatively 
smooth curve obtained with hydrogen.’ 

Nevertheless, and despite that more than 
100,000 pulses were counted in the nitrogen-filled 
chamber, the possibility of describing the data 
by a smooth curve, drawn through the experi- 
mental points, cannot be excluded entirely. 

Hansen,® in studying the range velocity rela- 
tion for boron, also investigated the reaction 
7N"(n, a)sB" with D+Be neutrons. From his 

data Hansen questions whether the maxima he 
observed have physical meaning or are caused 
by statistical fluctuations. As a check, he applied 
the x? test and found ‘that if the points were 

- represented by the smooth curve, then the 

chance that the deviations were as large as those 

7 See Fig. 8, reference 1 


= Hansen, Ph.D. thesis (Yale University, 1941), 
p. 
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Fic, 6. Upper curve is distribution with no absorber. 
(Same as upper curve, Fig. 1.) Dotted lower curve is dis- 
tribution with carbon slower placed between NH,NO; 
po ag and nitrogen ionization chamber. Curves at left, 

scale. 


observed was about 3 to 10.” He concludes that 
this does not prove that the smooth curves are 
not the best representations of the points. 
Hansen points out that although it is unlikely 
that an event which has three chances in ten of 
happening should happen in six consecutive runs, 
a possible interpretation might be that all six 
runs are equally bad. 

In this connection, it should be pointed out 
that the value of the nitrogen cross section, as 
deducted from experiment, is roughly 100- fold 
greater than the maximum theoretical cross 
section, though the character of the resonance is 
otherwise qualitatively in accord with the Bohr 
theory. 

The maximum theoretical value of the cross 
section for a reaction of a neutron with a nucleus 
is given by the larger of the two expressions, 
namely : 


x (geometrical radius of nucleus)’, 
or 
@ (neutron wave-length divided by 7)’. 


For 1-Mev neutrons on nitrogen, both these ex- 
pressions are of the order of magnitude of 
0.7 10-* cm?. From this value, the thickness of 
nitrogen absorber required to make the resonance 
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groups disappear is calculated as follows: 


thickness in g/cm? 
atomic weight of nitrogen 





~ cross section X Avogadro’s number 
= 14/0.7 X 10-**X0.6 X 10*4~ 33 g/cm?. 


In contrast to this theoretical value, the ex- 
perimental figure for the thickness required to 
observe the apparent disappearance of the 
resonance groups in question was about 0.38 
g/cm*. 

In view of the importance to theory of this 
question of interpretation, additional investiga- 
tion of the subject is indicated. 


IV. APPLICATIONS 


Resonance absorption can be applied to check 
the reality of observed groups. If the observed 
groups can be caused to fade out by inter- 
position of the appropriate resonance absorber, 
it is reasonably certain that the observed groups 
are real and are not the result of instrumentation. 

In the case of isotopes, the use of an enriched 
isotope resonance absorber might permit deter- 
mining precisely which isotopes are responsible 
for particular groups. 

The method of interposition of slower between 
resonance absorber and ionization chamber can 
be applied to estimate experimentally the widths 
of nuclear levels. To do this, the experiment 
would be repeated with varying thicknesses of 
scattering material between resonance absorber 
and detector. As the thickness of the scattering 
material is increased, the groups should begin to 
reappear, but not uniformly. The assumption 
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made is that narrower levels will yield groups 
more easily detected than broad levels under 
these conditions. In particular, the minimum 
thickness of scatterer at which a particular group 
reappears is noted. From a knowledge of the 
neutron-scatterer cross section, the width of the 
level can be estimated. However considerable 
difficulties due to the large number of corrections 
necessary are involved in this method. 

In addition, energy levels in nuclei could be 
measured by an extension of the method of 
resonance absorption by the following procedure. 
The element to serve as detector is placed into 
the ionization chamber. The experimental pro- 
cedure discussed in the text is followed, and the 
nuclear levels measured, each heavy particle 
group being associated with a nuclear level. The 
element whose nuclear levels are to be measured 
is then placed as absorber between detector and 
source and the energy distribution of the de- 
tector is redetermined. If any of the levels in 
the absorber overlap a level in the detector, then 
the group corresponding to that level should 
vanish. This is an extention of the method em- 
ployed by Salant, Horvath, and Zagor® for slow 
neutrons. A limitation in this application is that 


- the absorber must have overlapping levels in the 


same energy region as the detector or else 
resonance absorption will not be detected. How- 
ever, this might be overcome by the inter- 
position of various thicknesses of scatterers 
between absorber and detector in order to shift 
the neutron energy spectrum entering the 
detector into the proper energy region. 


*E. Salant, W. Horvath, and H. I. Zagor, Phys. Rev. 
55, 111 (1939). 
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A description is given of an alpha-ray spectrograph, con- 
sisting of the Princeton cyclotron magnet and a Plexiglas 
deflection chamber, in which the alpha-particles can be bent 
into a semi-circle of about 80-cm maximum diameter. Three 
different methods of detection have been employed ac- 
cording to the different strengths of the radioactive sources. 
They are the ordinary photographic method, the counting 
method, and the method of photographic tracks. The 
behavior of the spectrograph has been investigated with Po 
alpha-particles by the three methods of detection. The 
forms of the energy distribution have been determined, 
respectively, by these three methods and agree fairly well 
with one another. The half-width of the main line under 
favorable conditions is less than 4 mm, which is equivalent 
to about 0.01 Mev. The intensities a few millimeters away 
from the maximum are less than 0.1 percent of the main 
line intensity. The track method of detection has been used 
to study both the low and high energy regions of the main 
line of the Po alpha-particles. Microscopic examination 
reveals distinctly a series of weak groups in the low energy 
region, while in the high energy region no indication of any 


discrete group has been found. These experiments have 
been carefully repeated with three different sources and 
under different experimental conditions. In each case 
similar results have been obtained. The emission of these 
different alpha-particle groups from Po may, as in the case 
of ordinary fine structure, leave the product nucleus, Pb, 
in different excited states. Then, when the Pb nucleus falls 
from one of these quantum states to the normal state, a 
gamma-ray quantum is given off. As far as energy and 
intensity are concerned, these weak groups of alpha- 
particles are generally compatible with the gamma-ray 
lines from Po as measured by Bothe. However, their in- 
tensities are not in agreement with the current theory of 
alpha-decay, as can be seen from the abnormally large spin 
changes deduced from the theory for the different alpha- 
transformations, and also from the large deviation of the 
Geiger-Nuttall curve for these groups from the curve for 
the members of the Ra family. Possible explanations of 
these short-range alpha-particles from Po have been 
suggested and discussed. 





I. INTRODUCTION 


determine the energy of a particle coming 
out from a nucleus is one of the few ways 
to approach the problem of structure of the 
nucleus. From the energies of the particle-groups 
an energy level system can be deduced for the 
nucleus. It has been shown elsewhere that there 
are indications of regularities in the energy 
levels of both light! and heavy? nuclei. Since the 
data are rather insufficient at the present time, 
the conclusions so far drawn are only provisional. 
However, when the results become more com- 
plete and accurate, a more exact relationship 
between the energy levels and certain properties 
of the nucleus may be obtained, which will 
certainly serve as an empirical basis for theo- 
retical investigation. Therefore, it is important 
to determine accurately the spectra (alpha-, 
beta-, and gamma-rays) of radioactive nuclei as 
1M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 9, 
303 (1937); W. Y. Chang, Phys. Rev. A65, 352 (1944); 
K. M. Guggenheimer, Proc. Roy. Soc. A181, 172 (1942). 
? Rutherford and Ellis, Proc. Roy. Soc. A132, 667 (1931); 


Rosenblum, Comptes rendus 202, 943 (1936); Rosenblum 
and Guillot, Comptes rendus 204, 975, 1727 (1937). 


60 


well as the energy spectra of particles emitted 
in nuclear reactions. 

The best method so far known for measuring 
accurately the energies of alpha-particles is the 
semi-circular magnetic focusing method, the mo- 
mentum of the particle being calculated from 
the radius of the circular path and the value of 
the applied magnetic field. In order that the 
magnetic spectrograph may have a very high 
resolving power, the radius of the semi-circular 
path in the uniform magnetic field must be very 
large. This means that a huge magnet must be 
designed and installed. Fortunately it has been 
possible to make use of the Princeton cyclotron 
magnet for this purpose. We have been able to 
construct a Lucite deflection chamber in which 
the diameter of an alpha-particle path can be as 
large as about 80 cm. Since the gap between the 
two pole-pieces of a cyclotron magnet is usually 
greater than 10 cm, the solid angle of detection 
in the spectrograph so constructed is much 
larger than usual. In our case about one particle 
out of every 5000 particles emitted can reach 
the detection instrument. 
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Fic. 1. The schematic diagram of the deflection chamber. It shows the focusing action 


From some natural radioactive substances, one 
expects that alpha-particle groups of extremely 
small intensities may be emitted, for in one or 
two cases very weak gamma-rays have been 
observed. If an extremely sensitive method 
could be devised to measure the energy groups 
of these alpha-particles, additional information 
might be obtained about the actual mechanism 
of the alpha-decay process. During the last few 
years the photographic track method has been 
well developed and applied with certain success 


_ to the study of scattering, artificial disintegra- 


tion, and cosmic rays. This method registers 
every particle, and when it is used in conjunction 
with an alpha-ray magnetic spectrograph, there 
is practically no ambiguity in deciding the origin 
(and nature) of the heavy charged particles. 
This method and the use of Eastman fine-grain 
alpha-particle plates, have been shown in our 
case to be very satisfactory in recording the 
weak alpha-particle groups, especially in the 
low energy region with respect to the main 
group where the background is large. Using this 
method, we have found 12 very weak alpha- 


on two groups of alpha-particles, and the arrangement of the chamber slits. 





particle groups from Po, the intensities of which 
do not seem to be compatible with the current 
theory of alpha-decay. Their energies and in- 
tensities are, however, generally compatible with 
the gamma-ray lines from Po determined by 
Bothe. 

Since 1936 very little work has been done on 
the alpha-ray spectra. Rosenblum and his ce- 
workers* have redetermined the relative intensi- 
ties of the fine-structure groups of radioactinium 
and its disintegration products. Several new 
groups have been found. A 60° magnetic focusing 
spectrograph has been constructed by Roy 
Ringo.‘ It may prove useful for the study of 
particles emitted in artificial disintegration as 
well as for the investigation of alpha-ray spectra. 
It has the advantage of portability and low cost 
because of the use of a much smaller magnet. 
But for some reason or other, the energy dis- 
tributions obtained with this instrument do not 
seem to be satisfactory. Data on the alpha-ray 


3 Rosenblum, Guillot, and Percy, Comptes rendus 202, 
1274 (1936); 204, 175 (1937). 
* Roy Ringo, Phys. Rev. 58, 942 (1940). 
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Fic. 2. The vertical section of the deflection chamber between the two pole pieces of the cyclotron magnet. It shows 
the vacuum seals by neoprene strips and the detailed construction of the system containing the source-box and the 
support for the counter and camera. A—defining slit; B—quartz rod; C—lead screening plate; D—support for camera 
or counter ; E—connecting wires (7) for counter (each has a slide contact below the movable support “D’’) ; F—binding 


posts (7) for counter; G—turning head (one revolution = 50/1000”); H—Lucite chamber wall; 


neoprene strips; K—magnet pole pieces. 


spectra investigated prior to 1936 can be found 
in the books of Gamow and Rasetti, respectively, 
entitled “Structure of Atomic Nuclei and Nu- 
clear Transformations” and ‘Elements of Nu- 
clear Physics.”’ 


Il. THE SPECTROGRAPH 


The magnetic spectrograph consists of three 
main parts, namely, the deflection chamber, the 
magnet, and the detection instruments. 


A. The Deflection Chamber 


The deflection chamber of the spectrograph is 
made of transparent Plexiglas in which the alpha- 
particles can be bent by the magnetic field into a 
semi-circle of about 80-cm maximum diameter. 
Figure 1 shows the schematic diagram of the 
chamber. It has the form of a trapezoid with the 
legs having the curvature of the pole pieces of the 
cyclotron magnet. The over-all length (including 
the thickness of the walls) of the longer side is 
89 cm, of the shorter side 60 cm, and the over-all 
distance between these two sides is 53 cm. The 
thickness of each of the four walls is 5 cm and the 
height of the chamber 11 cm (cf. Fig. 2). These 
four walls were joined together with Lucite 
cement. A Lucite plate with a small hole at the 
center (so as to have equal pressures on both 
sides of the plate) was similarly cemented to 


—iron plates (2); J— 


the quadrilateral frame as the bottom of the 
chamber, and aluminum plates were used as 
the cover on the top. Particles from the source 
of the same velocity will complete semi-circles 
of the same radius at a line in the focal plane 
—the plate. In Fig. 1 the paths of two different 
groups of alpha-particles are indicated. It is 
seen from the paths of each group that the 
semi-circle (A or B) passing through the center of 
the defining slit hits the photographic plate, 
which is at the continuation of the defining slit, 
at a point farthest from the source, while all the 
other circles strike the plate at points nearer to 
the source with the former point as a limit. This 
gives rise to a sharp edge at the high energy 
side® of a fairly narrow line image. Each of the 
three chamber slits is also made of a Lucite plate 
which has a rectangular hole for the slit to allow 
the alpha-particles to pass through. The shapes 
of the edges of each hole are such that alpha- 
particles cannot hit the photographic plate after 
single scattering from the edges. These slits serve 
to reduce scattering from the top and the bottom 
as well as from the walls of the chamber. In fact, 
the primary advantage of the use of Plexiglas is 
that the effect of scattering is small. This was 
shown by the fact that the background of an 

5 For theoretical analysis, see K. T. Li, Proc. Camb. 


Phil. Soc. 33, 164 (1937); J. L. Lawson and A. W. Tyler, 
Rev. Sci. Inst. 11, 6 (1939). 
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ALPHA-PARTICLES FROM PO 63 


alpha-ray line as detected by the track method 
(cf. C) was very small and was not much in- 
creased when the chamber slits were not used. 
The defining slit and the source slit (not shown) 
are made of aluminum. Tests have also been 
made by using brass slits; the background then 
seemed to be larger. 

Figure 2 shows the vertical section through the 
deflection chamber which has been put between 
the two pole pieces, K, of the cyclotron magnet. 
Two iron plates, J, of about the same area as the 
pole pieces and about three cm in thickness are 
placed above and below the chamber, respect- 
ively. These plates are used to stand the force 
caused by the pressure difference produced when 
the chamber is evacuated. Between each of the 
iron plates and the wall of the chamber is placed 
a pair of neoprene loops, j, to make the system 
vacuum tight. The interspace between the two 
loops of each pair can be evacuated too. The 
source box and the detection instrument form a 
single unit which can be easily taken out through 
a hole in the chamber wall. A source in the form 
of a long narrow strip (maximum length about 
6.5 cm) can be mounted at a distance of 2.85 cm 
behind the defining slit. If the source is wide, it 
can be inclined firmly at any angle, so that the 
projection of its width on the defining slit is still 
only a few tenths of a millimeter. Thus a better 
line-shape can be secured. The source box is 
rigidly connected by means of the quartz rod, B 
(44 cm long and 2.7 cm in diameter), to the 
support, D, for a camera or a system of counters. 
The distance of the source relative to the detec- 
tion instrument then will not be changed ap- 
preciably when the room temperature varies 
owing to the small thermal coefficient of expan- 
sion of the quartz rod. The seven bronze wires, E 
(one for the common counter voltage), are used 
to connect six counters (cf. C) in a box set on the 
support D to the seven binding posts, F, outside 
the chamber. Since the support D has to be 
moved when a system of counters is used, each of 
the seven wires has a slide contact below D. The 
support, D, is moved through a distance 50/1000 
inch for each complete revolution of the turning 
head, G. The head, G, is turned through a Wilson 
seal® so as to prevent leaks into the chamber. 


*R. R. Wilson, Rev. Sci. Inst. 12, 91 (1941). 





This Wilson seal and the binding posts, F, are 
fixed by using rubber seals on a Lucite plate, 
which covers the hole in the chamber wall 
tightly by means of the neoprene loop, j. To 
prevent the hole from being collapsed by air 
pressure, a Lucite block is used which is just big 
enough to slide freely in and out of the hole. All 
the above parts form a single rigid system, which 
is confined to slide only on a pair of tracks parallel 
to its length, and can then be taken out easily 
through the hole. 

A large hole at the middle of the shorter side 
wall (cf. Fig. 1) has been provided (having a 
seal with a brass disk pressed against a neoprene 
loop) for measuring the magnetic field inside 
the chamber when necessary. Two small tubes 
on the opposite sides of this hole are inserted 
into the chamber wall, so that each tube can be 
used for pumping the air out of the interspace 
between the two neoprene loops of each pair. A 
small hole in the “‘leg’’ side wall near the detec- 
tion instrument is used to let the air into the 
chamber. This is done by turning a rod through 
the hole in order to loosen the Wilson seal inside. 
A system of three shutters over the source box is 
controlled through a Wilson seal by a turning 
head outside the deflection chamber. These three 
shutters are used to stop the beam of alpha- 
particles from coming out when necessary and are 
for three different positions of the source box. A 
large hole provided with a Wilson seal in the 
“‘leg’’ side wall near the shutters is connected to a 
high capacity oil diffusion pump. A Cenco- 
Hypervac 20 is used as a fore-pump. The pressure 
inside the chamber is always maintained lower 
than 10-* mm of Hg, so that scattering from the 
residual gas is negligible. Consequently the 
chance for the capture and loss of electrons’ by 
the alpha-particles in their paths is small too. 

The dimensions and hence the resolution of 
the present deflection chamber are larger than 
the one used by Rosenblum® before 1938 and 
about the same as the annular magnetic spec- 
trograph of Rutherford and his co-workers.’ 
However, in the case of the annular magnetic 
spectrograph, alpha-particles of different ve- 

7E. Rutherford, Phil. Mag. 47, 277 (1924). : 

® Rutherford, Wynn-Williams, Bowden, and Lewis, 
Proc. Roy. Soc. A139, 617 (1933); Rosenblum, Actualités 


scientifiques et industrielles (Exposés de Physique Thé- 
orique, edited by M. L. de Broglie), XX XI (1032). 





MAGNET 
CURRENT 


2az 0-200 we 


THE GAL, cilecurr IS PUT INA 
SEPERATE COPPER EF SHIELD 
Bok FOR THERMAL INSULATION, 





26x. 


1,000.1. 


Sv DRY CELL 


FIELD COIL OF LITTLE GENERATOR 
ARMATURE OF LITTLE GENERATOR 

















To cou. oF 
MAGNET 





ARMATURE OF LARGE 
GENERATOR 








FIELD COIL. OF LARGE 
© GENERATOR 


AUTOMATON 
SESe———————™==— 


Fic. 3. The automaton for controlling the magnet current. 
It shows the principle of operation. 


locities were brought into the ionization cham- 
ber by varying the magnetic field, so that at 
one field only one velocity could be studied. 
In our case, however, particles differing in energy 
by 3 Mev or more can be studied without vary- 
ing the magnetic field. Moreover, the solid angle 
of detection in our case is about 1 in 5000,° 
while in the annular case it was only about 1 in 
20,000. 


B. The Magnetic Field 


The description of the design and installation 
of the Princeton cyclotron magnet was published 


by Henderson and White.!° Here we only mention. 


a few things which are particularly important in 
our case. 


(1) The Mapping of the Magnetic Field 


The spatial distribution of the magnetic field 
at a value of about 11 kilogauss was mapped 
along two perpendicular diameters with a cylin- 
drical search coil and a fluxmeter. The search 
coil has a height about equal to the length of the 
gap between the two pole pieces and a diameter 


* Estimated by counting the number of alpha-particles 
from a known strength of Po source and also checked by 
geometrical calculation. 

10 M. C. Henderson and M. G. White, Rev. Sci. Inst. 9, 
19 (1938). 


of about 5 cm. The variations along the two per- 
pendicular directions are quite similar to one 
another and are only from 1 to 2.5 gauss, until 
a region beyond 37 cm from the center of the 
pole pieces is reached where the field drops 
rapidly toward the edge. This is consistent with 
the results of Henderson and White. In our case 
the radii of the alpha-particle paths are in most 
experiments smaller than 32 cm. Therefore, for 
this working region no effort is necessary to 
correct the variation of the edge field. 


(2) The Control of the Magnetic Field 


An automaton, which has a principle very 
similar to the automaton built by Henderson and 
White’ and which was constructed for the former 
cyclotron work by Mr. R. Kamm, to whom we 
are thankful, has been used with some modifica- 
tions to control the magnet current. Figure 3 
shows the arrangement. A small voltage (about 
200 mv) from a shunt in the field coil is applied 
to the galvanometer circuit and balanced by the 
voltage of the dry cell. The photo-cell-2A3 circuit 
is so arranged that when the magnet current 
increases relatively to the dry cell voltage, the 
galvanometer coil is deflected in such a direction 
that the light beam reflected from its mirror 
falls more into the upper photo-cell. This makes 
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the grid bias of the 2A3 more negative and hence 
the anode current smaller. Since this anode 
current passes through the field coil of the small 
generator, which supplies the field current of the 
large generator, its decrease thus reduces the 
magnet current. Similarly, when the magnet 
current decreases, more light falls into the lower 
photo-cell, and consequently it is increased again. 
The sensitivity of the twin photo-cell can be 
varied by the two extreme slide contacts, while 
the grid bias of the 2A3 can be independently 
changed by the grounded slide contact. The 
sensitivities of the galvanometer, of the twin 
photo-cell, and of the 2A3 have been so adjusted 
that a field of about 12 kilogauss, which is about 
the working field in the future, can be maintained 
constant within 1 in 10‘ for many hours as shown 
by the indicating coil (see below). Under such 
circumstances the range of control is only within 
about 80 gauss. This means that when, because of 
the variation of the field, a Po alpha-ray beam is 
fluctuating more than 5 mm" in the focal plane 
of the deflection chamber, the light beam will 
be thrown out of the twin photo-cell, and the 
control will be then out of action. When this 
happens (very rarely), the magnet current in- 
creases rapidly to its maximum—about 350 


amperes—and never goes back to its normal ° 


value again. This latter fact is very useful, how- 
ever. For, when too large fluctuation of the 
magnet current occurs, one is perfectly sure that 
the photographic plate has been ruined and must 
be discarded. Hence there is no possibility of 
being misled in analyzing the results of the 
spectra. A safety relay system™ has been used 


‘across the input voltage of the field coil of the 


magnet and breaks the field current of the small 
generator when the magnet current reaches a 
certain value. 


(3) The Measurement and the Indication 
of the Field 


A magnetic balance which is sensitive to about 
1 in 10‘ has been devised to measure the magnetic 
t Since the momentum of a given alpha-particle beam is 
constant, we have 
___ Hbp_ 1.0 10*x0.5 
6H= r 60 


"| am thankful to Mr. Thomas Coor for his help in 
building this relay system. 





=80 gauss. 
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field and has been described elsewhere.” Ordi- 
narily the field strength can be, however, more 
conveniently found from the values of Hp and p 
of some known alpha-ray beam. 

To indicate a small variation of the magnetic 
field, a non-linear indicating coil carrying a small 
constant current has been always used." It has 
been put in that part of the field which is not 
occupied by the deflection chamber. A deflection 
of 1 mm on a scale at 180 cm from the coil indi- 
cates a change of the field equal to about 2 
gauss at a total field strength of 12 kilogauss. 
It was intended at first to use this coil, instead 
of the galvanometer circuit, in conjunction with 
the photo-cell-2A3 circuit (cf. Fig. 3) to control 
directly the magnetic field. However, since the 
slow zero shift of the coil caused much trouble, 
this plan was not adopted. 


C. The Detection Instruments 


Three detection methods have been used de- 
pending on the different strengths of the sources 
or different intensities of the alpha-ray lines from 
a single source: 


(1) The Ordinary Photographic Method 


This method is used only when the source is 
strong, for it needs four or five thousand particles 
per square millimeter in order to make a line 
visible. The camera is made of aluminum so that 
scattering from its parts will be small, and can 
easily be slid on to the support, D. The back plate 
for holding the photographic plate in position is 
lined with a 5-mm thick lead sheet so as to absorb 
any gamma- or x-rays from outside. The emulsion 
layer of the plate is about in the focal plane and 
has the dimensions 6.5 X 19 cm. 

The Eastman “‘fine grain alpha-particle plates” 
have been used throughout the experiments. 
They are sensitive to alpha-particles but less 
sensitive to gamma-rays or light and give a very 
clear background. After being exposed to the 
alpha-rays, they are developed in D11 for 6 
minutes at 22°C. They are then dipped into 2.8 
percent acetic acid for about one minute in order 
to neutralize the alkaline D11 (or the background 
will be deep brown) and then fixed until clear 


4 W. Y. Chang and S. Rosenblum, Rev. Sci. Inst. 16, 
75 (1945). 
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(15-20 min.) in a fixing solution (four parts of 
hypo to one part of acid hardener). The emulsion 
layers usually break into pieces during exposure 
in vacuum; this can be cured at least partly by 
coating the edges with collodion solution. 


(2) The Counting Method 


This method is used when the source is weak, 
that is, when it takes too long a time to make an 
exposure in order to have a line visible. It is also 
a very convenient method when one wants to 
know instantly the homogeneity of an alpha- 
particle group. The counters used here have been 
specially developed for the purpose of investi- 
gating alpha-ray spectra and have been already 
described elsewhere."* It may consist of several 
fine wires stretched in front of a common brass 
plate, each wire then behaving as an alpha-ray 
counter. Hence in a single run several distribution 
curves for the same energy group can be obtained. 


(3) The Photographic Track Method 


Each counter of any type has a natural back- 
ground. If the source is so weak that the number 
of alpha-particles reaching the counter is below 
the background, the counting method is no 
longer applicable. In this case the photographic 
track method has been shown to be very satis- 
factory. Each alpha-particle produces a track in 
the emulsion layer of a plate. In order to get 
tracks of considerable length, the plate must be 
inclined to the incident beam. Eastman fine grain 
alpha-particle plates, which have been especially 
developed for this track method, 2 cm wide and 
19 cm long, have been used throughout all the 
experiments in this method. The width of a plate 
is inclined at an angle with the incident beam 
(the latter is actually in the form of a thin sheet, 
cf. Fig. 2), while its length is perpendicular to the 
beam, so that the whole plate is still more or less 
in the focal plane. The plate holder of the camera 
mentioned in (/) is so made that it can hold four 
plates at the same time side by side. Each of the 
two is inclined at an angle of 30°, while each of 
the other two is at 45° to the incident beam. 


(1945). Y. Chang and S. Rosenblum, Phys. Rev. 68, 222 

6 T, R. Wilkins, J. App. Phys. 11, 35 (1940); M. M. 
Shapiro, Rev. Mod. a 13, 58 (1941); Chadwick, May, 
Pickavance, and Powell, Proc. Roy. Soc. A183, 1 (1944); 
Champion and Powell, Proc. Roy. Soc. A183, 64 (1944). 


They are developed and fixed in exactly the same 
manner as in (/) except that they are developed 
in total darkness and for four minutes only in 
order to get clearer background. 

When a plate is examined with a microscope, 
it is illuminated with green light. The green light 
provides a good contrast between the tracks and 
the background and gives considerable comfort 
to the eyes. An intense beam of green light can be 
obtained from a mercury arc through an ordinary 
green glass filter. On the plate a number of 
equidistant lines (spacing about 0.9 cm) parallel 
to the main line of alpha-particles are drawn 
across the plate with a razor blade. These refer- 
ence lines are used for checking the measurement 
of the distance from the main line and hence the 
energy of the alpha-particle groups (cf. III B (3) 
and Fig. 7). If the plate being examined with a 
microscope is properly moved lengthwise and if 
the counting of the microscopic views along the 
plate is correct, then the number of views, in 
each of which the tracks are counted, must be 
the same in each interval between two successive 
lines. 

From the disposition of the photographic plate, 
it is clear that the tracks produced in the emul- 
sion layer caused by the alpha-particles, which 
have been coming directly from the source, must 
appear (1) parallel to one another, (2) perpen- 
dicular to the length of the plate, and (3) incident 
into the emulsion layer at about the same angles. 
It is this directional effect which enables us to 
discriminate between the alpha-particles coming 
directly from the source and the diffusely 
scattered alpha-particles. Certainly, particles 
scattered in their very early paths or coming from 
some contaminated parts near the source may 
also produce tracks approximately satisfying the 
above three conditions. However, experiments 
showed that the first effect is extremely small, 
while the contaminated parts (if any) can easily 
be cleaned or covered with thick aluminum foils. 
From our experience the number of irregular 
tracks, under favorable conditions, is negligibly 
small in the region near the main line and is only 
a small percentage of the corresponding regular 
tracks in the region further away. 

Therefore, the track method not only records 
every particle, but also gives information, respec- 
tively, from the orientation and appe&rance of 
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the tracks about the origin and nature of the 
particles. This is why it is very useful for count- 
ing particles in the presence of a relatively large 
general background (cf. the experiments below 
on the short range Po alpha-particles). 


Ill. EXPERIMENTS WITH Po ALPHA-PARTICLES 
A. The Main Group of Po Alpha-Particles 


The behavior of the above spectrograph has 
been investigated with Po alpha-particles by 
using the three methods of detection. Sources of 
different forms have been prepared, by spinning, 
in RaD solution, thin platinum wires, narrow 
strips of nickel foils, and the edges (about 0.2 mm 
wide) of nickel foils. All surfaces to be used for 
coating with the source have been well polished 
previously. In the latter case the two opposite 


_faces of the edge have been always well covered 


with two pieces of thin glass plates, so that during 
preparation Po is deposited only on the edge. 
Any one source can be covered with the source- 
slit so that only the central strip portion is used. 
It has been shown that a thin wire source does 
not give as sharp and narrow a line image as a 
narrow strip source or an edge source does. The 
third method of preparation gives a uniformly 
strong source and is the most economical one, for 
all of the Po is deposited only on the useful edge. 
In any case, the source must be fresh, uniformly 
thin, and show no traces of tarnish in order to get 
a sharp-narrow line image. 

The forms of the number-energy distribution 
curves have been determined, respectively, by 
the above three methods of detection and agree 
with one another fairly well. The width of the 
main line at half intensity is under favorable 
conditions less than } mm, i.e., about 0.01 Mev.’® 
According to the dimensions of the spectrograph, 
two lines separated by 0.01 Mev or less can thus 
be resolved. The intensities at 0.25 Mev or more 
below the energy of the main group are less than 
0.1 percent of the maximum intensity. This is 
about 100 times smaller than that found by Roy 
Ringo.‘ The distribution curves determined by 
the counting system have been published 
elsewhere." 

1* The relation between Hpand a, the energy of the alpha- 
particles in Mev, is given by Hp=1.436X 10°at without 


relativistic correction. Hence 23p= 1.436 10%3a/(Ha'), 
p being in cm and H in gauss. 
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The microscopic examination of the tracks 
distribution (as well as the microphotometric 
examination” of an ordinary photographic line or 
the distribution determined by the counters) does 
not reveal any fine structure of the main line. If 
there were fine structure of the main group, the 
separation should be less than 5 kev, for, as 
mentioned above, our spectrograph is able to 
resolve two lines separated by about 5 kev. 
Inconsistent splitting of the main line into two or 
three lines (separated by one or two millimeters, 
that is, equivalent to 0.02 or 0.04 Mev'*) was 
happening at first but was soon discovered to be 
caused by the improper adjustment of the 
automaton which controls the magnet current 
(cf. II (B)). In the following experiments, the 
form of the main line has been generally used, 
therefore, to obtain information about the be- 
havior of the magnet current during exposure of 
the plates. 


B. The Short Range Po Alpha-Particles 


By carefully studying the region of about 2 
Mev below the main line of Po, a series of weak 
groups of alpha-particles from Po have been 
found. For these experiments the sensitive track 
method of detection has been used in conjunction 
with the above large spectrograph. A preliminary 
report of these results was published recently.'* 
In the following we shall describe the experiments 
in some details: 


(1) Experiments with the Different Sources 


(t) The two and one-half months old source.—A 
source of about } mc was used at first. It was 
prepared by coating the edge of a nickel foil with 
Po as described in A. In order to have alpha- 
particles coming out only from the central strip 
portion, the source was covered with the source 
slit. The track method of detection as given in 
(II C.) was employed in recording the alpha- 
particles, the time of exposure being four hours. 
The plate was then examined with a microscope 
of 430 X (710 X for more highly populated plates). 
The number of alpha-particle tracks which 
satisfied the abpve three conditions was counted 
in a strip of 20 successive views across the plate. 


17] should like to thank Mr. J. R. Winckler for his 
help in some of these microphotometric measurements. 
18W. Y. Chang, Phys. Rev. 67, 267 (1945). 









NOV 18- 30° PLATE 
THE tAme SOURCE 


) 


T 


a 
NO OF 45 1N 10 VIEWS (DIA=26wm) 


No, OF 4 IN 10 Views (DIA-.a6c~, 
g 








sn 


me JAN I4- 30° PLATE 
Yarme SOURCE 


NO. OF 4s IN 20 VIEWS (DIA. =.26-men) 














©.76 cm FROM MAIN GROUP 


10 20 = 10 20 Q. 10 
0.84 cm FROM MAIN GROUP O.1t em FROM MAIN GROUP 


DISTANCE FROM MAIN LINE (0-10*5.2 wm) 


(a) (b) 


(c) 


Fic. 4. The distribution of the first two groups of short range Po alpha-particles as obtained from three different 
sources. The ordinate represents number of alpha-particles in 10 or 20 views across the plate. The abscissa repre- 
sents distance from the main line in terms of microscope views along the plate, distance of the origin O from the 


main line being indicated. 


Under favorable conditions, the number of ir- 
regular tracks is, as mentioned before, negligible 
in the region near the main line and is only 
a small percentage of the corresponding regular 
tracks in the region farther away. To save 
time, the counting was carried out only for 
every other strip across the plate (cf. Fig. 7). 
The number of tracks was then plotted against 
the distance from the main line. The number- 
energy distribution curve so obtained reveals 
distinctly twelve extremely weak groups in the 
low energy region from about 0.04 to about 
2 Mev below the main line. The first two groups 
stand out more clearly only on the curve obtained 
with the weakest source (cf. (3)), where the 
background near the main line is smaller. The 
number of tracks for each point on a curve varies 
from 50 to 250. The region lower than 2 Mev 
from the main line has not been examined, for in 
order to have this region on the plate, the main 
line, which is used to infer the behavior of the 
magnetic field, will be out of the plate. 

(1t) The condensed source-—The above source 
was already about two and one-half months old 
when it was used for these investigations. But it 
had been kept in vacuum, and according to 


Karlick and Rona,” Po nuclei diffuse extremely 
slowly into nickel. However, we were still afraid 
that some contamination of the surface or non- 
uniform distribution of Po might falsify the re- 
sults or at least might be responsible for the com- 
paratively large background. A fresh source was 
then prepared by vaporization and condensation 
of several old sources according to a procedure of 
Rona and Schmidt.?° Po vapor was condensed 
onto a well-polished palladium rod having a 
diameter of 1 mm and a length of 50 mm. The 
end of the rod further away from the old sources 
was hard-soldered to a long piece of copper wire, 
which served to cool the palladium rod so as to 
increase the efficiency of condensation. The 
palladium rod and the old sources were placed in 
a quartz tube which had a diameter of 1.5 mm 
and a length of about 250 mm. The old sources 
were carefully heated red hot with a very small 
oxygen flame when a stream of hydrogen gas 
was passed through the quartz tube. The rate of 
flow of the hydrogen gas through the tube was 
about 1 cc per sec. (higher than that recom- 


19 B, Karlick and E. Rona, Wien Ber. 143, 217 (1934); 
J. Schintlmeister, Wien Ber. 146, 389 (1937). 
20 E. Rona and Schmidt, Wien Ber. 137, 103 (1928). 
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Fic. 5. The distribution of the third to twelfth group of the short range Po alpha-particles. These were examined by a 
430X microscope, while those in Fig. 4 by a 710 microscope. The curve in Fig. 5 is one of five such distribution curves, 
which were obtained with different sources, different fields, and different arrangement of the chamber slits, etc. 


mended by Rona and Schmidt) in order to get a 
uniform distribution of Po along the long rod. 
The source prepared in this way had a strength 
of about 7s mc and looked shining and smooth. 
An exposure was made with this source also for a 
period of four hours. The energy distribution 
curve thus obtained is very similar to that 
obtained from the first one, except that the 
background is smaller. 

(tit) The fresh source——A stronger source of 
about } mc was freshly prepared by spinning the 
polished edge of a nickel foil in RaD solution. 
The distribution curve from this source is again 
similar to those above, and’ the background is 
smaller than that from the old source but larger 
than that from the condensed source. All the 
strengths given above are strengths actually 
used for exposures, 





(2) The General Results 


With these three sources, altogether five differ- 
ent exposures were made under different mag- 
netic fields. Three of them were done without the 
Lucite chamber slits and the other two with the 
chamber slits. In the former case three thin 
lead blocks lined with aluminum and having a 
height of 9 cm and a length of 18 cm were put 
into the chamber, so that they diverged out 
from the center of the line joining the defining 
slit and the nearer end of the photographic plate. 
After the first two exposures the material for 
making the defining and source slits was also 
changed from brass to aluminum; this resulted 
in a smaller background. The distribution curves 
obtained under these different experimental con- 
ditions are very similar to one another (except 
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Fic. 6. The last three groups of a curve like Fig. 5, obtained with the }-mc source. Inclination of the plate and its time 
of exposure are the same as in Fig. 5, but magnetic field is about 500 gauss more and other experimental conditions are 
different. After subtracting the general background, the integral intensities of these groups to that of the corresponding 
groups in Fig. 5 are about in the same ratio as the strengths of the two sources. 


the first curve, which was obtained by screening 
the photographic plate from stray light with one- 


cm air equivalent aluminum foil at a distance of 
about one centimeter away from the plate, so 
that the lines are broader). Plates inclined up- 
ward and downward with the incident beam 
gave the same results. The energy positions 
of the corresponding groups with respect to the 
main line in the different curves agree with one 
another to within about 0.04 Mev. The corre- 
sponding separations between these small groups 
themselves, in the different curves, however, 
agree with one another to within about 0.02 Mev 
only. This was chiefly caused by some errors in 
the measurement of the distance between the 
reference line and the main line so that in some 
curves the groups have been all shifted by the 
same amount from the main line (cf. Figs. 5 
and 6; also cf. (3) for discussion). The intensities 
of any one group in the different curves are, as 
will be also seen in the following, about in the 
same ratio as the relative strengths of the three 
sources (reduced to that actually used for ex- 
posure). The widths of these groups also increase 
with the width of the main line. 

Figure 4 shows the distribution of the first 
two groups. (a) was obtained with the 4-mc 


source having the plate inclined at 30° to the 
beam, (b) with the }-mc source (plate inclined 
at 30°), and (c) with the y)-mc source (plate 
inclined at 45°). They were examined with the 
710X microscope. It is seen that these two 
groups stand out more clearly in the curve from 
the s5-mc source, where the background near 
the main line is smaller. Figure 5 is one of the 
five curves mentioned above showing the dis- 
tribution of the groups from the third to twelfth. 
It was taken with the }-mc source having the 
photographic plate inclined at an angle 30° to 
the incident beam. In Fig. 6 are the last three 
groups of the a distribution curve, like Fig. 5, 
which was obtained with the 4-mc source, the 
inclination of the plate and the time of exposure 
being the same as those in Fig. 5. It is seen that 
the separations between these groups are prac- 
tically the same as those between the last three 
corresponding groups of Fig. 5, although they 
are all shifted toward the main line by a constant 
amount (cf. (3) for discussion). It is also interest- 
ing to note that, after the general background is 
subtracted, the ratios of the integral intensities 
of these groups to those of the corresponding 
groups in Fig. 5 are, respectively, 183:141 
~182:133-205:1494:4, which is about the 
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ratio of the strengths of the two sources. Similar 
relations have also been found for the other 
groups from these two sources and from the 
fy-mc source. This seems to indicate that these 
weak groups are related with the main group 
of the Po alpha-particles. 

The background of the curves must be chiefly 
caused by the straggling of the different groups 
of alpha-particles through the finite thickness of 
the source. By changing slightly the air pressure 
in the chamber or by performing experiments 
with and without the chamber slits (provided the 
top pole-piece was properly lined with alumi- 
num), the background was not affected very 
much. This shows that the effect of scattering 
from the residual gas as well as from the chamber 
was small. Moreover, the background due to 
radioactive contamination of the chamber as 
shown from blank experiments was found to be 
very small, about 2 percent of the third group 
intensity in the first part of the energy region 
mentioned above, and nearly zero in the latter 
half of the region. Certainly, the background 
(and also the number of irregular tracks which 
were not actually counted) became smaller, 
when the source- and the defining-slits were 
made of aluminum instead of brass, but the 
decrease was still too small to account for the 
total background. One experiment was also per- 
formed with a source in the chamber but covered 
with the shutter in order to test the proper 
operation of the shutter. The plate after being 
exposed for one hour in the chamber also showed 
a very small background. 


TABLE I. Energy of groups of Po alpha-particles. 


ALPHA-PARTICLES FROM PO 











Relative 
Group energy En =a0 —an integral 
Group in Mev in Mev intensity 
a 5.303 0 106 
a 5.113+0.02 0.190 (96) 
a2 5.065 +0.02 0.238 (126) 
a3 4.901+0.01 0.402 84 
a 4.838+0.01 0.465 49 
as 4.749+0.01 0.554 64 
ag 4.640+0.01 0.663 70 
ar 4.449+0.01 0.854 79 
ag 4.303+0.01 1.000 43 
ag 4.111+0.01 1.192 48 
10 4.016+0.01 1.287 19 
au 3.890+0.01 1.413 18 
a2 3.685+0.01 1.618 21 








(3)¥The Energy and the Intensity. 
Possible Error 


The energies of these groups have been calcu- 
lated by comparison with that of the main group, 
the latter being taken as 5.303 Mev. The value 
of p», the radius of curvature of the mth group, 
was obtained by measuring (1) the distance, /o, 
between the main group and the defining slit 
and (2) the separation, s, of this small mth 
group from the main line. s, was obtained from 
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Fic. 7. The scanning of a photographic plate by a micro- 
scope and the location of a small frou. relative to the 
source. The full circles represent the microscopic views 
in which the tracks are counted. Tracks in dotted circles 
are not counted. They are very much exaggerated. The 
lines 1, 2, --+ across the plate are the reference lines for 
the microscopic counting, etc. 


the distribution curve by calibrating the field of 
view of the microscope. The distance, k, of the 
source from the defining slit is 2.85 cm. These 
quantities are indicated in Fig. 7. The relativistic 
correction has been neglected. The energy, an, 
of the mth group is thus related to the energy, 
ao, of the main group by the following expression : 


(1) 
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2p0 PR? 

The counting of the tracks on the plate was 
usually started from some reference line near to 
the main line of alpha-particles ‘(cf. II C(3)). 
Therefore, s, equals the separation between the 
peak of the mth group and this reference line 
(from the distribution curve) plus the distance 
of this reference line from the main group. The 
latter quantity was measured up to the sharp 
edge at the high energy side of the main line. 
However, in the case where the main line is 
broad, this edge is not well defined, and this 
makes the measurement aflittle difficult. An 
error in the measurement of this quantity will 
cause a constant shift of all the groups in the 












distribution curve from the main line (cf. (2)). 
The energies so obtained are summarized in 
Table I. The energy value of each group in the 
table represents the average of three correspond- 
ing values from three different curves, which 
were obtained with the three Po sources. 

It may be useful to find how the errors in the 
measurements of J) and s, affect an, the energy 
of the mth group. This can be seen by differ- 
entiating a, with respect to J) and s,. Hence, 
letting po~/) and neglecting a term centaining 
1/po*, we have, 








~— (2) 


ban (= 1 Sndlo 5Sn, 
"Ape? 2po 


ao 2p0* 2po 


as the other two terms at most introduce errors, 
respectively, equal to 0.02és, percent and 0.0161 
percent (pp~60 cm, s,~15 cm maximum). 
Hence, the uncertainty of a, is chiefly deter- 
mined by the uncertainty in s,. The error, dsn, 
consists of (1) the error in measuring the distance 
of the reference line from the main line, (2) the 
error in locating the peak of the mth group due 
to different possible tracings of the distribution 
curve, and (3) the error in counting the number 
of the microscopic views along the plate (cf. 
Fig. 7 and II C. (3)). The first and the second 
errors are each about +0.5 mm. The number of 
views counted in each interval between two 
adjacent reference lines can fluctuate on the 
average by about +2 views (in about 30 views), 
i.e., about +0.6 mm. Hence, 


jam 3| 
2600 





x.303 }~ 0.008 Mev. 


In the last column of Table I are the relative 
integral intensities with respect to that of the 
main line. The area of each group in a distribu- 
tion curve was measured and then compared first 
with that of the third group (i.e., latter taken as 
unity), after the general background was sub- 
tracted accord{ng to its natural trend. Three 
series of such ratios were obtained from three of 
the distribution curves, each belonging to a 
different Po source. A series of average values 
was thus found. However, since the subtraction 
of the background is somewhat arbitrary and 
since the curves can be drawn differently, these 
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integral intensities can be off by a factor of 
magnitude ranging from about 0.75 to 1.5, 
especially in the earlier region where the arbij- 
trariness in determining the background is larger, 
It is to be noted that the integral intensities as 
well as the peak positions of the first two groups 
can hardly be accurately determined, because the 
background here is very large compared with the 
peak intensities of these two groups. The integral 
intensity of the main line cannot be determined 
from the same plate as the weak groups, because 
the population of tracks in this line is too high to 
permit counting. It was obtained under the same 
geometrical conditions by exposing a separate 
plate to the alpha-particles for a much shorter 
time (1 to 2 min.). The third group intensity was 
then compared to it. The intensity distribution of 
the main line was, as a check, also determined by 
an alpha-ray counter of thé type mentioned 
above and by using } of the 3-mc source (the 
other ? was covered with a metal foil). The ratio 
of the corresponding third group intensity to this 
main intensity was found to be about two times 
larger than the above ratio for which the track 
method was also employed for the main line. 
This is to be expected, because the number of 
alpha-particles counted by a counter must be 
smaller than that obtained by the track method, 
for this latter method, as mentioned above, 
registers every particle. Each value in the table 
represents the average of three corresponding 
values obtained from the three Po sources. The 
relative intensities reported before!® were relative 
peak values. ° 


(4) Comparison with Gamma-Ray Lines 


In 1935 Bothe* found five gamma-ray lines, 
namely, 0.202, (0.355), (0.433), 0.798, and 1.068 
Mev, the lines in the brackets being less certain. 
The emission of these weak groups of alpha- 
particles may lead to the formation of the Pb 
nucleus in different excited states as in the case of 
the ordinary fine-structure lines. If the effect of 
the recoil energy on the excitation energy is 
neglected, the’ differences between the main 
group energy and the individual group energies 
thus give the corresponding energy levels of 
the Pb™® nucleus.‘ These are shown in the 


21 Bothe, Zeits. f. Physik 96, 607 (1935). 
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Fic. 8. The distribution of Po alpha-particles in high energy region. The ordinate represents number of alphas 
in 20 views across the plate ; the figure 3.5 (opposite to 5) is the corresponding number of alphas when the plate is 
inclined at 30 degrees to the beam instead of 45. The abscissa represents distance from the main line in terms of 
microscope views along the plate; main line is about at the origin O and the curve extends to about 7.5 Mev. 


third column of Table I. When the Pb nucleus 
falls from an excited level to the ground level, 
a gamma-ray quantum is emitted. If this inter- 
pretation is correct, an exact correspondence 
is expected between the above gamma-ray lines 
and these alpha-particle groups. It is possible 
to find transitions in our level system which 
agree approximately with Bothe’s gamma-ray 
lines. But this comparison is rather arbitrary, 
as there are 78 possible transitions between the 13 
levels, and so more than one possible transition 
can be found to give the same gamma-ray line. 
A decision may possibly be made to find one (out 
of the two or more transitions) which actually 
corresponds to the gamma-ray line, if one knows 
the spins of the various levels and the selection 
rules for gamma-transitions (cf. Discussion be- 
low). The average number of quanta, as reported 
by Bothe, is about 7 in 10° of the alpha-particles. 
This ratio seems to be too small in comparison 
with the relative intensities of our alpha-particle 
groups. However, a smaller gamma-ray intensity 
relative to the main group intensity of alpha- 
particles is expected, as the method of detection 
he used for gamma-rays has an efficiency, which 
is neither high nor definitely known. In general 
we may say that, as far as energy and intensity 
are concerned, these gamma-ray lines are com- 
patible with the alpha-particle groups. 


(5) The High Energy Region of the Main Line 


An exposure of four hours to the alpha- 
particles from the strongest source was made 
with a higher magnetic field, so that the main 
line was near to the low energy end of the plate. 


Microscopic examination has been extended over 
the region from 0.01 to about 2.2 Mev above the 
main line. It shows a very small continuous 
background decreasing slowly to zero after about 
2 of the region mentioned above. Figure 8 shows 
such a distribution curve. It is seen that the 
distribution of the tracks varies only within the 
statistical fluctuation, and that it is only about 
1.5 percent of the third group intensity before it 
comes to zero. Therefore, if there is any group in 
this high energy region of the main line, its 
intensity must be smaller than 2 in 10° of the 
main group intensity. This small background was 
found to be partly because of contamination of 
and partly because of scattering from the 
chamber. 


IV. DISCUSSION ON THE SHORT-RANGE Po 
ALPHA-PARTICLES 
A. The Spin Changes from the 
Alpha-Decay Theory 


It has been mentioned above that to fit Bothe’s 
five gamma-ray lines into the level system of 
Pb?°* deduced from the alpha-particle groups is 
rather arbitrary, for there can be more than one 
possible transition to give the same gamma-ray 
line. However, if the spin changes, corresponding 
to the different alpha-transformations, could be 
estimated, the comparison might be much less 
ambiguous. The spin values of the various energy 
levels of the Pb nucleus may be deduced from the 
above spin changes. Hence one may then apply 
certain selection rules in conjunction with con- 
sideration of the intensities of the gamma-ray 
lines and the alpha-particle groups. A correct 
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transition corresponding to each of the five lines 
may thus be chosen. 

We have estimated” the spin changes by using 
the current alpha-decay theory of Gamow” and 
others: This theory is based on two main assump- 
tions, namely (1) the existence of the alpha- 
particles as subunits in a nucleus and (2) a 
simple rectangular hole for the auxiliary po- 
tential. The probability of penetration through 
the potential barrier is then calculated as a 
function of kinetic energy of the alpha-particle. 
The decay constant, Ax, corresponding to the 
a, transformation is connected with the velocity, 
vz, Of the a, group by the following equation: 


B 
logio Xe = A ——+CL (reer) }*, (3) 


U% 
where A, B, and C are constants for a given 
nucleus, and ), is given by 
M=ANi/Dd: Ni). (4) 
d is the total decay constant, and N;/>>; Ni is the 
intensity of the a, group relative to the total 
intensity, and is given in the last column of 
Table I. rere is defined as the effective radius of 
the residual nucleus and is given by 
(rett)x=Tol1—0.002j.(je+1) ], (S) 
where j; is the spin change in the a; transforma- 
tion, and ro is the true radius of the residual 
nucleus. fers is equal to ro when the change of 
nuclear angular momentum is zero. 


TABLE II. Spin changes in different Po 
alpha-transformations. 








Type of 
alpha- 


trans- 
formation As X10" sec.~ 


ao 5886 

a 0.5652 

a2 0.7419 

as 0.4943 
0.2887 
0.3769 
0.4120 
0.4651 
0.2531 
0.2826 
0.1119 
0.1059 
0.1236 


4) 
! ( 1 +753) 
x<10-° 
cm /sec. 


1.629 
1.600 
1.593 
1.566 
1.557 
1.542 
1.524 
1.492 
1.467 
1.435 
1.418 
1.395 
1.358 


Yett X10!2 
cm 
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#2, W. Y. Chang and Thomas Coor, unpublished. 


% G. Gamow, Constitution of Atomic Nuclei and Nuclear 
Transformations (Oxford University Press, 1937), second 
edition. 


Since A, B, and C are known and \, is.obtained 
from (4), rer is thus calculated from (3) using 
Table I. From (5) 7, can be then deduced. The 
various spin changes corresponding to the differ. 
ent alpha-transformations have been so esti- 
mated and are given in Table II. In order to 
obtain the value of ro, the value of jis corre- 
sponding to the maximum value of fer: (cf. 
Eq. (5)) has been arbitrarily set equal to zero. 

It is seen from the table that the alpha- 
transformations from the normal state of the 
ssPo”” nucleus to the excited states, lower than 
the seventh one, of the s2Pb?°* nucleus would 
require the spin changes of as high as 14 or 15 
units, while those to the states from the seventh 
to the eleventh result in changes of 11 to 7 units. 
These abnormally large spin changes are rather 
amazing, for in any transformation (particle or 
radiative) no such large spin changes have ever 
been observed. Moreover, since both the mass 
number and the charge number of either Po or 
Pb are even, the angular momenta of both nuclei 
in their normal states are expected to have zero 
values; accordingly the spin change in the ap 
transformation between the two normal states of 
the two nuclei should be also expected to be equal 
to zero. The value, 13, as calculated from the 
theory, for the ao transformation is far from the 
value, 0, expected. Certainly, we may assume 
that, since Po is an unstable nucleus, its normal 
state may have a spin as high as 13. But this 
assumption would also lead to the corresponding 
high spin values for the excited states of the Pb 
nucleus, which would then completely exclude 
the possibility of radiative transitions from high 
levels to low levels. This is in contradiction with 
Bothe’s observation of the gamma-rays from Po. 
Evidently, there are two possible alternative 
ways to account for this discrepancy; that is, 
either these discrete energy groups of the alpha- 
particles do not originate at the Po nucleus, or 
the alpha-decay theory needs some modification. 
However, considering the different experimental 
conditions described above, it is more probable 
that these discrete energy groups of the alpha- 
particles are from within the Po atoms than that 
they are caused by some external effects. Be- 
sides, as far as energy and intensity are concerned, 


_ these weak groups of the alpha-particles are 
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Fic. 9. The Geiger-Nuttall curves showing the deviation of the curve for the Po 


generally compatible with the gamma-ray lines 
from Po as determined by Bothe. 


B. The Geiger-Nuttall Relation 


The deviation of the experimental results from 
the theory can be more clearly revealed by a plot 
of the Geiger-Nuttall relation, i.e., logio A US. @n, @n 
being the energy of the alpha-particle in Mev. 
This is shown in Fig. 9. Curve I is the plot for 
these small groups of Po including the main 
group. Curve II is the one for the members of the 
radium family, while curve III is for the fine 
structure groups of ThC. It is supposed that 
curves II and III can be generally described by 
the above theory of alpha-decay. It is seen from 
curve I that logis \ drops very rapidly at first 
and then decreases very slowly with decreasing 
energy, a», of the alpha-particles. After this drop 
there is a large systematic departure from the 
curve for the Ra family. To be in line with the 
curve of the Ra family, a; to a3 on the one hand 
have too small intensities, while a7 to aj2 on the 
other hand have too large intensities. However, 
it may be interesting to mention again that, as 
reported in his papers, Bothe has not observed 
either any intense gamma-rays with intensities 


alpha-groups from that for the Ra family. 








which could be compatible with this high energy 
portion of the Ra family curve. 

Now let us consider the low energy portion of 
the curve. Take the eighth group, as, as an 
example. This group may correspond to Bothe’s 
1.068 Mev gamma-ray line. The observed partial 
decay constant of this group is about 3000 times 
larger than the value it should have if it would 
lie on the Ra family curve. In other words, if 
there should be any group of this energy, its 
intensity must be about 14 in 10° of the main 
group intensity of Po alpha-particles in order to 
be in line with the Geiger-Nuttall curve for the 
Ra family. Bothe’s gamma-ray intensity is about 
7 in 10° of the main intensity of Po alpha- 
particles. This is about 500 times larger than the 
above ratio required to fit in the Ra family 
curve. Therefore, comparing these experimental 
results with the Geiger-Nuttall curve for the 
Ra family, we are led to draw the following 
alternative conclusions: (1) Bothe, for some 
reason in his measurement, estimated the gamma- 
ray intensity 500 times too large (which seems 
rather unlikely, especially because the measured 
gamma-ray intensity is almost always smaller 
than its true value), and the discrete energy 
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groups of the alpha-particles do not originate in 
the nucleus of Po. (2) Both his gamma-ray lines 
and our discrete energy groups of alpha-particles 
do not have their nuclear origins in Pb and Po, 
respectively. Or (3) some modification is needed 
either in the development of the theory or in 
the application of the theory (cf. Professor 
Gamow’s suggestion in the next section). Since 
the intensities of the gamma-ray lines and these 
alpha-particle groups are of similar order of 
magnitude and since the gamma-ray energies are 
generally compatible with that of these alpha- 
particle groups, it is very likely that they are 
related with each other in the usual way. 


C. Some Possible Explanations 


A close examination of Eq. (3) immediately 
shows that the theoretical decay constant de- 
creases much more rapidly with decreasing 
velocity of the escaping alpha-particle than the 
observed decay constant. It seems, therefore, 
that there might be some other mechanism which 
would help the escape of the alpha-particle 
through the potential barrier, thus giving rise 
to an additional probability of decay to the 
probability of penetration through the statical 
potential barrier. It may be reasonable to as- 
sume that, during the alpha-disintegration, the 
nucleus might be vibrating, and this nuclear 
vibration might in some way or other affect the 
probability of penetration through the barrier. 
This was also generally suggested by Dr. J. M. 
Jauch** and seems to be fairly plausible. The 
idea of Goldhaber** that the Po nucleus may be 
in a state of pre-fission may be considered to be 
somehow connected to the idea of nuclear 
oscillation just mentioned. It will be interesting 
to determine if those easily fissionable nuclei 
have alpha-ray spectra similar to that of Po 
nucleus just discussed above. 

To fit the theory with the experiment, the 
term containing v, in Eq. (3) should be com- 
pensated somehow so that it would not change 
so rapidly with »,; otherwise, if it is left as it is, 
the value of res; and hence that of 7, have too 
large a variation as already shown above. We 
have tried to accomplish this by adding a term 
of the form D(ao—a,)* to the right-hand side of 
Eq. (3). D is a constant of about 3 if the alpha- 

* Jauch, Goldhaber, and Gamow, private communica- 


tions. I am grateful to them for their interest and dis- 
cussion about these difficulties. 
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energies ao, a, are expressed in Mev,:and a 
value of N~2 would suit the purpose roughly, 
For those groups which are near to the main line 
(e.g., the ordinary fine structure of alpha-par- 
ticles), this term is negligibly small, and hence 
the original form of the theory is still applicable. 
For those groups farther away from the main 
group, it becomes more important. Physically, 
since (a9—a,) represents the energy of excitation 
of the product nucleus, the above term so pro- 
posed must represent the contribution to the 
decay constant, 4, of some mechanism which 
plays a much more important role when the 
nucleus is being highly excited. If the ‘‘nuclear 
oscillation” during alpha-disintegration would at 
all contribute an additional probability of decay 
to that of penetration through the potential 
barrier, one might expect this additional proba- 
bility to be larger as the nucleus is being more 
highly excited. Therefore, it may be feasible that 
the above-proposed term may be connected in 
some way or other with the “nuclear oscillation.” 
However, the actual theoretical ground is com- 
pletely lacking at present. 

An alternative explanation as suggested by 
Dr. G. Gamow"™ is that each of the small 
observed groups might have “hyperfine struc- 
ture.” Therefore, if one could determine the 
intensities of these “hyperfine lines’ in each 
observed group and apply the above formula to 
these ‘“‘hyperfine lines’’ separately, the result 
might come out in agreement with the alpha- 
decay theory. It is interesting to note, as pointed 
out by Gamow, that when the formula of Hardy 
and Ramanujan* for the number of partitions 
of any integer n: 


1 
p(n) —— [w(2n/3)*] (6) 


is applied to our case (we assume basic levels of 
0.2 Mev so that n—~1.618/0.2-8 cf. Table I1), 
the calculated statistical weights, p(m), of higher 
levels (corresponding to the higher values of n 
up to 8) increase more or less in the same manner 
as the effective radii in Table II. This seems to 
suggest that the effective radii so calculated in 
Table II may perhaps correspond to the super- 
position of the ‘‘sub-partial” decay constants of 
the “hyperfine structure lines’’ in each observed 


% G. H. Hardy and S. Ramanujan, Proc. Lond. Math. 
Soc. (2) 42, 75 (1918). 
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group. The criterion in applying this equation to 
nuclear energy levels is that each energy of 
excitation may be regarded as to represent com- 
binations of a number of superimposed neutron- 
and proton-levels in the nucleus. It has been used 
by Bohr and Kalckar®* to calculate the density 
of nuclear energy levels for high excitation and 
has been shown to be satisfactory in explaining 
the experimental results of slow neutron capture. 
This formula is, of course, an asymptotic one, 
which holds better only for larger values of n. 
A much smaller basic level should be taken, in 
order to satisfy the above condition and hence to 
increase the calculated statistical weights so as 
to be more compatible with the observed in- 
tensities. However, the assumption of such a 
small basic level is subject to experimental 
verification. Unfortunately, it is not possible at 
present to measure the intensities of these ‘‘hy- 
perfine lines’’ of alpha-rays even if they would 
exist, for the straggling of the alpha-particles 
through the finite thickness of the radioactive 
source can easily cause an energy spread of as 
much as 0.05 Mev or even more. It would not 
be surprising if one or more cases similar to that 
of Po could be established in the near future. 

Another possible explanation of these short- 
range alpha-particles, as also suggested by Dr. 
M. Goldhaber,”* is the possible alpha-branching 
of RaD or RaE which may be present in the 
Po sources as impurities. However, the following 
experiment does not seem to support this expla- 
nation: 

The relative beta-activities of two different Po 
sources, which were also used for the work on 
the short-range alpha-particles, have been deter- 
mined. These beta-activities are presumably due 
to the impurities of RaD and RaE. After sub- 
tracting the background (about 15 per minute 
~S5 percent of the total counts) obtained when 
the sources were covered by a copper plate 
2 mm: thick, we have obtained the ratio: 


(11/I2)beta= 1.4 approx. (7) 

The ratio of the main alpha-activities of the two 

sources (after being reduced to that actually 
used for the exposures) is about 

(11/2) aipna = 0.75. (8) 

Since the intensities of any one of the weak 


26N. Bohr and Kalckar, Kgl. Danske Acad. Bind 14 
(1936-1937). 
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groups obtained from the two Po sources are 
about in the same ratio as the main alpha- 
activities of the two sources, i.e., about 0.75, it 
does not seem that these weak groups of the 
alpha-particles are related with the beta-par- 
ticles, or more specifically, with the amounts of 
RaD and RaE present in the Po sources. If they 
were caused by alpha-branching of RaD or RaE, 
the intensities of any one group from the two 
sources should be expected to have the ratio of 
the beta-activities, 1.4. One would come to the 
same conclusion with the following considera- 
tion: Counting all these weak groups together, 
their total intensity is about 1 in 1000 of the 
main Po alpha-particles. Now even if we assume 
RaD and RaE impurities in the Po sources as 
much as 10 percent (in terms of radioactivity) of 
the Po alpha-activity, we can only lower the 
branching ratio to 1 percent. Such a high 
branching ratio should have been detected be- 
fore. Besides, since the lifetime of RaD (RaE in 
equilibrium) is about 40 times longer than that 
of Po, 10 percent “activity-impurity” would 
mean that the actual amount of RaD present 
by weight in the Po source should be about 4 
times that of Po. This is very unlikely. In fact 
from the beta- and alpha-activities it has been 
roughly estimated that the “‘beta-activity im- 
purity” in the Po source cannot be greater than 
1 percent. 

The idea of constructing an alpha-ray spectro- 
graph by using the cyclotron magnet was ini- 
tiated by Dr. S. Rosenblum, who worked in 
this Laboratory for some months before I joined 
him and who left here in the spring of 1944 
for war work. Therefore, I am much indebted 
to him for the design of the Lucite chamber 
and for valuable information. I am also very 
grateful to Professor R. Ladenburg for his con- 
tinued interest and support throughout the 
progress of the work. I also wish to express 
my thanks to the former cyclotron group of 
Princeton for putting some of their apparatus 
at my disposal, to Dr. M. Blau for preparing 
some of the sources, to Mr. T. Coor for helping 
in some of the experiments, to my wife, C. S. 
Wang Chang, for her discussion and help in 
presenting this paper, and to Dr. E. C. Campbell 
for reading the manuscript. Thanks are also due 
Mr. C. L. Brecht for the machine work in con- 
struction of the deflection chamber. 
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The present paper gives the results of an analytical study of the directional intensity anomalies 
measured by Ribner and Cooper in Missouri in 1939-1940. The method used followsdirectly from 
the log vs. log graph found useful by previous workers. The straight line of Fig. 1 is then repre- 
sented by a least-squares equaiion of the form, 7 =A sec’ ¢. The anomalous points are shown to 
fall below the line and to be representable by a second line parallel to the first. The validity of the 
results is verified by direct comparison with the probable errors of measurement, the best test 
being Deming’s calculation of chi square. The use of this method brings to light the following 
results. The data of the two observers show an excellent agreement. The least-squares constants 
of the calculated equations show a real but rather irregular azimuthal variation. Special con- 
sideration of the reported time fluctuations of these fine-structure anomalies shows that in 
practically no case are they distinguishable from statistical fluctuations in counting. The 
anomaly close to the zenith is shown to be caused by magnetic effects, and not atmospheric ab- 
sorption. In the case of the other two anomalies, the importance of the secondary formation 
process is stressed. It must lead to an energy ratio S between the primary and mean secondary 
energies, the latter being concerned in the absorption. Comparison with the atom-annihilation 
results for the soft radiation suggests that the absorption anomalies, with S=2.2, may show 
agreement, but the magnetic angmaly, with S21, cannot agree. There is a slight suggestion 
here of two kinds of primaries, both involved in mesotron production. 





EVERAL observations have been made in 

recent years of the directional variation of 
cosmic-ray intensities at stations in Missouri and 
Mexico City.'~5 They afford an excellent set of 
data for use in an analytical representation if 
such can be found. The key to the form of 
representation is indeed available in earlier, less 
detailed measurements of the directional in- 
tensity variation. The present paper, Part I, 
presents the treatment of the data of Cooper! 
and Ribner* in Missouri. Part II, on Schremp and 
Banos’s Mexico City measurements,’ is nearly 
ready. 

The analytical form to be used to represent the 
fine structure may be inferred as follows. At sea 
level, the intensity-zenith angle curve has long 
been represented by the square of the cosine of 
the zenith angle.* Even for the fine structure,'? 
an M. Cooper, Phys. Rev. 55, 1272L (1939); 58, 288 
aad S. Ribner, Phys. Rev. 55, 1271L (1939); 56, 1069 
é' Ey. Schremp and A. Banos, Phys. Rev. 58, 662L 


(1940); 59, 614L (1941). 
*M. L. Yeater, Ph.D. Dissertation, Washington Uni- 


vers! 
WL. Yeater, Phys. Rev. 67, 74 (1945). 
*T. H. Johnson, Rev. Mod. Phys. 10, 193 (1938). 





this represents a helpful first approximation, 
from which anomalous departures can be calcu- 
lated. Though this is not even approximately 
true except at sea level, it does suggest the use of 
a double logarithmic plot, log intensity vs. log cos ¢ 
(¢ is the zenith angle), to represent the findings. 
A similar indication is found in Johnson’s plot? of 
intensity as a function of atmospheric depth. 
One of the best examples of all is in the work of 
Wilson,® whose logarithmic plot (Fig. 5) is much 
more readily interpreted than the semilogarithmic 
Fig. 4 of his paper. Still another excellent ex- 
ample of the use of the logarithmic graph is 
Greisen’s Fig. 2,8* in which it is used to interpret 
results from a considerable range of altitude and 
zenith angle. 

The first step in the present method therefore 
consists of the plotting of log j as a function of 
log cos ¢. The result is found to be represented by 
a straight line except for the fine structure, which 
naturally represents departures from the smooth 


7 Reference 6, Figs. 22, 23. Although the graph repre- 
sents total intensity, a simple application of the Gross 
transformation (ibid. p. 197) shows that the same power 
law applies to the vertical intensity. 

8 V. C. Wilson, Phys. Rev. 53, 337 (1938). 

8¢ K. Greisen, Phys. Rev. 61, 212 (1942). 
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line. Therefore the second step consists of the 
determination of the equation of the straight line, 
together with some means of allowing for the fine 
structure. In the Missouri data, the fine structure 
is found to be representable as a lower, parallel 
straight line, whose equation is also determined. 
In this way, a semi-analytical representation of 
the fine structure is obtained. 

Considerable care was also taken in testing the 
fit of the calculated curve. The best fit was shown 
by Gauss® to correspond to the smallest value of 
2, the mean square departure from the assumed 
“true” curve. In this case, where the accuracy of 
the individual measurements is known, the use of 
Q alone does not permit a close comparison with 
the expected departures, which are all of different 
magnitudes. This must be done in a different 
way, comparing the individual residuals with the 
corresponding probable errors or standard errors. 
This is expected to give a normal distribution of 
V;/o; with a standard deviation (n—m)*/n}, 
where n —m is the number of degrees of freedom 
and m is the number of observation points. 
Though this comparison cannot be made with a 
high degree of certainty, it does show that the 
sum of the squares of these relative residuals will 
have the same distribution as X?, a fact pointed 
out by Deming.'® This was found to be the most 
sensitive test applicable in the present case, 


probably because of the large number of degrees 


of freedom represented. The probability corre- 
sponding to the calculated value of X? (see ap- 
pendix) is taken from a table,*" and the result 
interpreted as the probability that the calculated 
curve agrees with the measurements within the 
experimentally determined uncertainties. 


APPLICATION TO COOPER’S DATA 


The data obtained by Cooper! in Columbia, 
Missouri are particularly complete and suitable 
for application of the method outlined above. 
The cycles used were such as to make the 
observations at each azimuth a complete, reliable 
unit, and the uncertainties in each measurement 

*A. G. Worthing and J. Geffner, Treatment of Experi- 
mental Data (John Wiley & Sons, Inc., New York, 1943), 
cf. appendix below for notation and formulas. 

10 W. E. Deming, Statistical Adjustment of Data (John 
Wiley & Sons, Inc., New York, 1943). 

"4G. U. Yule and M. G. Kendall, Introduction to the 


Theory of Statistics (Charles Griffin & Company, Ltd., 
London, 1937), eleventh edition. 











are given. The logarithmic graphs of his results 
in eight chief azimuths are given in the lower part 
of Fig. 1. The general linearity of the curves, and 
the suitability of the logarithmic representation 
are apparent. 

In the attempt to fit the different azimuths and 
fine-structure irregularities with an equation, 11 
different least-squares approximations were tried. 
The important results of these different treat- 
ments are summarized in Table I. Where the 
agreement is poor, the values of 2 are sufficient to 
determine the adequacy of the treatment; in the 
closer cases the more complete comparison with 
the errors of the measurements is necessary. Ap- 
proximations C, D, IJ, and K do not include all 
the points, but were found useful in interpreting 
the data. 

Of these, approximations C and D were used in 
studying the fine-structure irregularities. Careful 
scrutiny of the points plotted in Fig. 1 gives the 
impression that the straight, uniform line passes 
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Fic. 1. Logarithmic graph of the cosmic-ray intensities 
measured by Ribner (top) and Cooper. The ordinates are 
the logarithms of the relative intensities. The abscissae are 
the negative logarithms of the cosines of the zenith angles, 
which differ by a constant term from the logarithms of 
atmospheric path lengths. 
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TaBLe I. Approximation summary—Cooper. Success 
obtained in fitting Cooper's data with different approxi- 
mations. The names and constants characterizing the 
different analytical treatments are given. 2 is the mean- 
square deviation from the “true” curve of the assumed 
form, a smaller value indicates closer agreement. P(>X*) 


is the probability, based on the known probable errors, ° 


of obtaining a poorer agreement with the true curve than 
that observed. The expected value is 0.5. P(S.D.) is the 
probability that the iecteen contributing to X? have 
the standard deviation theoretically expected from the 
true curve. 











Approx. Name Constants Q P(>xX?*) P(S.D.) 
A Single line a,b 135 
for alla 
B Single line a(a), b(a) 104 
a Line high ay(a), b(a) 41.0 
D Line low ay(a), b(a) 69.3 
E Double line aHn,az,b 102 0.000 
for all a 
F Double line aH(a), az(a), 62.9 0.037 0.040 
Uniform points b(a) 
G Double line an(a), az(a), 29.4 0.999 0.0002 
Individual points b(a) 
H Non-|| lines an (a), by(a), 66.6 0.025 0.014 
Uniform points ar(a), bt (a) 
I Top line cut like F >60.9 
short 
J 5° Points like F $2.2 0.37 0.71 
asymmetrical 
K F, exc. omit like F 46.2 0.55 0.75 
35°E and 25°NE 








through the high points of the anomalies. Since 
Cooper speaks of the high points as the anomalous 
ones, this will bear closer examination. The 
method used is the comparison of approximations 
C (line high, i.e., through the high intensity 
points and the high ¢ points) and D (line low, i.e., 
through the low intensity points and the high ¢ 
points). The calculation verifies the visual con- 
clusion that the “line high” of approximation C 
fits onto the high ¢ points better than the “line 
low” of approximation D, since Q is 40 percent 
lower. It is for this reason that in Fig. 1 and in all 
later approximations the high line is extended 
through all the range. Indeed, it seemed better to 
draw only the one line, and let the anomalously 
low points of the fine structure be visible as 
separate points. : 
Approximations C and D were also used to test 
the reality of the fine-structure anomalies, as 
indicated in Table II. This was carried out by 
getting the weighted mean departure of the low 
points from the high line, and the high points 
from the low line, and comparing with the 


DANA T. 





WARREN 






standard error, o, of that departure. The proba- 
bility of obtaining such a departure, if we assume 
a normal distribution whose true value is zero, is 
also given. It is seen that all the departures fall 
below the 0.01 limit of probability except for 30° 
below the low line. This not only verifies Cooper's 
own conclusions that the fine structure is real and 
not the result of statistical fluctuations, but goes 
much further in showing how it may be expressed. 
From approximations C and D we may draw the 
following conclusions: (1) The fine structure is 
real, and the points cannot all be represented by 
one line ; (2) The points at 5°, 10°, and 20° belong 
on a high line; (3) The points at 0°, 15°, 25°, and 
30° belong on a low line; (4) Better results will be 
obtained by fitting the points at 35°, 40°, and 45° 
on the high line. Thus the smooth line is found to 
be that of the high points, with the anomalous 
fine structure consisting of the points that fall 
below it. 

Certain other important results can be seen 
from the comparison of approximations A, B, E, 
and F. A and B differ from E and F, respectively, 
by including only one line instead of both a high 
and low line. As would be expected from the dis- 
cussion of C and D, the departures, as measured 
by &, are considerably greater in A and B than 
when the line is double (Table I). Also, A and E 
are like B and F, respectively, as to nature of 
hi 


{ nigh lines—Cooper. Test of 


reality of fine structure. The weighted mean departures 
(V) of the measured points from the indicated line are 
compared with their probable errors (¢). Assuming a 
normal distribution of departures the probability of ob- 
taining such a large value of V if the “true” value is zero 
is P(O) given in the fourth column. The calculations are 
made by the formulas: 


The value of P(0) was determined from tables of the 
normal curve. 


TaBLe II. Summary of 











is 108 P(0) 
Line high 
0° — 8.1 2.9 0.0039 
15° — 13.3 2.8 0.00000 
25° — 13.3 3.1 0.00002 
30° — 7.7 2.9 0.008 
Line low 
Oe 11.3 2.8 0.00008 
10° 14.4 3.1 0.00001 
20° 8.3 3.1 0.007 
30° — §.2 3.2 0.093 
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curve, but in the latter the curves are fitted 
separately in each azimuth. Since this gives again 
a considerably better fit, we may conclude that 
the intensities are actually different in the 
different azimuths, a. This was also born out by 
a detailed calculation of the standard error of the 
constants of the equations for one azimuth in 
approximation F. Since the standard error was 
about one-fifth of the range of the constants be- 
tween azimuths, it is pretty certain that the 
constants of the curves are different in the 
different azimuths. 

Approximations F and H may be compared to 
determine whether the low points form a line 
which is actually parallel to the upper, regular 
line, or whether they constitute a line of a 
different direction. In this case, the indications 
are not so certain. The tests of Table I agree in 
indicating that F is slightly better, but the differ- 
ence is not very marked. The differences in the 
calculated constants for F and H are quite small, 
and there is no indication that the slope of the 
low line is consistently greater or less than that 
of the high line, since the means for the eight 
azimuths in H are identical with those for F, and 
also with E. It seems quite possible that it is 
largely accidental that the anomalously low 
points appear to fall at approximately the same 
distance below the high, regular line. 

In drawing the preliminary lines corresponding 
to Fig. 1, it was found that the west and south- 
east azimuths give an appearance of a break at 
the 35° point, with the line more nearly horizontal 
at higher zenith angles. To test the reality of this, 
the partial approximation J was carried out on 
these azimuths, by stopping the line at 35°. In 
one case a slight improvement, in the other a 
slighter worsening was found. It was concluded 
that the over-all agreement could not be im- 
proved in this fashion. 

So far, then, approximation F looks the best. 
It does not appear that any azimuthally uniform 
change in treatment could give an appreciable 
improvement. There remains the question of 
whether it is good enough, i.e., whether the 
residuals obtained in F agree with the values of 
the residuals expected from the known probable 
errors. This question is best treated. by con- 
sidering the probabilities of Table I and Table III. 
It is seen that the total probabilities fall in the 


COSMIC-RAY FINE STRUCTURE 








81 


“doubtful,” or “significant but not highly sig- 
nificant” range between 0.05 and 0.01, the most 
commonly used thresholds for determining the 
reality of an observed effect. From these values 
alone, we would be tempted to consider approxi- 
mation F as “probably not good enough,” i.e., 
the probability of obtaining the observed resid- 
uals by the mere random operation of the errors 
of determination of intensities is rather low. 


TABLE III. Calculated values—Cooper. The values of 
the parameters calculated by least squares for the three 
best approximations are listed in columns 3-6. The last 
three columns give the values of the quantities used to 
test the success of the approximations. The notation is 
explained in the appendix. Note that the values of ay and 
a, are the calculated values at 25°. 








X? P(>X®*) 





Approx. a —an —a@, an-ay —b Q 
F E 0.086 0.106 0.020 2.18 129.4 19.33 0.007 
SE .089 .106 017 2.40 69.4 8.76 0.27 
.089 093 .004 2.21 5.7 20 1.00 
SW .082 101 019 2.30 32.9 2.91 0.89 
W .084 .095 O11 2.23 30.1 2.65 0.91 
NW .090 .100 .010 2.36 §2.3 8.32 0.31 
.095 .107 012 2.43 48.7 9.49 0.16 
NE .096 .102 006 246 134.3 23.61 0.002 
Total .711 810 .099 18.57 502.8 75.97 0.037 
Mean .089 101 012 2.32 62.9 
s +0047 .0037 .0054 10 
G I 0.088 0.110 0.022 2.14 78.9 8.17 0.32 
089 .107 018 2.37 40.6 5.59 0.54 
.089 .094 .005 2.21 4.1 65 1.00 
SW .082 101 019 2.30 32.9 2.91 0.89 
V 084 .096 012 2.20 17.3 2.03 0.96 
NW .087 .103 016 2.42 19.1 2.50 0.92 
094 .110 016 2.45 28.7 4.50 0.72 
NE -085 .109 .024 2.58 13.6 1.37 0.98 
Total .698 .830 132 18.67 235.2 27.72 .999 
Mean .087 .104 017 2.33 29.4 
s .0035 .0059 0058 14 
J E 0.088 0.110 0.022 2.14 78.9 8.17 0.32 
.089 .106 017 2.40 69.4 8.76 0.27 
.089 .093 004 2.21 5.7 0.90 1.00 
SW -082 .101 019 2.30 32.9 2.91 0.89 
W .084 .095 O11 2.23 30.1 2.65 0.91 
NW .087 101 014 2.40 26.1 3.69 O81 
.095 .107 012 2.43 48.7 9.49 0.16 
NE .093 .103 010 2.49 125.6 22.51 0.002 
Total .707 816 109 18.60 417.4 59.08 0.37 
Mean .088 .102 014 2.33 52.2 
s .0040 .0055 0053 0.12 
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For this reason we consider approximation G. 
No uniform treatment of the lines in all azimuths 
can give improvement over F. But the fact is that 
neither Fig. 1 nor Cooper’s data show the fine 
structure as uniformly symmetrical. The high 
and low points do differ in the different azimuths, 
no two agreeing exactly, and only the southwest 
showing just the distribution determined from 
approximations C and D. We may, therefore, try 
approximation G, in which the form of the equa- 
tion is as in F, but the distribution of high and 
low points, instead of being uniformly that 
determined from approximations C and D (“uni- 
form points’), is arbitrarily made to coincide 
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with the observed distribution in each azimuth 
(“individual points”). This will necessarily de- 
crease the residuals, without introducing any 
change in the number of degrees of freedom. This 
means that 2 must necessarily decrease, and 
therefore cannot be used as a test of the suita- 
bility of the treatment. Instead, any justification 
for this arbitrary assignment must be sought in 
the necessity of fitting the observations within 
the limits expected from the known probable 
errors. 

Since great theoretical interest attaches to just 
this question of whether the fine structure is 
actually symmetrical or not, special effort was 
made in the comparison of approximations F and 
G. The results for these two approximations are 
given in Table III. Briefly, they may be summed 
up in the statement that where approximation F 
is perhaps not good enough, approximation G is 
highly significantly too good. For instance, we 
may consider the over-all probability from the 
chi-test. This gives us only 0.037 for F, but 0.999 
for G. Since the probability given is that of 
obtaining a value greater than that observed, we 
can say that we are 99.9 percent certain to get a 
worse agreement with the “true” calculated 
equation than is actually observed for approxi- 
mation G. Again, the standard deviation of the 
residuals is so small that there is only a 0.0002 
probability of obtaining such a close agreement. 

It is pointed out below, in the discussion of 
these anomalies, that there is a strong theoretical 
reason to believe that they are actually azi- 
muthally exactly symmetrical. At least, this is 
true of the anomalies at 15° and at 25°-30° 
(naming them according to the anomalous low 
points, instead of by the more regular “promi- 
nences”’ used by previous writers), but it is not 
true of the anomaly at the zenith. Not only is 
there no theoretical reason to think that this 
must be azimuthally symmetrical, but in actual 
fact, the observations of Cooper,! Ribner,? and 
Yeater® all agree in showing that it is not sym- 
metrical. The difference between the east and 
west is the most marked, with none of the curves 
given showing a significant prominence (i.e., 
above the cos? ¢ curve) at or near 5° in the east, 
whereas they all do in the west. Except for this 
exact agreement, the difference might be con- 
sidered uncertain, and has, indeed, been over- 


looked by all previous writers. The observations 
of Yeater® show an east-west asymmetry of 3-5 
times its standard error (the uncertainty of the 
multiplying factor makes the exact value inde- 
terminate) at 2° east and west. It seems certain 
that the asymmetry of this anomaly may be 
taken as real. 

Therefore we consider finally approximation J. 
Cooper’s data show the 5° points low in the 
northeast and northwest azimuths. Unfortunately 
there are no other observations in these two 
directions, but the completeness of the verifica- 
tion of his results in the east and west may give us 
confidence that they are reliable. These two 5° 
points are therefore called low in approximation 
J, though they do not make a great deal of 
difference in the interpretation, as will be seen 
from a comparison of approximations F and J in 
Table III. Approximation J differs from F only 
in the east (zenith on the high line) and northeast 
and northwest (5° on the low line). It is seen from 
Table I that it shows a very good over-all 
agreement with the observations within the 
limits of the uncertainties of the individual 
measurements. However, Table III shows that 
even approximation J does not fit well in the 
northeast. Unfortunately, the irregularities (20° 
low, 25° high) of this curve cannot be checked 
against any other observations. The single point 
at 25° shows a 0.0004 probability of agreement in 
approximation J, so that even out of 80 total 
points the expectation of such a poor agreement 
anywhere is low. This presents some evidence of 
asymmetry in the outer anomalies as well. 

The situation, then, is this. We may insist on 
strict azimuthal symmetry, and select approxi- 
mation F. In this case, we must explain the low 
over-all agreement, and the low agreement in 
east and northeast. Or, secondly, we may 
abandon insistence on symmetry in the zenith 
anomaly and take approximation J. In this case, 
the only difficulty to be explained is in the 
northeast. Finally, we may abandon symmetry 
altogether, adopt approximation G, and attempt 
to explain the fact that the resulting agreement is 
much too good to be expected from the known 
statistical fluctuations. These alternatives will be 
considered in order of increasing probability. 

The greatest difficulty is found with approxi- 
mation G. There is, however, a possible explana- 














tion of the low values of X? in the time variations 
of these anomalies. These have been mentioned 
by all the writers on the fine structure, but it is 
shown below that there is very good reason to 
consider them spurious appearances from the 
random counting fluctuations. If, however, they 
could be supposed to increase the variance of the 
indiv:'dual measurements by 50 percent (an im- 
possibly high value), then the value of X? would 
be reduced to 3 of the expected value, as is seen 
from the formula for X? in the appendix. This is 
the value found in approximation G. This ex- 
planation also permits us to picture the instan- 
taneous fine structure as always symmetrical, 
though variations in the time fluctuations would 
make the results of different series average to 
asymmetrical totals. However, any simultaneous 
series should show strict symmetry. There is some 
evidence in Yeater’s observations‘ that this is not 
the case. In any event, the dubious nature of the 
time fluctuations and the great magnitude re- 
quired leave approximation G the weakest of the 
three. Sos 
Careful examination of approximation F shows 
that the difficulty with it is localized in the east 
and northeast azimuths, and indeed in two bad 
points on each azimuth. An attempt was there- 
fore made, in approximation K of Table I, to 
eliminate the difficulty by omitting two bad 
points. The attempt was very successful, X? 
having practically its most probable value. It 
seems very dubious, therefore, to rule out ap- 
proximation F on account of only two bad points 
out of 80, especially since the over-all probability 
is not decisively low anyway. 

However, these arguments apply still more 
strongly to approximation J. In this case, the 
omission of only one point in the northeast would 
give an even closer fit. And the over-all fit of J is 
already quite acceptable. There remains only the 
doubt of the poor fit in the northeast, which is 
certainly better than is the case with F. It may be 
ascribed either to fluctuations of the fine struc- 
ture or to some special disturbance of the 25° 
point. Altogether, the weight of the evidence 
seems to favor approximation J considerably 
over the others. 

Table III lists also the values of the constants 
calculated for the least squares fit to the two 
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parallel lines in approximations F, G, and J. The. 
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calculations were all made with uniform weights 
for the logarithms of the intensities, since trial 
calculations showed no appreciable improvement 
in the calculated value of X? (theoretically im- 
proved by using weights) when the considerably 
more laborious process was carried out. One 
complete treatment by the method of Deming'® 
was carried out for approximation F in the north 
azimuth, and the standard errors for the calcu- 
lated constants were determined from the recip- 
rocal matrix. The results are 


CaH = 0.0025, 
GaL= 0.0030, 
o, = 0.046. 


-* — 


Since the probable errors and the corresponding 
weights do not differ greatly from one azimuth to 
another, these values may be taken as reasonable 
estimates of the standard errors of the parameters 
listed in Table III. The values show considerable 
azimuthal variation. An attempt to account for 
these variations met with little success. Since the 
low line includes the origin ({=0, log 7=0), the 
value of az is largely determined by that of b. 
The most regular azimuthal variation is that of 
ay, whose Fourrier expansion shows a consider- 
able degree of convergence in approximations F 
and J, with a maximum of the first- harmonic in 
the southwest. The quantities 6 and ag—az 
showed poor convergence and not much agree- 
ment as to the phases in the different ap- 
proximations. 


APPLICATION TO RIBNER’S DATA 


The data of Ribner? are far less complete than 
those of Cooper,' and are not so directly suitable 
for the present form of analysis. The smaller 
number of cycles, and the necessity of superposing 
the results of different runs in each azimuth make 
the measured values subject to rather greater 
uncertainties than Cooper’s. The general agree- 
ment, however, is excellent, and Ribner’s values 
furnish very satisfactory corroboration of the 
deductions from Cooper’s data. The general 
nature of this agreement is shown in Fig. 1 and 
Table IV. The parallelism is evident, although 
there are four cases in which the difference is 
real. Three of these cases are covered by the 
rather large uncertainties in Ribner’s data, as 
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was found by fitting Ribner’s measurements to 
the equations determined from Cooper’s data. 
The exception is the value of a; in approximation 
F in the east. Such differences might be expected 
as a result of small differences in latitude, alti- 
tude, or transition and absorption in the counter 
train, none of whose effects on the fine structure 
have yet been investigated. There is also a slight 
indication of azimuthal asymmetry in the 15° and 
25°-30° lows when compared with Cooper’s 
approximation G, but these again are covered by 
the uncertainties in Ribner’s data. 

There is one very significant additional piece of 
data given by Ribner and not by Cooper. In 
addition to the probable errors calculated in the 
same manner as those quoted by Cooper, Ribner 
gives the values of the probable errors based on 
the number of counts alone. The errors by 
residuals include all variations, such as those 
caused by instrumental sensitivity or time fluctu- 
ations of the fine-structure pattern, as well as the 
statistical counting irregularities: The striking 
fact is that the errors by residuals are actually 
less than those caused by counting, though only 
by an amount neatly covered by the probable 
errors of the probable errors. However, there is 
here no evidence whatsoever of real fluctuations 
in sensitivity or in fine-structure pattern. This 
does not prove, of course, that very much higher 
counts could not reveal the reality of such 
fluctuations. But it does make it quite certain 
that such fluctuations cannot amount to 50 
percent of Cooper’s counting errors and so ex- 
plain the very high probability found in Cooper’s 
approximation G. It also makes it clear that the 
fluctuations in pattern reported by Ribner to 
Cooper cannot be distinguished from statistical 
fluctuations in counting. 


DISCUSSION 


Much of the treatment presented above is clear 
enough. Qualitatively, it agrees with previous 
treatments, except for the identification of the 
low points as anomalies, instead of the promi- 
nences. Quantitatively it goes farther and presents 
a semi-analytical representation of the entire fine 
structure. There are, however, a number of 
matters connected with the interpretation of the 
anomalies that will bear considerable elaboration. 


WARREN 


TaBLe IV. Summary—Ribner. Comparison between the 
results of Cooper’s and Ribner’s data. The columns com. 

re the qimined values with those listed in Table III. 

e significance of the column headings is the same as 
in Table III, and is given in full in the appendix. 














Approx. Name a —-@H —a@p @H-ay —b Q P(>x 
B Single line ‘ 
Ribner E 0.090 2.14 Iii 0.27 
Cooper E 0.095 2.11 210 
Ribner W_ 0.090 2.12 118 0.13 
Cooper W_ 0.089 2.19 59.0 
F Double line 
Uniform 
points 
Ribner E 0.085 0.097 0.012 2.18 82.3 0.35 
Cooper E 0.086 0.106 0.020 2.1 129.4 0.907 
Ribner W 0.083 0.096 0.013 2.15 27.7 0.70 
Cooper W 0.084 0.095 0.011 2.23 30.1 0.91 
G Double line 
Individual 
points 
Ribner E 0.087 0.104 0.017 2.1 53.8 0.51 
Coo: E 0.088 0.110 0.022 2.14 78.9 0.32 
Ribner W 0.084 0.099 0.015 2.14 16.8 0.90 
Cooper W 0.084 0.096 0.012 2.2 17.3 0.96 








Let us consider first the matter of time fluctu- 
ations of the fine structure. Ribner? mentioned 
such fluctuations in his observations, but it has 
just been pointed out that they cannot be dis- 
tinguished in his data from the expected varia- 
tions inherent in the counting process. Cooper! 
not only mentions the observation, but gives a 
detailed report of an extreme case occurring in 
the south azimuth, under the heading of “8 
anomalous cycles.” For these, he quotes the 
observations, together with their probable error 
determined from the number of counts. Two 
statistical tests of these eight cycles: have now 
been made. The first was the direct measure of 
the probability of observing such a series of eight 
cycles, by the chi-test comparison of the de- 
partures from the average curve with their 
errors. The result gave a probability of 0.21 or 
about $. Since there were 8} groups of eight cycles 
in the south, there is only a 0.14 chance that at 
least one such group would not show as much 
departure as this. In other words, the known 
limits of error prove that departures as great as 
these would be expected, though that is not to 
say that the next such irregularity would be again 
of the nature of a reversal of observed structure. 

This group of eight anomalous cycles was also 
investigated for its effect on the size of the prob- 
able errors. This can only be done roughly, with 
the assumption that the other cycles had total 
counts the same as those where the effect of the 
counting is known. But the result was entirely 

















negative. The probable errors listed for the total 
of 68 cycles are not increased by the presence of 
this apparently extreme case of fluctuation. 
Even if this does not prove the non-existence of 
the time fluctuations, it certainly shows that 
their effects cannot amount to the 50 percent (in 
the variance) required to make approximation G 
agree with expectation. 

Again, Yeater,*® with apparatus and experi- 
mental method designed specifically to determine 
these time fluctuations, was unable to establish 
any such effect. He further points out that 
Schremp’s earlier apparatus” was also inadequate 
for the purpose. From the comparison of all these 
lines of evidence, then, one may conclude that 
most, and probably all, supposed cases of fluctu- 
ations can be better explained as caused by 
statistical counting variations. Certainly the 
effect cannot be large enough to explain approxi- 
mation G. 

Now wé may consider the probable interpre- 
tation of the fine structure. From the earliest days 
of its prediction and consideration" the symmetry 
it shows has been taken as the key to its interpre- 
tation. Irregularities caused by atmospheric ab- 
sorption should depend only on the length of 
atmospheric path traversed, and therefore be 
independent of azimuth. The effect of the earth’s 
magnetic field, though imperfectly known" at 
these high latitudes, is not expected to show any 
such symmetry. It may show an east-west 
symmetry if there are oppositely charged particles 
of identical total energy in the primaries, and 
except to the extreme north and east show little 
variation from northwest to southwest or north- 
east to southeast. 

The consequence of this reasoning is that these 
fine-structure irregularities, with their at least 
approximate azimuthal symmetry, have always 
been ascribed to the effects of atmospheric ab- 
sorption. To show such irregularities the absorbed 
particles must possess certain characteristics. In 
the first place, their absorption must be of the 
nature of a range phenomenon, their effects 
ceasing rather suddenly after a certain amount of 


2 E. J. Schremp, Phys. Rev. 57, 1061A (1940). 
4% E. J. Schremp and H. S. Ribner, Rev. Mod. Phys. 11, 
149 (1939), 
“4 E. J. Schremp, Phys. Rev. 54, 153, 158 (1938). 


COSMIC-RAY FINE STRUCTURE 





85 





absorption. The theoretical work'’ seems to show 
quite conclusively that this is characteristic of 
the penetrating, mesotron component of second- 
aries, but not of the soft electron showers. In the 
second place, the energies of these absorbed 
secondaries must fall into definite, well-separated 
groups, whose suddenly completed absorption is 
responsible for the sharp drop in intensity. It is 
often overlooked that this must be the charac- 
teristic of the absorbed secondaries, whose ener- 
gies cannot be expected to be nearly so sharply 
defined as the primaries from which they are 
derived. Yeater,‘ for instance, elaborates many 
possible forms of absorption curves, all based on 
the assumption that the secondaries have the 
same energies as the primaries. The strength of 
the present method is that it allows for the 
straggling in secondary energies by the empirical 
selection of a baseline determined by the actual 
energy distribution. 

There is an additional very important charac- 
teristic of anomalies caused by absorption. It is 
clear that any increase in absorbing path, i.e., in 
zenith angle, can only decrease the intensity, 
never increase it. This is observed in the case of 
the present anomalies at 15° and 25°-30°. But it 
is not true of the anomalies at the zenith. Both 
Fig. 1 and Yeater’s curve’ show a ring of higher 
intensity all around the zenith, except that in the 
east it first falls lower and then rises to the same 
height. Such an anomaly, whether azimuthally 
symmetrical or not, could not possibly be caused 
by atmospheric absorption, and there is therefore 
no reason to suppose that its approximate sym- 
metry must be exact. In view of the close 
agreement of all the observations that it is not 
symmetrical, it seems better to abandon the 
assumption that it is, as was done in approxi- 
mation J. It is therefore probably ascribable to 
magnetic effects. Such magnetic effects are also 
caused by a discontinuous primary energy spec- 
trum (or possibly by peculiarities of the penum- 
bra, which are not looked for at this latitude.") 

Weare left, then, with two edges of atmospheric 
absorption and one magnetic anomaly caused by 
bands of energy in the primary distribution. It is 
very desirable, therefore, to compare these bands 
with other work on the energy bands of the 


% Cf. e.g., B. Rossi and K. Greisen, Rev. Mod. Phys. 13, 
240 (1941). 
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primaries, together with their theory. The most 
extensive material of this nature is the atom- 
annihilation theory of Millikan, Neher, and 
Pickering.'® Although this theory has been tested 
chiefly on the soft, shower component, recent 
evidence!’ indicates that the mesotron com- 
ponent may show the same primary energy 
spectrum. The pertinent portion of the primary 
spectrum is the four lightest atom-annihilation 
bands, named by Millikan, Neher, and Pickering 
the helium (1.9 Bev), carbon (5.6 Bev), nitrogen 
(6.6 Bev), and oxygen (7.5 Bev) bands. Of these, 
the helium band lies well below the limit (2.5 
Bev) for observation at sea level. The other three 
bands may be compared roughly with the two 
atmospheric absorption anomalies at 15° and 
25°-30° zenith angle. Unfortunately, the shape of 
the anomalies is not very exactly known from 
these rather widely-spaced measurements. But it 
was found by repeated trials that if the significant 
zenith angle was taken at the succeeding high 
points, 20° and 35°, a fairly good agreement with 
the carbon and nitrogen bands could be obtained, 
together with a prediction of another high point 
at 44°. On Yeater’s more detailed curves, a 
symmetrical prominence is readily seen at 42°, 
though it naturally does not show up on the 
present curves which do not include measure- 
ments at this point. 

It seems therefore quite possible that the atom- 
annihilation bands can be used to account for the 
present atmospheric absorption fine structure. To 
do so requires the use of a secondary energy 
factor (cf. appendix) S=2.16. This compares 
favorably with the factor 2.20 deduced from the 
data of Gill'* on the assumption that the sea level 
cosmic-ray shelf edge is determined by the atom- 
annihilation value 5.5 Bev. 

In the case of the magnetic anomaly, the situa- 
tion is different. The energy determined in this 
case is the primary energy. At this latitude, 
magnetic anomalies at the zenith should be 
caused by primaries with rigidities* between the 
limits 0.215 stoermer (Stoermer cone) and 0.23 


R.A. Millikan, H. V. Neher, and W. H. Pickering, 
Phys. Rev. 61, 397 (1942); 63, 234 (1943); 66, 295 (1944). 

17 Es. M. Schein, W. P. Jesse, and E. O. Wollan, 
Phys. Rev. 59, 615L (1941); J. Hamilton, W. Heitler, and 
H. W. Peng, Phys. Rev. 64, 78 (1943). 

18 P, S. Gill, Phys. Rev. 55, 1151 (1939). 
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stoermer (Lemaitre and Vallarta full light cone’), 
The corresponding energies are slightly less than 
2.9 Bev and 3.2 Bev because of the weakness of 
the earth’s magnetic field on this side. Such a 
band of primaries could evidently only make its 
presence felt at sea level with a secondary pro- 
duction factor close to one, instead of over two. 
It appears that the intensity of this anomaly, 3 
percent in the vertical, may be so small when 
averaged over the entire sphere (a rough calcula- 
tion gives § percent) as not to affect Gill's 
observation'® of the cosmic-ray shelf. 

Such a primary energy band, however, does 
not fit into the atom-annihilation spectrum. Nor 
is the situation improved by the assumption that 
this anomaly is caused, like the others, by ab- 
sorption, for the energy difference from the next 
one is then far too small to match the spacing of 
the atom-annihilation energies. The discrepancy 
is further emphasized when it is noted that the 
present observations lie well within the plateau” 
observed between St. George and Pocatello in 
vertically incident soft radiation. 

Similar difficulties beset the attempt to account 
for the magnetic anomalies observed by Yeater‘ 
at 20° east and west of the zenith. He ascribes to 
them a rigidity of 0.20 stoermer, but such par- 
ticles could not give rise to secondaries capable of 
penetrating that thickness of air. It seems more 
probable that the fluctuations observed are in the 
edge of this same zenith anomaly, either as the 
fluctuation of the main cone to the west (for 
positive particles) or the Stoermer cone to the 
east. The interpretation given by Yeater, of the 
Stoermer cone in the west, leads to an unob- 
servably low energy. Whether the variations are 
to the east or west of the zenith (i.e., whether the 
energy is higher or lower than for the zenith 
position of the same cone) can probably be 
determined from: the fine structure of the east- 
west asymmetry, since the inclination of the 
magnetic meridian plane to the east will cause the 
edges to lie at different zenith angles, as was 
actually observed. 

The following tentative conclusions emerge. 
(1) The absorption edges may be obtained from 
atom-annihilation energy bands, with an S-factor 

19G, Lemaitre and M. S. Vallarta, Phys. Rev. 50, 493 


(1936). 
20 Reference 16, Fig. 7 on page 241. 
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of over two. (2) The magnetic anomaly, at the 
zenith, and with edges at 20° east and west, does 
not occur in atom-annihilation bands, but has an 
S-factor close to one. It seems possible there is 
here evidence for the dual nature of mesotron 
primaries, with only one kind able to generate 
soft electron showers. 


APPENDIX 
Notation and Formulas 


a=azimuth; ¢=zenith angle; 2=most probable mean 
square residual from “true” curve; 2=ZV;,*/(n—m), 
where V;=residual at point i, »=number of points fitted, 
and m=number of parameters adjusted from data. 
X?=2(V;/oi)*, where o;=standard error of measurement 
whose residual is V;. I.e., o;=(1/0.675) X probable error. 
When the measurement is a true frequency, this reduces 
to the more common definition X*= = V;,?/f., where f, = theo- 
retical frequency at point 7. P(>X*)=probability of ob- 
taining a value of X* greater than that observed. For the 
least-squares calculations by which lines were fitted to 
the logarithmic plot, the following formulas and definitions 
hold: j=relative intensity from data, x = log cos {>—log cos ¢ 
({9=25° in references 1 and 2). The equation fitted has the 
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form log j=a+bx, where a=intercept of calculated line at 
to; b=slope of logarithmic line = exponent of generalized hy- 
perbola whose equation is j = A (cos ¢)~*. When two lines are 
used to represent the fine structure, @y is for the higher line; 
a, is for the lower; oa#=(standard error in ag) = (prob- 
able error in ay)/0.6745; oa, =standard error in az; and 
o,=standard error in 6. In Table III, the standard devia- 
tions of these quantities are given for the azimuthal 


t 
variation s»)= {720 -5 * When a single constant term is 
Ux*Zy—Exdxy_, _ NExy—Zxdy 
N2x?—(2x)* ’) Nax?—(2x)?” 
When two constant terms are used, 


used, a= where y =log j. 





1 1 
—Yixti1y+ —Z2x2y —_ 21, 2xy 
ny, ne 





be 


J (Ein) + (Sax)* 2? 
a= Jey— (b/m) 2.x, 
my 


a2= + yy — (b/m2)Zox. 
Ne 


In these equations, the subscripts 1 and 2 distinguish the 
points assigned to the high and low lines, respectively. 
S=Epri/Eswc=ratio of primary energy to most probable 
energy of the secondaries it generates. 
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This paper discusses the possibility that radiation losses because of the high radial accelera- 


tions experienced by the electrons in an induction electron accelerator may introduce limi- 
tations in the design of accelerators for energies above 100 million electron volts. The effects 
of radiation losses on the electron orbits are calculated, and it is shown that not only should 
the orbit shift pulse necessary to bring electrons to a target inside the equilibrium orbit fall 
below the value expected in the absence of radiation, but also electrons should eventually 
arrive at the target with no orbit shift pulse whatever, at a phase of the field wave predictable 
from the theory. Both effects have been observed in the General Electric 100-Mev unit in a 
‘manner consistent with the predictions of the theory. The radiation itself has not yet been 


detected. 


1. INTRODUCTION* 


N the induction electron accelerator, the elec- 
trons are subjected continually to radial 
accelerations of the order of 10'’ meters per 


* Symbols:—Unrationalized m.k.s. units will be used 
throughout: The following symbols will be employed: 
A=peak value of applied magnetic flux density at the 
equilibrium orbit (webers per sq. m) 
A’ =peak value of magnetic flux in orbit shrinking pulse 
at the equilibrium orbit (webers per sq. m) 
Bo=applied magnetic flux density at the equilibrium 
orbit dodbers per sq. m) 





second per second. It has been pointed out by 


B, and B, are components of magnetic flux density 
(webers per sq. m) 
c=velocity of light = 3.00 x 10* m per sec. 
e=charge on the electron = 1.602 x 10~* Coulomb 
E, and E, are components of electric field (volts per m) 
fn and f; are normal and tangential components of the 
acceleration vector f (m per sec. per sec.) 
F(wt) = (wt/sin wt)—cos wt— (2/3) sin? wt cos wt 
h=Planck’s constant = 6.624 X 10-* joule sec. 
H, and H, are components of magnetic field 
I =beam current (amperes) 
mo=rest mass of the electron =9.107 X 10- kg 
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Iwanenko and Pomeranchuk' that radiation from 
these electrons will be appreciable and may set 
an upper limit to the energy attainable unless 
due attention is paid to new choices of operating 
parameters. The present paper discusses this 
conclusion and outlines the effects of radiation 
on the electron orbits. 

The discussion will be based on the General 
Electric 100-Mev unit which will be described 
in detail in a forthcoming paper in the Journal 
of Applied Physics.2 The radius of the ‘‘equi- 
librium orbit” in the G.E. accelerator is 0.833 
meter. Magnetic fields of the order of 0.4 weber 
per square meter are used to attain the highest 
electron energies. Focusing is achieved by shaping 
the poles of the magnet so that the field varies 
inversely with the ? power of the radius in the 
neighborhood of the equilibrium orbit. The target 
to be bombarded by the high energy electrons is 
located 0.07 m in from the equilibrium orbit. 
At the appropriate moment, the orbit diameter 
is shrunk by a short magnetic field pulse applied 
by an auxiliary system of coils so that the beam 
hits the target. 

If the electrons lose energy .continually by 
radiation, they will tend to spiral inward. The 
strength of the orbit shrinking pulse necessary 
to bring the electrons to the target will then be 
less than that expected from a calculation neg- 
lecting radiation effects. If sufficient energy has 
been lost by radiation, the beam may reach the 
target without any orbit shrinking pulse. Both 
effects have been observed in our accelerator. 

N=number of electrons cu. m in the beam 

n=index of variation of applied magnetic field with 
radius = — (3/4) in our case - 

Ro=radius of equilibrium orbit =0.833 m in our case 


U,=rate of radiation from single electron (watts) 
U>o=radiation energy density in the beam (joules per 


cu. m ‘ 
v=electron velocity (m per sec.) 
V=injector volts 
w=maximum width of beam in plane of orbit (m) 


W=electron energy (joules) 
€9= dielectric constant of free space = 1.11 x 10~!° farad 


per m 
wo= permeability of free space = 10~" henry per m 
p=charge density in beam (Coulombs per cu. m) 
¢= magnetic flux density (webers) 
w= 2x X60=377 
w’ = x/(duration of orbit shrinking pulse) 
(Note: 1 weber per square meter = 10,000 gauss.) 
— and Pomeranchuk, Phys. Rev. 65, 343 
1944). 
2 W. F. Westendorp and E. E. Charlton, ‘‘A 100-million 
volt induction electron accelerator,” J. App. Phys. 16, 
581 (1945). 


Distortions and time delays in the field and flux 
waves which might cause such effects have been 
investigated and found to be too low by an order 
of magnitude. Since no other explanation has 
been offered, the possible results of radiation 
have been calculated and seem to offer an ade- 
quate description of the observed phenomena, 
The nature of the radiation is discussed in 
Section 3 below. The radiation itself has not yet 
been detected experimentally. 


2. RADIATION FROM A SINGLE ELECTRON 


It is shown by Page and Adams that the rate 
of radiation per solid angle dw along the vector ¢ 
from an electron having velocity v and accelera- 
tion f is given by: 

10-"e? 
c-v\> 
arc'(1 —-— 


c? 


dU, 





X {(fX(c—v)) Xc}* watts. (1) 


If v is approximately equal to c and if f is normal 
to v, the radiation has a strong maximum in the 
direction of v. The half intensity contour will 
be a rough cone whose vertical angle is less than 
0.3° for energies of the order of 100 Mev. 

The total rate of radiation from the electron is’ 
given by integrating (1) to obtain 


2X10-e? fi 
_ 3c (1—v?/c?)8 
fn? 
pemccmmmnianens 
(1—v?/c?)? 





atts. (2) 


This is Page and Adams’ formula 75-6‘ expressed 
in m.k.s. units. As v approaches c, the ratio of the 
first term in Eq. (2) to the second term becomes 
so small that the first term is entirely negligible. 
The second term can be rewritten by use of the 
relations: 


Mofn _ 

(1—v2/c2)t 

fr=v"/r, (4) 

3L. Page and N. Adams, Electrodynamics (D. Van 


Nostrand and Company, New York, 1940), Chap. 7. 
* See reference 3, p. 328. 
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whence 


Mv 
————— =erB 
(1—v?/c*)# 


(5) 


We substitute from (4) and (5) for f, and 
(1—v?/c*) in (2) and obtain: 
2x 10-7e°B,*r? 
U,= watts. (6) 
3cmo* 





Equation (6) is essentially the relation used by 
Iwanenko and Pomeranchuk. If we insert the 
values of the fundamental constants in Eq. (6), 
we obtain: 

U,=5.44« 10-°B, 'r? watts. (7) 


The question may be raised as to the adequacy 
of special relativity in dealing with this problem. 
The equations of motion of a charged particle 
will, however, be found to be unchanged by any 
considerations introduced in the general theory 
of relativity.® 


3. NATURE OF THE RADIATION 


When the total radiation from a group of NV 
closely spaced electrons is evaluated by the 
techniques which led to Eq. (1), the net energy 
radiated proves to be proportional to N? since 
the fields are summed up and the energy is pro- 
portional to the square of the net field. On the 
other hand, the radiation fields due to the ele- 
ments of a continuous distribution of charge 
travelling in a beam are self-canceling, so that 
the energy radiated by a continuous beam adds 
up to zero. The phenomena observed in our 
accelerator and described below, however, are 
not affected by changing the current in the beam 
and are adequately described by the single 
electron radiation theory. It would appear that 
the fields radiated by the various electrons are 
not coherent. 

The existence of radiation is probably attribu- 
table to statistical fluctuations in the beam 
density. The average density in a random fluctua- 
tion should be proportional to the square root 
of N, the electron density, so that the energy 
radiated would be proportional to N. If this is 
the case, the single electron calculations will be 

* Cf. Tolman, Relativity, Thermodynamics and Cosmology 


Ceendon Press, Oxford, England, 1934), Section 103, 
p. 259. 


adequate to describe the behavior of the whole 
beam. 

The energy radiated should be distributed 
among the harmonics of the rotation frequency 
of the electron in its orbit, which, in our case, 
is about 57 megacycles. It would appear at first 
glance that the distribution of energy between 
these harmonics should be deducible from a 
Fourier analysis of the radiation pattern de- 
scribed by Eq. (1). Since the radiation cone has 
a width of about 27/1200, we might expect the 
greater part of the energy to be found distributed 
more or less uniformly between the first thousand 
or so harmonics. This argument, however, has 
been shown to be fallacious by Schwinger for 
reasons which he discusses at length in a forth- 
coming paper in The Physical Review.* Schwinger 
has demonstrated that the energy is distributed 
among more than 10’ harmonics, and that the 
energy distribution has its maximum in the near 
infra-red or in the visible spectrum. 

While we were under the impression that the 
energy distribution included only a thousand or 
so harmonics, we made a search over the micro- 
wave range for the expected radiation. From the 
geometry and operating parameters of the ac- 
celerator and from the calculations outlined in 
Section 4 below, it seems that a total radiation 
power of the order of one watt should be available 
for detection. The range from 50 to 1000 mega- 
cycles was searched with receivers capable of 
detecting less than 10 microwatts, but no radia- 
tion associated with the beam was detected. 
This is easily understandable in the light of 
Schwinger’s calculations which indicate that the 
power in a microwave harmonic is only about one 
part in 10° of the total energy radiated. 


4. CHARGE DENSITY IN THE BEAM AND EFFECT 
OF OTHER ELECTRONS ON SINGLE 
ELECTRON RADIATION 


We must consider three fields in evaluating the 
shape and charge density of the beam; the 
applied magnetic field whose components are 
B, and B,, the electric field due to Coulomb 
forces whose components we shall call EZ,’ and 
E,', and the magnetic field due to the current in 


6 J. S. Schwinger, “On Radiation by Electrons in a 
Betatron,” to be submitted to The Physical Review. 
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the beam whose components we shall call H,’ 
and H,’. 
In the neighborhood of the beam the applied 
field 
(8) 


to a high degree of approximation. From Max- 
well’s equations 


B,= —{Boz/Ro, 
also to a high degree of approximation. 


From Maxwell’s equations it is easy to show 
that 


B,=B(r/Ro)-, 


(9) 


H, = veoE,’, 
Hy,’ = veoE,’, 


(10) 
(11) 


for small values of (r— Ro) and z. 

When the beam is in equilibrium, we shall 
have *=2Z=0. Taking this fact into account we 
write Newton’s force equations and make the 
appropriate substitutions from Eqs. (8) to (11). 
The centrifugal force term 


mov? 
r(1—v?/c?)! 
. . . . . 
in the radial force equation we write in terms of 


the centrifugal force at the center of the beam 
(i.e., on the equilibrium orbit) thus: 


2 





—-1 
Ro(1 ry a . 
Since the fields due to the beam will disappear at 
its center, we can substitute from Eq. (3) so 
that the centrifugal force term reduces to 
—evB,(r/Ro)~. With these substitutions, in the 
neighborhood of r= Ro, z=0, the radial and axial 
force equations reduce to: 


: vBo 
E,’ (1 — eouov®) = E,’(1 —v*/c?) =——(r—Ro), (12) 
4Ro 


——- SuBoz 
E,/ (1 — eouov”) = E,'(1—v?/c*) = (13) 


0 


From Poisson’s equation and Eq. (5), the charge 
density : 


eoBov €9¢*RoBo* 
p = = 
4rRo(1—v?/c?)  4armo*v 
= 775By*c/v Coulombs per 





cubic meter. (14) 
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The number of electrons per cubic meter 


N= p/e=4.84 X 107 By%c/v. (15) 


These densities are considerably higher than 
densities achieved in conventional electronic 
devices but are still materially lower than the 
densities of 10° Coulombs per cu. m found in 
metals. 

If we know the current in the beam, our 
knowledge of the charge density enables us to 
calculate the dimensions of the beam. The 
boundary of the beam will be an equipotential. 
But from Eqs. (12) and (13), the equipotentials 
have the form 


(r — Ro)? +32" =a constant. (16) 


Evidently the beam is elliptical in cross section. 
If w represents the width of the beam in the 
plane of the orbit, then its width in the z direction 
will be w/4/3. The cross-sectional area of the 
beam will be rw?/4\/3=0.45w*. The beam cur- 
rent will be given by 


I=0.45w* pv amperes, (17) 


whence, if we substitute from (14) for p, 


w= 3.06 X 10-*J'By meters. (18) 


For a current of one ampere, at the peak field of 
0.4 weber per square meter, the beam has the 
surprisingly small radial extension of 10~* meter. 

Since Eq. (14) can be useful in estimating the 
amount of charge which can be started in the 
beam by an injector voltage V, it has been re- 
written in terms of V instead of Bo and v. The 
substitutions follow from the ordinary relativistic 
form of the force equations. Equation (14) 
becomes : 


p=4.96 X 10-8 V(V+5.11X 10°) 


x (V+10.22 105) Coul./cu.m. (19) 


For example, for 50 kv injection voltage, 
p=1.5X10~-* Coul./cu. m. 


The volume of our vacuum envelope is of the 
order of 0.1 cu. m so if the injector fills the whole 
tube with charge, the total charge present will 
be about 1.5 10-7 Coulomb. 

If this quantity of charge is injected 60 times 
per second, the average input current will be 
about 10 microamperes. If this charge were all 
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accelerated to 100 Mev, the output of the 
machine would be about 1000 watts. When the 
charge has reached approximately the velocity 
of light, the circulating current in the accelerator 
will be of the order of 10 amperes. Unfortunately, 
no technique is available for measuring this 
current, and its actual value is probably ma- 
terially lower than 10 amperes because of the 
initial transients which result in collection on the 
tube walls and because of interception of elec- 
trons by the electron source. 

We now attack the question as to the effect of 
the radiation field caused by all of the electrons 
on the radiation from a single electron. We shall 
show that this effect is so small as to be negligible. 
First we evaluate the energy density U» in the 
beam. 

In Section 2 above, it was shown that the 
radiation is mostly confined to a narrow beam 
in the direction of the electron velocity vector. 
On the average the beam of radiation will travel 
a distance of the order of (Row)! before it 
emerges from the electron beam. Hence, at any 
one time, the energy in the beam caused by 
the radiation from one electron will be about 
U,(Row)'/c. The energy density in the beam 
because of the radiation from all the electrons 
will be obtained by multiplying this quantity by 
the total number of electrons in the beam and 
dividing by the volume of the beam. This pro- 
cedure yields the relation 


Uy=NU,(Row)*/c. (20) 

We make the appropriate substitutions from 
Eqs. (15) and (18) and obtain 

Uo=1.29 X 10" J"/4By*/4U, joules/cu.m. (21) 


This energy distribution may be considered to 
be associated with an electric field. We shall 
assume the worst possible case in which this 
field acts along the radius. In this case, the field 
E,= (8% Uo/€0)! (22) 
= (2.91 X 107'J"/4B,*/4U,)! volts per meter. (23) 
To obtain the radial acceleration for the radiation 

formula, we must now change Eq. (3) to 
Mof n 


—_——— =cE,—eB,. (24) 
(1—v°/c*)} 


Squaring : 
fn? = (€/mo)*(1 —v*/c?)(E,—vB,)* 
= (e/mo)*(1—v?/c*) (E,?+0°B,’), (25) 
since £, is an incoherent a.c. field, and so its 
product with B, will average out to zero. 


We now substitute in the radiation formula 
which, from Section 2, has the form 


2X 10-7e*f,,” 


= watts. (26) 
3c(1—v?/c?)? 





When we substitute from Eq. (25), this becomes 
" 2X 10-7e2(e/my)*(E,? +0°B,”) 


A watts. (27) 
3c(1 —v?/c?) 





From this relation and Eqs. (5) and (23), if 
r= Ry and v=c, 


U,= 3.80 X 10-°By?(3.24 XK 1047 "/4*Bo*/*U,+ Bo’). 


Therefore 


U,(1—1.23 K 10-47 '/4Bo*!*) = 3.80 K10-°Bo*. (28) 


For currents as high as one ampere and mag- 
netic fields up to 0.4 weber per sq. m, the cor- 
rection term in Eq. (28) makes a difference in 
the radiation of less than one part in 30,000. 
It seems legitimate to neglect the effects of 
other electrons on the single electron radiation. 


5. DECREASE IN ORBIT SIZE CAUSED 
BY RADIATION 


The rate of change of electron energy W plus 
the rate of radiation of energy U, must equal the 
power fed into the electron by the changing 
magnetic field, so that: 

ve d¢ 


2nr dt 


Also, from Eq. (5) 
myc? 


W=———— = er B c*/0 
(1—v?/c?)! 


=erB,c approximately. (30) 


In the neighborhood of the equilibrium orbit, 
B, is given by — 


B,=Bo(r/Ro)", (31) 
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THE FUNCTION F 
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Fic. 1. The function F(wt). 


where 
Bo=A sinwt, n= —#? in our case. 


The accelerator magnet is so designed that the 
flux is given by 


r 


2nrB.dr, (32) 


g=2rRo?Bot 
: me 


n+ 
= 2nRe Bo 


1 1 
+—(/R webers. (33) 
n+2 n+2 


From (29) 
W+ U,=veRoBo(r/Ro)— 
n+1 1 
| + (r/ Ran. (34) 
n+2 n+2 
From (30) and (31) 
W =ceBoRo(r/Ro)"*!+(n+1)ceBor(r/Ro)". (35) 


In the range over which radiation will be of 
importance, we can set v=c in Eq. (34). We 
subtract (35) from (34) and rewrite to obtain: 


n+1 ecR,? d 
U,=—_ —{Bo(1—(r/Ro)"**)} watts. (36) 
mt+2rd. 
But from Eq. (7) 
U,=5.44 10-°B,*r? 
= 5.44 10-°Bo*Ro 7442, (37) 


From (36) and (37) 


d 
—{Bo(1—(r/Ro)"**)} 
at 
n+2 


= 113 Botr* +8 Ro (an +2) | (38) 
n+1 


Since By=A sin wt webers per sq. m and n= —3, 
Eq. (38) becomes: 


d 
<1 sin wt(1 = (r/Ro)*!*) } 
t 
=565A‘ sin‘ wtRo. (39) 


We now integrate Eq. (39) from ¢=0 to t=¢y and 
rearrange terms to obtain: 


212A*R,p 
(r/Ro)*/4= 1 —-—— (wtp — sin wlo COS wo 
w SIN who 


— 2 sin® who cos who) 


212A®RoF(wto) 


= ° 


Ww 





where 

F(wt) =wt/sin wt —cos wi — 3 sin? wt cos wt. 
Substitution in (40) of Ro=0.833 m, w=377, 
gives: 


(r/Ro)*!4=1—0.473A®F(wto). (41) 


If ér (=r—Ro) is small compared with Ro, 
Eq. (38) can be integrated before the substitution 
is made for n, to give 


42R,°A 3 F(wto) 
(n+1)w 





ér= meters. (42) 


With the appropriate substitutions for Ro, 1, 
and w, Eq. (42) becomes 


br = —0.315A*F(wto). (43) 


The function F is plotted in Fig. 1. This graph, 
used in conjunction with Eq. (41) or (43) gives a 
complete picture of the decrease in r caused by 
radiation losses as a function of time for our 


accelerator. For example, at 90° phase, for 


A =0.4 weber per sq. m, 6r =0.032 meter. 
In Fig. 2 we have plotted dr as a function of 
phase for A =0.4 weber per sq. m. 
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6. CALCULATION OF NECESSARY ORBIT 
SHIFT PULSE 


We have already mentioned in Section 1 the 
fact that radiation losses will make the strength 
of the pulse necessary to bring the electrons to 
the target less than that expected from a calcu- 
lation neglecting radiation effects. 

We shall make the simplifying assumption 
that the pulse is so short in duration that we 
can neglect radiation effects during the pulse 
and can assume that the value of By does not 
change during the pulse. This should be a valid 
assumption, since the first half of the pulse, 
which is the active part, lasts less than 100 
microseconds. 

The pulse will have two effects on the beam, 
first, the electron energy will be increased slightly 


because of the slight change in the flux enclosed . 


by the beam, and second, the orbit radius will 
be changed. The second effect, of course, repre- 
sents the primary function of the pulse; the first 
effect is incidental, but will be included in the 
calculations. 

The orbit shift system consists of two coils, 
one of radius 0.71 m and one of radius 0.96 m. 
The same current flows through the two coils, 
but its direction in the inner coil is opposite to 
its direction in the outer coil. Because of the 
shape of the iron core, this arrangement results 
in no field inside the inner coil, but does result in 
a field which varies with the inverse three- 
quarters power of the radius between the coils. 
We represent this field by A’ sin w’t’(r/Ro)-, 
where ?’ is measured from the beginning of the 
pulse. At the beginning of the pulse the orbit 
will have a radius r which is less than Rp because 
of radiation losses and whose value is given by 
Eq. (41). The pulse now shrinks the orbit to 
the circle on which the target lies and whose 
radius is 0.763 m. The flux due to the orbit 
shift is 


. : 
f 2rrA’ sin w't'(r/Ry)-*/*dr 
0.71 


= (8/5)2A’ sin w’t/Ro?((r/Ro)*/*—0.820). (44) 


Therefore, the rate of change of energy in the 


beam due to the pulse (cf. Section 5) 
W =(4/5)evA'w’ cos w't'(Ro*/r) 


X ((r/Ro)**—0.820). (45) 


But 


W=ecRo(Bot+A’ sin w’t’)(r/Ro)'. — (46) 


Since we assume By to remain constant during 
the pulse, differentiation of (46) yields 


W =ecR,A'w’ cos w't'(r/Ro)! 


+(1/4)ec(Bo+A’ sin w't’)(r/Ro)-**. (47) 


We eliminate W between Eqs. (45) and (47) 
letting v=c as we did in Section 5. The resulting 
equation is easily integrated between the limits 
r=r to r=0.763 m and ?t/=0 to t/=2/2. The 
result of this integration is 


Bo{ (r/Ro)*!/*+3.280} =4.176(Bo+A’). (48) 


Substitution from (41) for (r/Ro) and rearrange- 
ment of terms gives 


A’ = B,(0.0249 —0.1132A*F(wt)) 


=A sin wt(0.0249 —0.1132A*F(wt)) *erers/sa-m, 
(49) 


Figure 3 shows the orbit shift pulse peak value in 
webers per sq. m for peak magnetic fields A of 
0.2, 0.3, and 0.4 webers per sq. m. 

The observations on the orbit shift pulse have 
been made on the current flowing through the 
orbit shift coils. Direct comparison of experiment 
and theory will involve a conversion factor from 
amperes to webers per sq. m which has not yet 
been measured. The general form of the experi- 
mental curves agrees with that of the theoretical 
curves plotted from Eq. (49) ; the maxima of the 
theoretical curves occur at about the right phase 
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Fic. 2. Change in orbit radius as a function of 60-cycle 
phase for a peak magnetic fiéld of 0.4 weber per sq. m. 
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Fic. 3. Magnitude of orbit shift pulse required to move 
beam to target, as a function of 60-cycle phase. 


and the shapes of the theoretical curves agree 
at least qualitatively with the shapes of the 
experimental curves. 





7. CALCULATION OF PHASE OF COLLECTION IN 
ABSENCE OF ORBIT SHIFT PULSE 


It is evident from either Eq. (41) or Eq. (49) 
that in time the beam will arrive at the target 
even if no orbit shift pulse is applied. The phase 
of collection may be deduced by setting r = 0.763 
in Eq. (41) or by setting A’=0 in Eq. (49). In 
either case we obtain 


A*=0.220/ F(at), 
or 


A=0.604F-4. (50) 


There are no undetermined parameters in this 
relation, and we can compare it directly with 
experiment. In Fig. 4 we have plotted wt against 
A from Eq. (50). A set of experimental points are 
included on the same graph. For low values of A 
the agreement is good. At higher values the 
experimental values diverge slightly from the 
theoretical curve until a maximum disagreement 
is reached of about 12 percent in A or about 15° 
in phase. If this divergence is real, it is in such 
a direction as to indicate a slightly greater loss 
in energy than that given by radiation alone. 
It is quite possible, however, that experimental 
errors in phase measurements, magnetic field 
measurement, determination of position of the 
equilibrium orbit, and assignment of the value 
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— ? to the index of magnetic field variation with 
radius could add up to an error of this order. 
Recent experiments on the range of orbit 
shift pulse magnitudes which result in electron 
collection at the target indicate that the beam 
does not have the very small extension indicated 
by Eq. (18) above. The diameter of the beam 
would appear from these experiments to be of 
the order of two centimeters. The reason for 
this discrepancy between experiment and theory 
is not yet understood. If the beam is as large as 
this, its first interception by the target will 
occur before the time indicated by Eq. (50) by 
just about the difference between the theoretical 
and experimental points of Fig. 4. If the diver- 
gence between the experimental results and the 
predictions of Eq. (50) is real, the large extension 
of the beam evidently provides an adequate 


- explanation. 


8. TRUE ENERGY OF THE COLLECTED ELECTRONS 


Although we have shown that changes in orbit 
radius of the order of 10 percent occur in the 
range of operation of our accelerator, the associ- 
ated energy losses will generally be less than 
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Fic. 4. Phase of collection with no orbit shift pulse, as a 
function of peak magnetic field. 


2 percent. The true energy of the electrons 
striking the target will be given by (cf. Eq. (30)): 
ecRo(A sin wt+A’)(r/Ro)! joules. (51) 


If we substitute r=0.763 m, Ro=0.833 m and 
express this energy in Mev, we obtain 


True energy = 250 
X0.978(A sin wt+A’) Mev. (52) 
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The value of A’ is given by Eq. (49). Substitution 
of this value in (52) gives 


True energy = 250A sin wt 


X (1.003 —0.111A*F(wt)) Mev. (53) 


If we collect at 90° phase, F(wto) = 2/2 and 


True energy = 250A (1.003 —0.175A*) Mev. (54) 


The increase of 0.3 percent in the first term in 
the bracket represents the amount by which the 
energy would be increased because of the orbit 
shift pulse flux linked by the beam if no energy 
were lost by radiation. The second term repre- 
sents the radiation loss. As an example, we 
consider the ‘‘100-Mev” case (A=0.4 weber 
per sq. m). In this case the true energy turns 
out to be 99.2 Mev. 


9. PROCEDURE FOR ATTAINMENT OF 
HIGHER ENERGIES 


If the radiation hypothesis is definitely estab- 
lished, it will be necessary to make some changes 
in operating parameters of the accelerator before 
energies are achieved which are appreciably 
higher than those already reached. The primary 
objective will be to speed the electrons up faster 
so that they do not have so much time in which 
to lose energy by radiation. We shall also gain 
by increasing the radius of the orbit so that the 
radial accelerations are not so high. Analytically 
these factors follow from the considerations of 
Section 5. From £q. (40) we find that if the 
ratio of the target radius to Ro is to be kept at 
its present value, the ratio RoA*/w must be kept 
constant or must be less than its present value. 
Since the energy W is proportional to the product 


AR, it follows that if the energy is to be in- 
creased by increasing the magnetic flux density, 
w must be increased as the cube of the energy. 
This procedure does not look promising, both 
because of the fast increase of w and because the 
central part of the magnet is already being 
operated uncomfortably close to saturation. It 
would be far preferable to attain the higher 
energy by increasing Ry in which case w must 
increase linearly with the energy. 

Further gains will be effected by moving the 
target farther in from the equilibrium orbit and 
by letting the magnetic field fall off more slowly 
than according to the inverse three-quarters 
power of the radius (cf. Eq. (42)). Neither of 
these changes will decrease the energy loss 
caused by radiation. They will simply delay the 
time at which the beam spirals into the target 
without the assistance of the orbit shift pulse. 

One feature of the present operating technique 
which is undesirable from the point of view of 
radiation losses is the procedure of collecting at 
90° phase. During the period from 65° to 90° the 
beam energy is increased by less than 10 percent, 
yet more than half of the radiation energy loss 
takes place during this time. It would seem to 
be desirable to raise the peak magnetic flux by 
10 percent and collect at about 65° phase. 
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Means were evolved to observe the cathode dark space 
of the mercury arc, the object being to measure the thick- 
ness and evaluate therefrom the voltage gradient at the 
cathode to distinguish between field and emission theories 
of electron liberation. A magnetic field transverse to the 
arc drives it in the opposite direction to the force involved. 
This wrong way motion made it possible to race the arc 
spot over smooth mercury while ions, electrons, and vapor 
were blown rearward. Photomicrographs showed a negative 
glow, its image in the mercury, and a space between; 


(Received September 12, 1945) 


evidently twice the dark space. A one-ampere arc had a 
dark space of 0.001 cm; a hundred times too large for the 
field theory, and causing excessive space charge limitation 
of current unless compensating ionization occurs through. 
out said space. Phenomena in the negative glow must 
cause the needed ionization and also intensive electronic 
bombardment of the cathode. This is assumed to ‘cause 
cumulative excitation of the liquid resulting in emission 
of electrons and light. A continuous spectrum originates— 
within the limits of measurement—at the liquid. 











HE mercury arc cathode spot displays in- 
tensive electronic phenomena apparently 
not explained by inter-relationships of well- 
known factors. Outstanding problems are briefly 
stated. The low pressure arc goes the wrong way"” 
(opposite to the electromagnetic push) in any 
transverse magnetic field less than about 5000 
oersteds. Electrons are liberated from the liquid 
at the cathode-spot in sufficient numbers to 
supply the arc current*® of 4000 amperes per 
sq. cm. Evidence from studies of vapor pressure‘ 
yielding a surface temperature less than 200°C 
and the possibility of extinguishing the arc’ 
within 10-* sec. clearly indicate that the electrons 
are not liberated by ordinary thermionic means; 
and if the present work proves valid, the electrons 
are not extracted by electrostatic force, because 
the field is too weak. 

It is the object of the present work to view the 
free cathode spot with a microscope to find the 
dark space, if possible, and the negative glow; 
and to measure the thickness of the dark space 
and to evaluate therefrom the approximate 
voltage gradient at the cathode surface. A further 
object is to discuss the implications of the 
findings and point out in particular that the 
inferred voltage gradient at the cathode is en- 


1C. G. Smith, Phys. Rev. 62, 48 (1942); Phys. Rev. 64, 
40 (1943). 

2 N. Minorski, J. de phys. et rad. 9, 127 (1928). 

3A. Guntherschulze, Zeits. f. Physik 11, 74 (1922); 
L. Tonks, Physics 6, 294 (1935). 

4L. Tonks, Phys. Rev. 54, 634 (1938); J. V. Issendorf, 
Physik. Zeits. 29, 857 (1938). ; 

5G. Mierdel, Zeits. f. tech. Physik 17, 452 (1936). 
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tirely too small for the field theory® of electron 
liberation, and also that a powerful ionizing 
agent must be operative right at the liquid 
surface. 

The dark space has not heretofore been ob- 
served because of gross turbulence at the arc 
spot caused by blasting vapor, spraying mercury, 
electromagnetic forces, and an arc pressure’ of 
80,000 dynes per sq. cm exerted upon the liquid. 

Photo-micrographic studies of the free cathode 
spot are here made possible by exploiting two 
properties of the low pressure arc in a transverse 
magnetic field. First, the arc moves at right 
angles to the horizontal component of such field 
in exactly the opposite direction to that expected. 
This retrograde motion is pronounced at low 
pressure and persists or even increases as the 
anode is brought closer to the cathode. Speeds 








Fic. 1. Arc tube. Arc E between mercury cathode C 
and iron anode A, driven magnetically in retrograde 
direction starts at B, races past microscope lens L for 
photographing, and goes out at M. 


* 1. Langmuir, Gen. Elec. Rev. 26, 731 (1923). 
7E. Kobel, Phys. Rev. 36, 1636 (1930). 
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CATHODE DARK SPACE 


exceeding one hundred meters per second are 
attainable. The second usable property shows 
itself when the magnetic field has a vertical 
component. Besides the rapid motion due to the 
transverse component of the field, the arc moves 
slowly away from the obtuse angle between field 
and mercury surface and toward the acute angle. 
By employing a magnetic field in the form of a 
long narrow, shallow arch the lines of which 
cut the liquid surface in two points, the arc can 
be made to race along in an accurately pre- 
determined straight line under the middle of 
the magnetic archway. At a sufficient speed the 
arc cannot sensibly depress or otherwise disturb 
the serenity of the mercury surface until it has 
made its transit past a microscope. The retro- 
grade feature of the motion is particularly 
valuable because the average electromagnetic 
force upon the ions and electrons results in a 
mechanical action driving the vapor from the arc 
spot to the rear while the arc glides forward 
over smooth mercury free from disturbances 
that might otherwise have preceded it. 


APPARATUS 


Two sections of the arc tube are shown in 
Fig. 1. G is a glass tube 2.8 cm in external 
diameter, containing mercury cathode C and 
iron anode A. Pieces F and M of non-magnetic 
stainless steel are fastened together rigidly by 
two rods of the same material, imbedded in the 
mercury and shown at Z and D in the upper 
section. The anode A is supported on insulators 
K and I inserted into F and M, respectively. 
A hole through A permits passage of a tungsten 
rod T down to the mercury surface. T is fastened 
to a small iron plunger P that can be raised by a 
coil W to break contact with the mercury at B. 
V is connected to a mercury diffusion pump. 
H connects with a bulb containing mercury for 
distillation into tube G. NS represent pole pieces 
of an electromagnet each pole face of which 
extends from F to M. The anode A is wedge 
shaped, and being of iron distorts the magnetic 
field making it into an arch concave toward the 
mercury just below the anode. This arched 
field holds the arc in a line as it races along under 
the anode. Without the arched field the arc may 
wander about free from definite control by the 
position of the anode. 


a 


Fic. 2. Photomicrographs of one-ampere mercury arc 
cathode spot racing once past a microscope, with velocity 
of 70 meters per sec. 0=negative glow. /=image of 0 in 
surface of cathode. 2D=twice the cathode dark space. 
A is a direct enlargement of negative. B is enlargement 
made while sliding paper horizontally. 





A small are energized. by an auxiliary circuit 
is started at the tungsten mercury contact B 
when the iron plunger P is raised by the coil W. 
The current at B is limited to about one-half 
ampere and the voltage on open circuit is about 
40. The anode A takes over the arc current as 
the discharge at B goes out. The arc spot E 
races to the right with positive column trailing. 
The energy for the arc is supplied by a condenser 
Q of 160 microfarads. charged to 440 volts. 
Current from Q flows through R, a limiting 
resistance. The charge on the condenser is 
quickly exhausted, and the arc goes out soon 
after reaching the end of its track. This quick 
cessation of current prevents excessive distilla- 
tion of mercury. After things have calmed down 
another shot can be made. Lens L used in 
making photomicrographs is an archromat of 
3.5-cm focal length. 


PHOTOGRAPHY OF NEGATIVE GLOW 
AND DARK SPACE 


Photomicrographs were taken with a magnifi- 
cation of 7.9 and utilized a single traverse of the 
arc spot racing at 70 meters per sec. past the 
microscope. This showed two bright streaks 0 
and J, the negative glow and its image in the 
mercury, respectively. The space between the 
two is twice the dark space. Figure 2 with its 
attendant scale shows the situation for a one- 
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ampere arc. The cathode dark space is about 
0.001 cm, and the glow is roughly 0.004 cm 
thick. The light along the axis of the camera 
made an angle of 3° with the horizontal. Lens 
aperture was 4 mm. The total arc drop was 15 
volts as determined by an oscillograph. The 
magnetic field at the mercury surface is in- 
ferred from the velocity of the spot to be about 
1500 oersteds. When the axis of the camera 
made an angle of more than 5° with the hori- 
zontal, it was not possible to find a double line 
represented by 0 and J. One line was found. 
At the higher angle the reflected light must come 
from a slightly distorted mercury surface too 
near the arc spot. At the small angle of 3° the 
reflected component comes from a point some 
distance nearer the camera and sufficiently flat 
to yield a good image. 

The drawing of Fig. 3 is a pictorial representa- 
tion of the observations. The arc moves in the 
direction of the arrow S. O is the negative glow, 
and J its image in the mercury surface M. 
A possible source of light at the liquid surface is 
omitted from the drawing. D represents the 
cathode dark space and P the positive column. 

The brilliant and surprisingly thin negative 
glow must be a region of potential maximum® 
(by analogy with similar situations). It is as- 
sumed to be a region wherein a relatively high 
electronic temperature is established because it 
is the recipient of much of the enormous energy 
input of the primary electrons leaving the 
cathode surface and gaining approximately ten 
electron volts of energy before plowing into the 











Fic. 3. Pictorial sketch of racing cathode spot. 0 = nega- 
tive glow. J = image of 0 in cathode surface M. P = positive 
column. D=cathode dark space. S shows direction of 
motion. 





*K. T. Compton and C. Eckart, Phys. Rev. 25, 139 
(1925); M. J. Druyvesteyn and F. M. Penning, Rev. 
Mod. Phys. 12, 151 (1940). 
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negative glow. The essentially mechanical energy 
of these primary electrons is mostly converted 
into electronic thermal energy and radiation. 
The majority of the electrons in the negative 
glow must have a temperature, and consequently 
some of these electrons have energies in excess 
of the ten electron volts of the primaries. There- 
fore, electronic diffusion back to the cathode 
against the potential drop is to be expected. 
The magnitude of this diffusion seems to be 
about one-third of the rate at which primaries 
leave, if the electron temperature in the glow is 
45,000°K, a value found to be applicable for an 
energy balance at the cathode as previously 
considered. 

lonization in the dark space is required to 
neutralize a part or all of the negative space 
charge that would otherwise impose limitations 
upon the current density at the cathode. Without 
any ions the electronic current 7 across a space d 
impelled by a voltage drop V is 1=2.33Vd~ 
X<10-* ampere sq cm. When d=0.001 cm and 
V=10 volts, 1=74 amperes, an insignificant 
value compared to the are current of approxi- 
mately 4000 amperes per sq. cm. If now we 
assume uniform rate’ of ionic production through- 
out the dark space at the optimum rate to 
reduce space charge, we may multiply the 74 
amperes by 2.86, and we are still far below 4000. 
If more ions are generated than the optimum, 
we merely destroy the upper end of the dark 
space by creating a surplus of ions that are 
relatively slowly swept down to the cathode. 
Therefore, the ionization needed to permit the 
observed current across the dark space is appar- 
ently improbable by ordinary means though the 
ionization is evidently present. A remaining 
recourse in the dilemna is to fall back upon the 
processes that create and maintain an intense 
plasma such as the negative glow. It is the 
author’s belief that a satisfactory understanding 
of a plasma will involve more than our present 
knowledge of electronic ballistics and ion forma- 
tion. For example, there are matters such as 
the saturation current!® of a gaseous conductor, 
the inadequacy of accepted theory in dealing 
with electronic distribution in a magnetic field 
transverse® to a positive column, and the rapid 


° 1. Langmuir, Phys. Rev. 33, 954 (1929). 
© A. W. Hull, Elec. Eng. (Nov., 1934). 
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Fic. 4. Photomicrographic spectral observations 
of moving arc. Horizontal line is the continuous 
spectrum at surface of the mercury. The line spec- 
trum shows above the mercury surface and also 
below, due to reflection in the mercury. Number of 
traverses of arc past microscope for a = 2700, b = 900, 


c=300, d=100 


establishment" of an electronic temperature. 
These major effects may require additional 
knowledge for even an approximate under- 
standing. The net charge in the dark space would 
seem to be the difference between relatively 
large numbers of ions and electrons, respectively, 
which in turn are there partly as a result of 
processes of ionization and recombination in 
the negative glow, and the extension of these 
processes downward into the dark space, the 
understanding of which must involve a better 
basic understanding of a plasma. 

The liquid cathode must suffer an intense 
electronic bombardment since ions are being 
created at its very surface. Such bombardment 
was previously postulated to result in cumulative 
electronic excitation of the liquid surface, the 
liquid then becoming an emitter of electrons and 
also of light—presumably a continuous spectrum. 


SPECTRA 


A much-widened slit of a spectroscope was 
placed where the moving magnified image of the 
arc spot might cut across it. Light from the vapor 
of the arc and the reflection in the mercury 
surface give rise to the broad spectrum lines 
shown in Fig. 4, where @ was made by 2700 
traverses of the arc past the microscope, )—900, 
¢—300, and d—100, respectively. The horizontal 


"1 T, Langmuir, Phys. Rev. 26, 585 (1925); A. F. Dittmer, 
Phys. Rev. 28, 507 (1926). 
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streak is the continuous spectrum originating 
apparently at the very surface of the mercury 
and only there. The emergent light from the 
arc made an angle of about 10° with the hori- 
zontal. In order to record the line spectrum it 
was necessary to employ ten traverses of the arc 
spot past the microscope lens. The red portion 
of the continuum was observable after about a 
hundred traverses. Numerous traverses of the 
arc spot give rise to poor definition of the image, 
due probably to slight differences in level of the 
mercury and small variations in the arc track 
along the surface of the mercury. Even so, the 
line of the continuum is remarkably sharp, 
except where faint lines of HgII and Hg! super- 
posed upon it make slight bulges. The identity 
of these was determined by other observations 
not published. There is a long portion at the red 
side of \4358 that is quite narrow. Measurement 
of the vertical depth of this region together with 
the total magnification of 7.2, including that of 
the spectroscope, shows that the apparent verti- 
cal dimension of the source giving the continuum 
is about .001 cm. However, it may be that the 
horizontal width of the arc streak is really causing 
the apparent vertical dimension; and actually 
the continuum may originate in the liquid surface 
itself. If we wish to reduce the effect of the 
horizontal width, we should view the arc spot 
at a smaller angle to the horizontal, but this 
seems to reduce the total light of the continuum, 
and the necessary labor for this observation was 
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not carried out, though it should be. A more 
direct way of fixing the spatial origin of this 
light may consist in a study of the polarization 
of the continuum near grazing angle of emission, 
since a bright emitter is expected to show light 
strongly polarized” with electric vector vertical. 
The continuous spectrum does not terminate 
near \4358 as indicated in Fig. 4. Films more 
sensitive in the violet show an extension into 
the ultraviolet. 


CONCLUSION 


The moving mercury arc cathode spot driven 
rapidly in the characteristic retrograde direction 
by a transverse magnetic field has been studied 
microscopically. A thin, brilliant disk-like nega- 
tive glow, the apparent heart of the arc problem, 
has been found. It is mostly within 0.004 cm of 
the liquid surface. This finding shows that we 
have hitherto been handicapped in studies of the 
arc either with probes or by spectral means. 
The probes fail because the negative glow is too 
close to the liquid to admit of insertion of a 
probe. Spectral studies of this region have been 
frustrated by light from the relatively enormous, 
though intrinsically less brilliant, adjacent regions 
of luminous vapor. In short, both electrically and 
optically there has heretofore been little or no 
direct study of the important negative glow. 
The glow is evidently a region of potential 
maximum and this means that our previous esti- 


#2 A. G. Worthing, J. Opt. Soc. Am. 13, 636 (1926). 
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mate of about 10 volts for the cathode potential 
drop is too low. 

The observation of a cathode dark space of 
0.001 cm yields a gradient at the cathode of 
about 10,000 volts per cm. This low value makes 
the field theory of electron extraction untenable 
for the low pressure arc; and this seems to be 
true even if one should postulate pointed areas 
over the liquid, because space charge limitations 
that are already too great would be still further 
augmented, if the current distribution were con- 
strained by pointed regions. The large dark space 
would seem to make other theories untenable 
unless ions are produced throughout the dark 
space and especially at its lower boundary, to 
reduce negative space charge. The ionizing agent 
in this region is not explicitly understood. 

The micro studies have shown that the con- 
tinuous spectrum comes either from a region 
about 0.001 cm thick at the cathode surface or 
else from the cathode itself. This latter origin 
does not necessarily imply that the liquid is hot 
enough to be incandescent. The concept that 
the liquid is excited in a cumulative way! re- 
sulting from bombardment by electrons diffusing 
down from the negative glow was outlined 
previously. The present study lends support to 
that concept. 

My obligations to Harvard University in 
whose physical laboratory the work was done 
are very great. To Professor T. Lyman, the 
writer owes the opportunity to work and the 
inspiration to carry on in this field. 
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The Dielectric Constants of Eight Gases 
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The dielectric constants of eight gases have been measured by a heterodyne beat frequency 
method. Special precautions were taken to improve minor defects in the equipment which are 
believed to have contributed principally to the difference in average values obtainéd by other 
observers where these differences were incompatible with the self-consistency of the individuals’ 
data. Data on helium, neon, argon, hydrogen, oxygen, nitrogen, carbon dioxide, and air are 
given and compared with the results of numerous other observers. 





INTRODUCTION 


ALUES of the dielectric constants of gases 
have been of interest to both physicists and 
chemists over a long period of years, and many 
experimental methods for making such measure- 
ments. have been devised. These methods of 
measurement may be roughly classified as static 
methods and high frequency methods. With 
either method the dielectric constant of the gas is 
determined from a change in capacity of a con- 
denser which may be immersed in the gas at 
various pressures. In the case of high frequency 
measurements the technique usually employed is 
to determine the change in frequency of a 
resonant circuit when the gas pressure sur- 
rounding a condenser is changed, or to determine 
the change in value of an auxiliary condenser in a 
resonant circuit to maintain the frequency con- 
stant as the gas pressure is changed. The usual 
procedure for making either of these types of 
measurements is to make the condenser immersed 
in the gas a part of the resonant circuit of one 
oscillator in a beat frequency circuit. The obvious 
advantage of this system is that it permits the 
comparison of two high frequencies so that a 
small difference between two large numbers is 
readily observed. 

In nearly all cases of the measurement of the 
dielectric constants of gases the experimenters 
started by measuring air, and later they turned 
their attention to the particular gas or gases in 
which they were basically interested. Conse- 
quently, more measurements on the value of air 
are in the literature than for any other gas. 

The senior autHor started a study of the 


* Now with National Union Radio Corporation, Newark, 
New Jersey. 


measurements on dielectric constants primarily 
to determine whether there was a marked differ- 
ence between values obtained by static techniques 
and by high frequency techniques. No such 
significant difference was observed, but it was 
noticed that the self-consistency in most indi- 
viduals’ sets of data was good, whereas the 
differences in average values between observers 
was unexplainably large. It was also noticed that 
much of the work done at high frequencies did 
not employ entirely modern techniques in the 
design of the oscillating circuits and in the 
isolation of the circuits so that they could not 
couple. 

Consequently, experimental work was started 
by the senior author and Howard L. Schultz to 
see if the application of good high frequency 
technique would produce more reliable results. 
Data obtained by these observers were published 
in 1936 (11).T 

Later the experimental work was re-started 
with the aid of Carl Seddon. A number of sup- 
posed refinements were introduced in the circuits 
and in the technique employed but the results 
obtained on air differed from those obtained with 
Schultz. In other words, the data obtained with 
Schultz and that obtained with Seddon were like 
that of previous observers—very self-consistent, 
but in disagreement between the two sets of data. 

The senior author then started work with 
Donald L. Woernley and carried through ex- 
tensive modifications of details of the equipment 
and technique. In general, the data obtained by 
Woernley were in agreement with those obtained 
by Seddon. The improvement in consistency be- 


t References are given at the end of this article in the 
Bibliography. 
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tween the data of these observers indicates that 
the desired reliability is obtained by special care 
with regard to several minor details of the 
equipment and the procedure. 


EXPERIMENTAL TECHNIQUE 


The general technique employed by Schultz 
(11) was used throughout. The standard high 
frequency was obtained by the demodulation of 
a radio broadcast signal after the manner de- 
scribed by Schultz (30). The oscillator circuit in 
which the variable condenser was placed was a 
so-called electron coupled oscillator in which the 
_ resonant circuit uses the screen grid of a pentode 
as the plate and the output signal is taken from 
the plate itself. This resonant circuit was adjusted 
to zero beat with a standard frequency with the 
test condenser evacuated. The difference of fre- 
quencies resulting from the admission of gas to 
the test condenser was measured by recording the 
beat frequency with a string oscillograph. A 
standard audiofrequency was recorded on a 
second string of this oscillograph simultaneously. 

The following items are the principal ones 
about which excessive care appears to make the 
necessary difference regarding the reliability of 
the data. 

It is desirable to make the case of the condenser 
of very heavy material to avoid minute changes 
in dimensions caused by pressure changes inside 
of this case as the gas is admitted. Changes 
caused by this effect were further reduced by 
using a perforated shield inside of the outer case. 

It is also necessary to take precautions to keep 
the temperature of the test condenser and its case 
constant, not only to simplify the problem of 
pressure and temperature measurements, but 
also to prevent changes in mechanical dimensions 


L. WOERNLEY 


of the condenser and its case. This problem was 
solved by building the test condenser external to 
the shield enclosing the oscillator since tempera- 
ture difficulties would otherwise arise caused by 
the stabilization of the oscillator temperature 
above that of the room in which the gas supply is 
kept. 

Another problem arises from the fact that it is 
difficult to build a condenser in which no part of 
the capacity is so located that it cannot be 
affected by the change in pressure of the gas, 
This problem is inherent since it is necessary to 
use solid insulating material to support the live 
plates and to surround the lead from the live 
plates as it comes through the shield and outer 
casing. It is possible, however, to take data under 
sufficiently varied conditions so that the un- 
changing part of the capacity can be eliminated 
from the calculations. It is, nevertheless, desirable 
to keep this unchangeable capacity as small as 
possible in order to improve the over-all reliability 
of the data. 

Figure 1 shows the exterior appearance of the 
apparatus directly associated with the variable 
frequency oscillator in the latest version of this 
equipment. At the left we see a steel cvlinder 
with a wall 2?” thick and with 1” ends. It was 
tooled from a solid piece of steel. Inside this 
cylinder is a second cylinder, this one being of 
brass with perforated walls and ends. Gas pres- 
sure is changed inside of the steel cylinder and so 
changes on both sides of the perforated brass 
cylinder. The latter, then, should not experience 
serious structural changes as the gas pressure 
changes. It serves ‘as an unchanging electrostatic 
shield for the condenser located inside of it. The 
live plates of the multiple plate condenser are 
mounted on three isolantite pillars. The plates 
are arranged so that various numbers of them can 
be used without changing the amount of the 
capacity in the solid dielectric. Data taken with 
two sizes of the condenser permit us to determine 
the value of the non-effective part of the test 
condenser capacity—in our case it is approxi- 
mately 8/10 of 1 micro-microfarad in comparison 
to a total capacity of approximately 647 micro- 
microfarads. All the parts of this condenser are 
gold-plated to establish fixed surfaces of high 
conductivity. 

The large section in the center of the apparatus 




















is a shield box made of §” aluminum. It contains 
vernier condensers and an inductance coil wound 
on an isolantite form. The small thermionic tube 
mounted on this section is a single-ended radio 
frequency pentode with the characteristics of 
sharp cut-off plate current. This type of tube 
permits all of the wiring to be in the shield 
beneath the tube, and at the same time all of the 
heat developed by the tube is external to the 
shield. . 

The final section of the apparatus is a buffer 
amplifier which effectively prevents coupling be- 
tween the variable frequency oscillator and the 
source of fixed frequency. The buffer is tuned to 
the second harmonic of the oscillator to further 
reduce coupling and to double the beat note. The 
temperature of the entire room was held constant 
to +0.1°C by automatic thermostatic control 
during the taking of a set of data. 





DISCUSSION OF DATA 


Table I is a summary of most of the values for 
air that are available. The value 1.000600 at the 
beginning was obtained by a high frequency 
method as was the value 1.000540 at the end. 
The value 1.000590 was found by Boltzmann (1) in 
1874 and the second value of 1.000590 is at- 
tributed to Schultz as measured in 1934 and 
published in 1936. The closest agreement to the 
results of Seddon and Woernley is that obtained 
by Zahn (8), although a static value obtained by 
Tangl (12) is also very close. 

Table I] shows a summary of the data on 
nitrogen. It may be noted that except for the one 
very low value there is much better agreement 
among observers for nitrogen than for air. 


TABLE I. 14 values for dielectric constant of air at N.T.P. 
arranged in order of magnitude. 














Value Year Ref. Value Year Ref. 
1.000600 1922 6 1.000586 1885 2 
594 1927 3 579 1929 23 
592 1931 4 576 1931 18 
591 1934 5 573 1908 12 
590 1874 1 572 1924 8 
590 1936 11 547 1928 10 
589 1928 9 540 1924 7 
1.000580 Av. 
1.000569 Hector and Seddon 
1.000567 Hector and Woernley 








DIELECTRIC CONSTANTS 


103 








The most recently published data for nitrogen 
are those obtained by Fox and Ryan (20) in 1939. 
Their value agrees closely with ours. Again it may 
be noted that the value of Zahn (8) lies close to 
these values. 

Table III shows values for oxygen. Our value 
for oxygen is higher than three of the four values 
shown here, whereas our values for air and nitro- 
gen were lower than those of many observers. 

It may be of some interest to compute the 


TABLE II. 10 values for the dielectric constant of mitrogen 
arranged in order of magnitude. 











Value Year Ref. Value Year Ref. 


1.000590 1931 14 1.000581 1908 12 








589 1931 16 581 1924 8 
587 1913 13 579 1939 20 
587 1931 18 573 1932 17 
582 1934 27 555 1924 7 

1.000580 Av. 






Hector and Seddon 
Hector and Woernley 


1.000581 
1.000580 








value of air from the values of nitrogen and 
oxygen already shown and that of argon given in 
Table IX. Table IV shows the results of such a 
computation. It can be seen that the value for air 
obtained by this computation lies within 5 points 
in the 4th place of the value obtained by direct 
measurement. This agreement is within the 
uncertainty which we place on any of these data. 

Table V shows a summary of published values 
for the gas hydrogen together with the values 
obtained by us. 

Table VI is a set of data for carbon dioxide. 
Measurements by Fox and Ryan in 1939 indicate 
that the dielectric constant of carbon dioxide 
shows an anomolous rise in the temperature 
region 25 to 40 degrees. We have, therefore, 
shown here only such data of other observers as 
were taken near the temperature employed by us. 

Table VII shows values for the gas helium. The 
first four are theoretical calculations. Three of 
the four are in fairly close agreement among 
themselves and also with the value obtained by 
extrapolating optical data (22). Our value is in 
fair agreement with the theoretical values and 
with the one obtained from the optical data. For 
a discussion of the theoretical values in Table 
VII see reference (31). 
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TABLE III. Dielectric constant of oxygen at N.T.P. 
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TABLE VII. Dielectric constant of helium at N wy 2 





ae 








—= 





Value Year 


1.000507 
507 1901 
518 
547 


1.000520 


1.000525 
1.000523 





Hector and Seddon 
Hector and Woernley 








TABLE IV. Air calculated from the measured values of 
its components (CO, and H,0 free). 








Contribution 


5796 4525 
5233 1098 
5451 $1 


1.0005674 
1.0005669 
+10 


Percent Value 


78.06 
21.00 
0.94 


Substance 


Nitrogen 


Oxygen 
Argon 





Measured value 








TABLE V. Dielectric constant of hydrogen (He). 


Theoretical values 


1.0000653 by Atanasoff 
691 Hassé 
1.0000790 Hassé 
715 Slater and Kirkwood 


712 Av. 


Experimental values—Ref. 
1.0000693 From optical data 
1.0000740 1908, 21 

728 1934 26 
720 Av. 


1.0000684 Hector and Woernley 








TABLE VIII. Dielectric constant of neon at N.T.P. 








Year Ref. 


1929 24 (based on air at 
1.000589) 


Value 


1.000148 
( 139) 
134 1934 26 





1.000127 Hector and Woernley 








Value Year 


1.000263 
264 
265 
273 
275 


1.000268 


1.000271 
1.000272 





Hector and Seddon 
Hector and Woernley 











Value 


Year Ref. 
1.000575 1927 25 (based on Oz at 
( 550) 


1.000547) 
574 1929 24 (based on air at 
( 553) 


1.000589) 
550 1934 26 





1.000545 Hector and Woernley 





TABLE VI. Dielectric constant of carbon dioxide at N.T.P. 
(data taken at temperatures indicated). 








Value Year 


1.000982 1926 
983 1931 
987 1934 
994 1939 
996 1939 


1.000988 
1.000988 





23.9 Hector and Woernley 








Table VIII is a set of data for neon. Our value 
is somewhat lower than that of other observers. 

Table IX is a set of values for argon. Again our 
value is relatively low. 

The values in parentheses in Tables VIII and 
1X are from data given in references (24) and (25) 
after correction to our values for air and oxygen. 





The more common impurities to be found in 
each of the rare gases are the other rare gases and 
nitrogen. Such impurities would, of course, tend 
to increase the measured values, especially for 
helium and neon. The three rare gases used by us 
were specially prepared and purified for these 
measurements by the Linde Air Company re- 
search laboratories, through the courtesy of Dr. 
Leo I. Dana. 

Table X is a résumé of our data on all of the 
eight gases which we have so far measured. The 
second column shows the number of observations 
and the third column the temperature range at 
which the data were taken. The fourth column 
shows the smallest and the largest values ob- 
tained and so gives the maximum variations. 

The reliability factors which we have attached 
to the final values are determined arbitrarily. 
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They are approximately twice the computed 
uncertainties as determined from the variations 
in the data and from the reliability with which 
individual constants in the apparatus could be 
determined. 
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Tempera- 
Number ture 
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23.85 
23.9 
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Argon 
Hydrogen 
Oxygen 
Nitrogen 
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The paper deals with the electromagnetic effect of 
motions in the earth’s core, considered as a fluid metallic 
sphere. On the basis of simple estimates the electric con- 
ductivity of the core is assumed of the same order of mag- 
nitude as that of common metals. The mathematical 
treatment follows Hansen and Stratton: three independent 
vector solutions of the vectorial wave equation are intro- 
duced; two of these have vanishing divergence, and they 
are designated as toroidal and poloidal vector fields. The 
vector potential and electric field are toroidal, whereas the 
magnetic field is poloidal. These vectors, expressed in terms 
of spherical harmonics and Bessel functions, possess some 
notable properties of orthogonality which are briefly dis- 
cussed. The theory of the free, exponentially decaying 
current modes is then given, leading to decay periods of 
the order of some tens of thousands of years. Next, the 
field equations in the presence of mechanical motions of 
the conducting fluid are set up. The field is developed in a 
series of the fundamental, orthogonal vectors, and the 
field equations are transformed into a system of ordinary 
differential equations for the coefficients of this develop- 
ment. The behavior of the solutions depends on the sym- 
metry of the “coupling matrix” that arises from the term 
of the field equations expressing the induction effects. In 


order to evaluate this matrix the velocity field is developed 
into a series of the fundamental vectors similar to the 
series for the electromagnetic field. It is then shown that 
when the velocity is a toroidal vector field the coupling 
matrix is antisymmetrical. When the velocity field is 
poloidal, the coupling matrix is neither purely symmetrical 
nor purely antisymmetrical. For stationary fluid motion 
the linear differential equations can be integrated in closed 
form by a transformation to new normal modes, whenever 
the matrix of the system is either symmetrical or antisym- 
metrical. In the latter case the eigenvalues are purely 
imaginary and the coefficients of the new normal modes 
are harmonic functions of time, representing oscillatory 
changes in amplitude of the field components. For a sym- 
metrical matrix the eigenvalues are real and the time 
factors of the new normal modes are real exponentials 
representing amplification or de-amplification as the case 
may be, depending on the sign of the velocity. For a 
matrix without specific symmetry, normal modes do not, 
as a rule, exist but similar, somewhat less stringent results 
can be derived in special cases. In the case of toroidal flow, 
in particular, the oscillatory changes of the field com- 
ponents are superposed upon the slow exponential decay 
characteristic of the free modes. 





HE existence of a fluid metallic core in the 
interior of the earth seems sufficiently well 
established! as a result of seismic observations to 
be used as a starting point for the explanation of 
some of the more outstanding phenomena of 
terrestrial magnetism. From this viewpoint the 
metallic core is the place where the electric 
currents flow that are the sources of the field. 
The secular variation of the field is interpreted 
as a modification of the current system caused 
by inductive interaction between mechanical 
motions of the fluid and the magnetic field. This 
interaction is expressed mathematically by the 
differential Eq. (32), below, which can be solved 
by the methods of electromagnetic wave theory. 


* The completion of the work presented here has been 
delayed owing to several years of war research work. 

** Now with RCA Laboratory, Princeton, New Jersey. 

1B. ne, ed., The Internal Constitution of the 
Earth (McGraw-Hill Book Company, Inc., New York, 
1939). 


One need not make any particular assumptions 
about the magnitude and form of the fluid 
motions but can try to determine these as far as 
possible from the analysis of the secular variation 
of the field. This attitude is no doubt somewhat 
unsatisfactory as one might expect that informa- 
tion about the character of the motions could be 
derived from general dynamical principles. Un- 
fortunately, physical hydrodynamics is still a 
rudimentary subject in many respects. A few 
remarks on this topic will be found in the ap- 
pendix immediately following Part II of this 
paper. We shall also abstain from speculations 
about the cause of the motions. The magnitude 
of the velocities derived from the secular varia- 
tion is indicated below; the corresponding kinetic 
energy is an exceedingly small quantity when 
compared to any thermodynamical energy of 
physical interest. The rate at which this energy 
is converted into heat by turbulence is probably 
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extremely slow and the minute power needed to 
maintain such motions can be supplied by a 
number of possible processes (radioactive or 
thermodynamical). 

Among the various forces acting upon the 
fluid in the earth’s core there is one that cannot 
be overlooked; that is the mechanical reaction 
(magneto-mechanical force) of the electric cur- 
rents upon the fluid. If, as we shall see later, 
amplification of the field by the fluid motion can 
occur then, with sufficiently strong amplification, 
the field must increase exponentially in time. The 
“braking” action necessary to slow down and 
eventually stop such an increase does not arise 
from the inductive process itself, but will be 
brought about by the magneto-mechanical forces. 

Although the electric conductivity of the core 
is probably a function of depth, it does not seem 
plausible that this variation introduces phe- 
nomena that modify the theory in a fundamental 
way, and for the sake of simplicity the con- 
ductivity will be assumed constant throughout 
the core. We shall use a value, c=10° mhos/ 
meter which is one-tenth of the conductivity of 
iron under ordinary laboratory conditions. Geo- 
chemists have brought forth a considerable 
amount of plausible argument to the effect that 
most of the matter in the core is metallic iron. 
The high temperature decreases, but the com- 
pression increases the conductivity, as will be 
shown in the appendix following Part II where 
a justification of the numerical value adopted 
will be given. The results of the theory presented 
here are not substantially changed, however, 
should the conductivity be several times smaller 
or larger than the value assumed. 

The magnetic susceptibility of the core will be 
taken as u=yo, the susceptibility of space. 

We might remark at this point that in the 
matter of units and dimensions we have con- 
sistently adopted a rationalized m.k.s. or Giorgi 
system, following Stratton.” 


THE FUNDAMENTAL VECTORS 


Maxwell’s equations, neglecting the displace- 
ment current, are 


VxE+0B/at=0, wxB—yoE=0. (1) 
(McGraw-Hill 


2J. A. Stratton, Electromagnetic Theor 
hap. 1. 


Book Company, Inc., New York, 1941), 


MAGNETISM 


On introducing the vector potential, A, by 
B=y xA, “E=-—0A/dt, veA=0, (2) 
and using the abbreviation? 


(3) 


one obtains in the usual way the differential 
equation 
(4) 


Outside the conducting sphere this reduces to 
v*A=0. (5) 


Numerically, we have in the m.k.s. system with 
the value of o given before 


v’A=Vv-A-—V xv xA, 


v*A—yo 0A/dt=0. 


uo = 1.26 sec. /meter?. 


It is a well-known fact that freely decaying 
currents can exist in a conducting sphere,* and 
these will be derived below. On putting 


A=Aye™*, (6) 


Ayo = k?, 
Eq. (4) reduces to 


V7Aot+k*Ao=0, (7) 


which is identical in form with the familiar wave 
equation. 

In order to obtain solutions of (7) we apply a 
method due, in this form, to Hansen.** In the 
field of terrestrial magnetism similar methods 
have been used to treat induction in the earth's 
crust.* Let yw be a scalar subject to the wave 


equation 
Vytky=0. (8) 


There are three solutions of the vectorial wave 
Eq. (7) corresponding to any given solution of 
(8). They can be constructed as follows. The 
first solution is 

(9) 


where R is a constant of the dimension of a 
length which will be taken later as the radius of 


U=Rvy, 


3 W.R. Smythe, Static and Dynamic ee ae 


Hill Book Company, Inc., New York, 1939), Chap. 11. 

*W. W. Hansen, Phys. Rev. 47, 139 (1935). 

5 See reference 2, Chap. 7. 

6S. Chapman and T. r Whitehead, Trans. Camb, Phil. 
Soc. 22, 463 (1923); S. Chapman and A. T. Price, Phil. 
Trans. 229, 427 (1930); A. T. Price, Proc. London Math. 
Soc. 31, 217 (1930); 33, 233 (1932); B. N. Lahiri and A. T. 
Price, Phil. Trans. 237, 509 (1938). 
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the earth’s core. Since (9) is an irrotational vector 
it is not represented in the magnetic field. The 
second solution is 


T=v xry=vy xr, (10) 


where r designates the radius vector from the 
origin. The third solution, finally, is 


S=Ry xv xry. (11) 

The last two solutions are connected by the 
identities® 

S=Ryv xT, T=(Rk’)"'yvxS. (12) 


It should be noted that vectors of the types 
(9)-(11) may be derived from an arbitrary 
scalar, not necessarily fulfilling the wave Eq. (8). 
The second relation (12) no longer holds in the 
general case. We shall, however, not use com- 
pletely general expressions for y but shall always 
assume the dependence upon the polar angles in 
the form of a spherical harmonic. Fundamental 
vectors of this type will be used to represent the 
fluid motion. 

We shall now introduce names for these three 
types of vectors. They will be designated as 
scaloidal (U), toroidal (T), and poloidal (S) vector 
fields. The electric field and vector potential 
pertaining to a poloidal magnetic field are 
toroidal, and vice versa. 

On writing down components in spherical coor- 
dinates one finds 


Ti =0, Tw =(sin @)—! dp/dg, 
Tw =- aan 06, 


(13) 


1 oy? 


=n(n+1)Ry/r 14) 
(for Y~nth order harmoni¢), 
Sim = (R/r)d*(rp) /0rd8, 
Sc») = (R/r sin 0)0?(np) /drd¢. 


The function y from which these vectors are 
derived fulfills, in the case of the electromagnetic 
field vectors, the wave equation (8). Accordingly 
we introduce® a system of fundamental solutions 
of (8) whose component functions are Bessel 
functions of r and spherical harmonic functions 
of 6 and ¢. 
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Von™ =const -7~4Tns s(Ren?) Yn™(0, ¢). (15) 


As will be shown in a later section, the k,, are 
determined by the condition 


In—y(RenR) =0, 


where R is the radius of the sphere considered 
(radius of the earth’s core). 


(16) 


ORTHOGONALITY AND NORMALIZATION. 
COMPLEX VECTORS 


This section contains some mathematical en- 
largements on the subject of the orthogonality 
of the fundamental vectors and some related 
problems which will be found useful later on. 
The normalization of the toroidal vectors is given 
in formulas (21)-(24). Thereafter the more 
physical aspects of the problem will be resumed 
in the next section. 


Orthogonality 


The scalar functions (15) form an orthogonal 
set for the interior of the sphere of radius R. The 
vectors U, T, S derived from y in the way de- 
scribed also have notable properties of ortho- 
gonality.® We shall in the first line be interested 
in orthogonality on integration over the polar 
angles @ and ¢. 

It is convenient to introduce temporarily an 
operator V, that is the projection of the operator 
Vv upon the surface of the unit sphere. Thus, if 
Y is a scalar function of @ and ¢ 


(V2¥Y)@=9Y/00, (V2¥))=(sin 0) Y/d¢, 


and if X is a two-dimensional vector in the surface 
of the sphere 


V2°X = (sin 6)—'0(sin 0X) /d0 
+ (sin 0)—'dX.~)/d¢g 


and identically for any vector X, on integrating 


over the sphere 


(17) 


J vexdo =0 


which is the two-dimensional form of Gauss’ 
theorem applied to a closed surface. From (17) 
we can derive two-dimensional cases of Green’s 
formulas, for instance (where @ and @ stand for 
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two sets of indices) 


0= f v2*(Y(a)¥2¥(8))do 


= f v:¥(a)*w2¥(8)do 


— 14(t09+1) f V(a)¥(8)de. (18) 


Similarly one obtains the identity which will be 
used later 


f Y(a)¥2¥(8)*v2¥(y)do 
+ f V(8)¥2¥(a)*W2¥(y)do 
= thy (tty-+1) f V(a) ¥(8)¥(y)de. (19) 


Later on we shall make use of complex spherical 
harmonics as defined in (22). These are ortho- 
gonal in the Hermitian sense; i.e., when the 
asterisk designates the conjugate complex quan- 
tity 


f Y(a)¥*(8)do=0, if ax. 
On applying the identity (18) one can show that 


any two different toroidal vectors are orthogonal. 
One finds readily from (13) and (15) 


f T(a)*T*(8)de = f(r) f v:¥(a)*¥2¥*(8)do 


= f(r)np(g+1) f Y(a)¥*(8)do, (20) 


where f(r) stands as an abbreviation for the 
product of the two radial functions. 

The orthogonality of the radial functions 
follows from the formula proved in the theory 
of Bessel functions, 


(ka* — kg*) f J nati (Rat) Jng+s(ker)rdr 


= ker Jna+i(Rar)Jng—s(Rer) 
— karIng—}(Rat)Jng+i (Rar) 


in conjunction with (16). 












The orthogonality of any two different poloidal 
vectors, S, with respect to integration over the 
angles @ and ¢ can be proved by a procedure 
analogous to the one just applied. Orthogonality 
of the radial functions of S will not be of im- 
portance later on so that we need not enter into 
an analysis of the behavior of these functions 
which is somewhat more complicated. The 
orthogonality of two different scaloidal vectors, 
U, with respect to integration over the angles @ 
and ¢ may also readily be demonstrated. 

Any vector §S is orthogonal to any vector T; 
from (13) and (14) the integral over the sphere 
has the form 


f S(a)*T*(8)de 





dY(a) a¥*(6) _9¥(a) a¥*(6)] de 
“Se ls de 00 |sin 6 


and the integral on the right-hand side may be 
shown to vanish by means of integrations by 
parts. There is, however, no full orthogonality® 
of the vectors S and U, but we shall not be 
concerned with this question later on. 


Normalization 


The toroidal vectors, T, will be normalized 
with respect to integration over the interior of 
the conducting sphere: 


frrev=1. (21) 


The vectors S are derived from T by (12), and 
are therefore in general not normalized simul- 
taneously with the T. We shall use complex 
spherical harmonics and set: 





2n+1 (n—m)!}} 
4an(n+1) (n+m) ] 
xXP,™(cos @)e™*, (22) 


¥."(8, ¢) -| 


Here, m goes from —n to +n. The normalization 
factor appearing here is slightly different from 
the one used to normalize the scalar functions, in 
order to take account of the factor m(m+1) in 
(20). The radial functions, normalized to unity, 
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are 
Zn(ker) =2*R-“ [Ing a(keR) Sng a(har), 
and thus the normalized function y becomes 


Ven™ =Z,(kr) Y,"(4, ¢) . 


(23) 


(24) 


It may be noted here that for nm =0 the vectors 
T and S vanish, so that the dipole, »=1 is the 
lowest spherical harmonic that appears in these 
vectors. 


Complex Vectors; Symmetry 


The use of complex quantities in electromag- 
netic theory is of course commonplace, but in 
view of the special use which will be made of 
them below, some remarks about Hermitian 
symmetry might seem justified. The relation 
between a complex vector and its real counter- 
parts may be written symbolically 


T,*" =a T,”°+iT,™, 


where the upper indices c and s refer to the 
cosine and sine functions, respectively. If an 
arbitrary toroidal vector is developed into a 
series, >> c,"T,”, the relation between the real 
and complex coefficients is 


Ca>™ = Cr F164". 


We shall, in particular, need developments of 
vectors of the form QT,,”’ where Q isa (vectorial) 
-operator. The coefficients of development are of 
the form 


c(n, m3 n'm!) =F f Y,*"QY,."de, 


where F designates the integral over the radial 
functions. These coefficients are linear combina- 
tions of four coefficients of the corresponding real 
vectors, which are readily derived from the 
preceding formulas. These coefficients c form a 
matrix, and the symmetry properties of such a 
matrix within the Hermitian symmetry will turn 
out to be of considerable importance. Within the 
Hermitian symmetry a matrix M will be called 
symmetrical when 


M (a, 8) = M*(8, a) 
and antisymmetrical when 


M(a, 8) = — M*(B, @), 
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where a stands for n,m and 8 for n’,m’. To any 
antisymmetrical matrix there is related a sym. 
metrical matrix by 


Muym =1M aot. (25) 


The eigenvalues of a symmetrical matrix are al] 
real numbers; consequently the eigenvalues of 
an antisymmetrical matrix are all purely imag- 
inary numbers. 


FREE MODES 


We shall now determine the freely decaying 
current modes in a conducting sphere.* The 
vector potential is a toroidal vector, and ac- 
cordingly we set in the interior of the sphere, for 
an individual mode 


(Asn™) =Con™(Ten™) exp (—Arnt), (26) 


where, by (6) 


Aon = Ren?/yco. (27) 


The generating scalar is given by (24); the c.,” 
are constants. Similarly we set in the space ex- 
ternal to the sphere 


(Aan™) © = Con™(Tan™) © exp (—Aent), (28) 
where T® is derived from a scalar 
ne =" Y,"(9, ¢). 


The magnetic and electric field are given by (2); 
for the internal field they are (on dropping for 
simplicity the index 7) 


Bun™ = (Con™/R)Sen™ Exp (— Asn) 


an = Con™Aca hea” exp ( = Ant) 


and similar expressions for the external fields. 
The boundary conditions require continuity of 


Biwy, Bo/u, Bi/u, Ew, Ew 


at the boundary of the sphere, r = R. E,,) vanishes 
according to (13) so that no boundary condition 
for this component is required. We assume that 
there is no discontinuity of u at the boundary. 
On writing down the relations expressing equality 
of the internal and external components of the 
field at the boundary, one is led to the condition 


0/Orl Tina y(Ren?) |rar —MP Ing i(RenR) =0. (29) 
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This is the characteristic equation whose roots 
are the k,n; after some transformation it can be 
put into the simple form (16). Thus we have, 


using (27), 
Jn—\(Xen) =O, Ren =Xen/R, 


Bow = Xen?/ Ryo. 


(30) 


For n=1, the roots of (30) are x,=sr. The 
longest period of decay for the earth’s core 
(R=3.5 X10* meters) is 


(Ay,)~'~ 50,000 years 


with the value of « adopted here. The length of 
this period is proportional to «. An appropriate 
physically significant average over the lowest 
decay periods will be several times smaller than 
this value and will therefore be of the order of 
some tens of thousands of years. As will be shown 
in Part II, the intervals of time characteristic of 
the secular variation of the earth’s field (more 
precisely, the periods of its predominant Fourier 
components) are of the order of magnitude of a 
few hundred years. Hence the periods of free 
decay are (for the lower spherical harmonics) 
usually large compared to the periods of the 
secular variation. 

The solutions given here, representing a slowly 
decaying magnetic field, are distinct from the 
well-known ‘‘magnetic’”’ modes of electromag- 
netic oscillations of a sphere. The difference 
resides not so much in the fact that the displace- 
ment current is neglected here, but in the 
boundary conditions at infinity which require 
that the field vanishes in our case, whereas in 
the case of oscillations an outgoing wave is 
required. The lowest oscillatory modes corre- 
spond to wave-lengths that are comparable with 
the diameter of the core; the frequencies are of 
the order of some ten cycles per second. 

A second system of solutions may formally be 
constructed by assuming B toroidal, A and E 
poloidal (corresponding to the ‘‘electric’’ modes 
of oscillations). If one tries to satisfy the bound- 
ary conditions, however, it appears that they 
cannot be simultaneously fulfilled. One is prob- 
ably justified in concluding that exponentially 
damped, free toroidal modes cannot exist inside 
a conducting sphere. It is perhaps premature to 
infer that no magnetic field of the toroidal type 
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can exist or can be generated in the sphere by 
agents operating inside it. In this paper the pos- 
sibility of a toroidal magnetic field is not con- 
sidered, but the subject seems to be open to 
further investigation. 


THE DIFFERENTIAL EQUATION OF INDUCTION 


We shall now proceed to the analysis of the 
effects of fluid motion upon the magnetic field. 
If the motion is described by a velocity vector v, 
the second of Maxwell’s Eqs. (1) is replaced by’ 


v xB—yucE =yov xB. (31) 


On retaining (2) for the vector potential we have 


v’A—yuo 0A/dt=—povx(yxA). (32) 


Before entering into the mathematical analysis 
we may obtain an estimate of the order of mag- 
nitude of the velocities required for an effective 
inductive action. A figure representing a lower 
limit may be gained by setting the term on the 
right-hand side of (32) equal to the value that 
each of the terms on the left-hand side has for 
the lowest free mode. This gives for the order of 
magnitude 

v~LlAu, 


where L has the dimension of a length. On 
putting L = 10° meters, 


v~6.5X10-7 m/sec. ~0.2 mm/hour. 


Since, however, the periods of the secular vari- 
ation of the magnetic field are, in general, several 
hundred times shorter than the lowest period of 
free decay, the actual velocities prevailing in the 
core will be several hundred times larger than 
this value, or of the order of several centimeters 
per hour. 

In order to integrate Eq. (32) we develop the 
vector potential 


A=), ¢,(é)T,, 


where the symbol y stands as an abbreviation 
for the triple s, m, m. On introducing this into 
(32), multiplying with T,* and integrating over 
the interior of the conducting sphere there 
follows by (21) and (27), since T fulfills the wave 


(33) 


7 The field equations for this case may be derived from 
the field equations for moving systems (Minkowski's 
equations) to be found in textbooks on relativity. 
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Eq. (7), 
A,c,+dc,/dt= >, afv x (v7 x Ts)eT,*dV. 








This will be re-arranged to read 
dc,/dt+A,c,= > 2 oR f vB, *xT,*dV. (34) 


The physical behavior of the solutions is deter- 
mined by the character of the coupling matrix on 
the right-hand side. In general the velocity v will 
be a function of time and the integration of the 
system of differential Eqs. (34) is then quite 
difficult. If the fluid motion is stationary so that 
v is independent of time, the matrix elements 
become constants. The integration can then be 
performed by introducing normal coordinates, 
i.e., linear combinations of the c’s which trans- 
form the coupling matrix into a diagonal matrix 
while retaining the diagonal form of the left-hand 
side. The transformation to normal coordinates 
will be taken up in the last section. 


THE COUPLING MATRIX 


In studying the coupling matrix we will have 
reference to symmetry of the Hermitian type as 
explained previously. We now develop the 
velocity field in a series 


V= Dea (VaS(a) +WaT(a)+uUaU(a)). (35) 


The components U are carried along here for the 
sake of completeness; they do not vanish since 
the fluid in the core of the earth is known to be 
compressible. According to seismic observations 
the density p of the core increases appreciably 
from the boundary downwards. The equation 
of continuity, V«(pv) =0 gives readily 


Vv = — V7 9p/ pdr. 


Since ¥VeT = ¥V*S =0, this shows that the U-com- 
ponents actually appear in the series (35). We 





[T(a)*S(8) x T*(y) ]=[S(6)-T*(y) x T() ] 
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shall, however, assume that the effects. of com- 
pressibility are small and shall confine our atten- 
tion primarily to the T and S components of the 
velocity. 

Now introduce the abbreviation 


f wSXTdV=[veSxT*]. (36) 


There are three types of elements of the coupling 
matrix, namely, 


[T(a)-S(6) x T*(y) ], 
[S(a)-S(8) x T*(y) J, (37) 
[U(a)-S(6) x T*(y) ], 


which will be considered in turn. Here and later, 
(a) refers to the fluid motion, (8) to the “pri- 
mary” magnetic field, and (y) to the “‘secondary” 
or induced magnetic field. 

It is convenient to use uniform mathematical 
expressions for the three vectors involved in each 
of the matrix elements (37); thus all vectors will 
be derived from generating functions of the form 
(24). This imposes an unnecessary restriction 
upon the velocity vector since the scalar (24) is 
subject to the wave Eq. (8), whereas no such 
condition is required for the velocity field. It is 
convenient to use the same spherical harmonics 
for the angle functions of the velocity vector as 
for the electromagnetic field, but the radial 
functions of the velocity vector are not subject 
to the Bessel differential equation. Notwith- 
standing their outward similarity to the field 
components these radial functions, Z(a), are to 
be considered as arbitrary functions, vanishing 
at the boundary. 


Toroidal Flow 


We consider the first of the three types of 
matrix elements in which the velocity is a 
toroidal vector. On writing this element out in 
spherical coordinates we find from (13) and (14) 





Y(a) dY*(a) _9¥(a) d¥*(a)) do 
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where 
F=R f Z(a)Z(8)Z(y)r*dr. (39) 


Now we must have 


mM.+mg =mMy,, (40) 


otherwise the integrals vanish (‘‘selection rule’ 
for the index m). The first integral in (38) may 
be transformed by means of an integration by 
parts; if x =cos @ we have 


dY(a) 


fy ve 


It is seen from (38) that the integral considered 
only appears when m,+0; in this case Y*(y) =0 
at the boundaries, so that the contribution of 
the boundaries in the last formula vanishes. By 
use of this formula and (40) the matrix element 
can be written 


[T(a)*S(8) x T*(y) ]=na(mp+1)F 


(y)dx = | ¥(a) ¥(8) ¥*(y) | 23 


P va(y)de— f ¥(a) ¥(8) as 


dni f Y(a)— OPE) vec de 


ro f ¥(a) rte) (41) 


The factor mg(ms+1) which disturbs the sym- 
metry of this formula can be made symmetrical 
by the substitution, in the equations of motion 
(34), of new coefficients — 


d,=(n,(ny+1) }ic,, 
whereby this factor changes to 
[(mp+1)n,(n,+1) }}. 


The matrix composed of the elements (41) then 
becomes antisymmetrical in the Hermitian sense: 
when the indices 8 and y of the primary and 
secondary field are interchanged the matrix 
elements (which are purely imaginary) remain 
unchanged, i.e., they change into the negative 
of their conjugate complex. We thus have the 
result: For toroidal fluid motion the coupling 
matrix ts antisymmetrical. 

The matrix element (41) has some other prop- 
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erties of interest. It vanishes unless the inte- 
grands are even functions of x =cos 6. Since Y,” 
is a polynomial in x of order .(m—m), the inte- 
grands are polynomials in x of the order 


Natngtn, 


By virtue of (40) this order is even or odd 
according to whether ”.+ +n, is odd or even. 
Hence: In the matrix elements of toroidal flow 
either two of the indices n are even and one is odd, 
or all three are odd. 

There are other selection rules for the indices 
n which we shall mention without going into 
details. These rules may be illustrated by the 
special case where Y(a) = P;(x) =x. It can then 
be shown from the recurrence relations of the 
spherical harmonics, and by taking account of 
the rule just stated, that the only non-vanishing 
matrix elements are those where »s=n,. Simi- 
larly, when for instance Y(a)=P2(x), the only 
non-vanishing elements are those where ng=n, 
+1. We need not elaborate on these and similar 
relations here. 


—M_—Mg—m,—1. 


Poloidal Flow 


We consider now the second type of the matrix 
elements (37) where the velocity vector is 
poloidal. On writing the element out in spherical 
coordinates we find by (13) and (14) 


[S(a)+S(8) x T*(y)] 
= =the thet 1)Gp f ¥(a)¥2¥(8)*w2¥*(y)do 
+(e +1)Ge f ¥(B)¥2¥(a)-v2¥*(v)da, (42) 


where 


G,=R? : Z Z d 
=f Fs (| ()Z(y)dr, (43) 


and similarly for Gs. Applying the identity (19) 
to the second integral in (42) it is found that the 
matrix element can be written 


[S(a)+S(8) x T*(y)]= —[ma(ta+1)Gp 
+-19(1g+1)Ga] f V(a)¥2¥(8)*w2¥*(y)de 


+ 1(thp +1) ty (tty +1) f V(a)¥(8)¥*(y)de. (44) 
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Here, everything is symmetrical with respect 
to the interchange of the indices 6 and y except 
the bracket multiplying the first integral over 
the angles. Some additional manipulation shows 
that neither the symmetrical nor the antisym- 
metrical part of this bracket is, in general, zero. 
Hence: For poloidal fluid motion the coupling 
matrix ts in general neither purely symmetrical nor 
purely antisymmetrical. 

Clearly, the selection rule (40) applies again. 
A further rule is obtained, as above, from the 
consideration that the integrals in (44) vanish 
unless the integrand is an even function of cos @. 
By a similar argument as before it is found that 
the integrand is even when n,+ng+n, is even. 
Hence: In the matrix elements of poloidal flow 
either two of the indices n are odd and one is even, 
or all three are even. 

Other selection rules for the indices m can be 
derived analogous to those given for toroidal 
flow. If for instance Y(a)=P;(x)=x, the only 
non-vanishing matrix elements are those where 
Np=n,+1. Again when Y(a)=P2(x) the non- 
vanishing matrix elements have either »g=n, or 
ng=n +2. We need not elaborate on those rules 


here. 
Scaloidal Flow 

The third type of matrix element, finally, 
where the velocity is a scaloidal vector is very 
similar to the type where the velocity is poloidal. 
The (6, ¢) functions of a scaloidal vector are 
identical with those of the corresponding poloidal 
vector, only the radial components are different. 
Hence the matrix elements of U(a) contain the 
same integrals over spherical harmonics as (42) 
or (44); only the radial integrals G, and Gz are 
replaced by somewhat different expressions. 
Again, the matrix elements are neither purely 
symmetrical nor purely antisymmetrical. The 
selection rules, depending only on the spherical 
harmonic functions, are the same as in the 
poloidal case. 

A table of values of some of the lower matrix 
elements will be found in Part II. 


INTEGRATION OF EQS. (34) 


The discussion will be limited to the case 
where the fluid motion is stationary so that the 
coupling matrix on the right-hand side of (34) 
reduces to a set of constants. 


WALTER M. 













ELSASSER 


Consider the system of differential equations 
dc,/dt+ >is M(B, y)cs=0. (45) 


The integration of this system offers no difficyl- 
ties when the matrix M is either symmetrical or 
antisymmetrical. In this case a suitable linear 
transformation of the c’s will make M diagonal, 
i.e., will lead to a system of new “normal modes” 
of the problem. It is furthermore known that in 
the Hermitian symmetry the transformation can 
always be assumed unitary. The unitary character 
may be expressed as a relationship between a 
transformation and its inverse, thus 


Bs= Diy Mayly, Cy = Des May * gs. (46) 


If by this substitution the system (45) is trans- 
formed to normal coordinates, the eigenvalues of 
the matrix M being \3, we have 


dgs/dt+ags=0, gs=ge(0) exp (—Ast). (47) 


One must now distinguish between the cases 
where the matrix M is symmetrical and where it 
is antisymmetrical. In the Hermitian symmetry 
there is associated to every antisymmetrical 
matrix a symmetrical one by the relation (25). 
The eigenvalues of a symmetrical (Hermitian) 
matrix are real; hence the eigenvalues of an 
antisymmetrical matrix are purely imaginary. 

It follows that in the case of an antisym- 
metrical matrix M in (45) the coefficients of the 
new normal modes are purely harmonic functions 
of time, so that 


| gs|?=const. 


It will be useful to derive a related result in a 
more direct way. If one multiplies (45) with c,* 
and adds to the formula obtained its conjugate 
complex there comes 


d/dt|c,|*= Yip M(B, y)cgcy* +p M*(B, y)cycp* 
and when M is antisymmetrical, i.e., 


M(8, 7) - —M*(y, B), 


one finds, by summing over y, that 
d/dt(>, | C,| 2]=0, (48) 


(convergence of all the bilinear forms involved 
being assumed). While this requirement is much 
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racter 
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weaker than the preceding one, it is also more 
general in that it does not presuppose a trans- 
formation to principal axes. 

The equations of motion (34) are of the form 


(45) with 


where 6(68,7)=1 for B=y and 6(8, y)=0 for 
#7; and where K is the coupling matrix dis- 
cussed in the last section. We have seen that 
when the velocity field is a toroidal vector (or a 
sum of toroidal vectors), K is an antisymmetrical 
matrix. If now the velocity is large, K becomes 
large and the term with A in (49) may be 
neglected in a first approximation. Then M is 
purely antisymmetrical and the results stated 
above apply. The amplitudes of the new normal 
modes carry out harmonic oscillations in time. 
If, from the complex, Hermitian representation of 
the problem we go to a representation where the 
fundamental vectors are real, the corresponding 
real amplitudes carry out sinusoidal oscillations. 

If now we admit a finite value for the A’s in 
(49), the matrix M is no longer either antisym- 
metrical or symmetrical. It is known that in 
this case the existence of a transformation of the 
matrix M to diagonal form, is not assured. In the 
case of a matrix without special symmetry it can 
therefore not, as a rule, be assumed that normal 
modes exist. The integration of the field equa- 
tions in a simple, closed form is then not possible 
and one must be content with less stringent 
results. One such result may be derived as 
follows. By the same procedure that previously 
led to Eq. (48) it is found that the bilinear form 
corresponding to the antisymmetrical matrix K 
vanishes, and so 


d/dt(dy|\¢y|*J=—Ly Ay|e,|*<0, (SO) 


since the A’s are essentially positive. Hence the 
over-all amplitude of the field decreases con- 
tinually. In physical terms the result may be 
stated as follows, assuming the A’s to be small: 
For toroidal flow the induction effect consists in 
oscillatory changes of the field amplitudes super- 
posed upon the slow, general decay of the field. 
We now consider the case where the coupling 
matrix K is symmetrical. In this case the matrix 
M of the Eqs. (45) is also symmetrical. The eigen- 
values \; of M are all real. Hence the coefficients 


of the new normal modes given by (47) are now 
real exponential functions of the time. The eigen- 
values may be either positive or negative. When 
the coupling matrix K is small, the eigenvalues 
will be near the ‘‘free’’ ones which are A and the 
exponents (47) will then all be negative. When, 
however, A is large the eigenvalues are near 
to those of K alone. Since all elements of K are 
proportional to the magnitude of v, all eigen- 
values change their sign when the direction of 
the velocity vector is reversed. Physically, the 
positive exponentials represent amplification and 
the negative exponentials de-amplification of the 
corresponding normal mode. Thus, the normal 
modes corresponding to negative \’s (positive 
exponents) increase without bound while the 
modes corresponding to positive \’s decay to 
zero. 

Unfortunately there is no simple type of fluid 
motion for which the coupling matrix is purely 
symmetrical. As we have seen, the matrix of 
poloidal flow contains both symmetrical and anti- 
symmetrical components. One can think of 
eliminating the antisymmetrical part by forming 
a linear combination of vectors of different 
types, but such a velocity field would have little 
physical significance and be of a very complicated 
geometrical form. 

An interesting illustration of the fact that the 
coupling matrix of poloidal flow is not simply 
symmetrical is found in a paper by Cowling.* 
This author investigates the question of whether 
the fundamental equations of induction, (32) or 
(34), can have stationary solutions. In the 
stationary case our problem would reduce to the 
set of equations 


Ly M(B, y)c, =0. 


This is a system of homogeneous algebraic 
equations for the coefficients c,. Expressed 
otherwise and in more physical terms, the 
problem is a boundary value problem for the 
field Eqs. (32). Cowlings presuppose that the 
fluid motion is incompressible, has rotational 
symmetry, and that the velocity vector is con- 
tained in the r—@ plane. In our terminology this 
corresponds to poloidal flow where only zonal 
harmonics are involved (m=0). Cowling has 


* T. G. Cowling, M. N. R. A. S. 94, 39 (1934). 
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given a formal proof to the effect that the 
problem in this form does not possess charac- 
teristic solutions. If the coupling matrix were 
symmetrical in the poloidal case, one might 
indeed be led to expect that characteristic solu- 
tions are possible. This goes to show that, even 
apart from the fact that fluid motion on the 
rotating earth is intrinsically turbulent, a sta- 
tionary model of induction of simple geometrical 
symmetry cannot exist, on purely electromagnetic 
grounds. 

When the coupling matrix is neither symmet- 
rical nor antisymmetrical, normal modes do not, 
as a rule, exist. The theory of integration of the 
differential equations then becomes far more com- 
plicated than in the case of symmetry. One would 
still expect that amplification or de-amplification 
of the field components is possible when the 
coupling matrix has a symmetrical part and 
when the velocities are large enough. If Msym 
designates the symmetrical part of M we have, 


WALTER M. ELSASSER 


in generalization of (48) and (50), 


d/dtlDy|¢y|*]= Ley Meym(B, v)cacy*. (51) 
It may be noted that the choice of the coefficients 
c in this relation is quite arbitrary, provided only 
that all the bilinear forms converge. Thus Eq, 
(51) applies to any magnetic fiéld. When Vv is large 
enough, M is proportional to v, and by the 
suitable choice of a numerical factor multiplying 
v the right-hand side of (51) can be made to 
have either positive or negative sign and any 
prescribed magnitude. Also, when the form (51) 
is sufficiently large at a given instant, a finite 
time interval must pass before it can change its 
sign. The results may be summarized as follows: 
To any given magnetic field at a given instant a 
fluid motion can be found which amplifies or de- 
amplifies this field at a prescribed rate, and con- 
tinues to amplify or de-amplify it over a finite 
length of time. 
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Relativistic Interaction of Electrons on Podolsky’s Generalized Electrodynamics 


, D. J. MonTGOMERY* 
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The wave equation for a system of particles is derived on the basis of Podolsky's generalized 
electrodynamics. An extension of some work of Fock leads to a representation in terms of a 
series of functionals. With this formalism the matrix element for the relativistic interaction of 
two electrons is determined, and is seen to be a generalization of Mgller’s formula. 





1. INTRODUCTION 


N aseries of papers Podolsky' has formulated the basis of a generalized electrodynamics involving 
higher derivatives in the field equations, and Podolsky and Kikuchi*® have developed the theory 
to include quantum electrodynamics. Here we extend the formalism, basing the treatment on the 
work of Fock,* and apply the results to the determination of the relativistic interaction of two 
electrons. 
2. WAVE EQUATION FOR A SYSTEM OF PARTICLES 


Derivation of Wave Equation for a System of Particles 
According to Dirac, Fock, and Podolsky,‘ and GE II and GE III, the Dirac wave equation for 


a system of particles and field, together with their interaction, is 


(H,+H;+H,,)¥ =ihav/aT, (2.1) 
with 


W=W(r,---r,; A(k), A*(k), o(k), o*(k); A(x), A*(k), (k), 6*(k); 7), (2.2) 
where 


H,=> (Cas-Pst+m,c*B,) ; (2.3) 


e=1 


A,= f [A* (lk) - A (Ic) —* (le) (lk) + (Ic) + A* (ke) — 9 (ke) *(Ie) Je2dk 


és f CA* (ke) - A (de) —3* (ie) (ke) +A (ke) -A* (kt) —$(e)$*(e) Jedke; (2.4) 


Hyy@> ef 6(t,, T)—e-A(t., T)]. (2.5) 


s=1 
When the single Eq. (2.1) with common time T is replaced by the set of equations with separate 
times ¢, in accordance with DFP, and after several transformations and the use of auxiliary con- 
ditions, it is shown in GE III that the equations to be solved are 


(ca,-P,’+m,c*B,)Q=T,'Q, (2.6) 
P,’ =p, —(e,/c)D(r,, t,) —(€,/2c)V.U,, (2.7) 
T,' =thd/dt, — (€./2c)dU,/dt, — (€,2/82a). (2.8) 


* Now at Princeton University, Princeton, New Jersey. 

1B. Podolsky, Phys. Rev. 62, 68 (1942). This paper has been called GE I. ; 

2B. Aa sa and C, Kikuchi, Phys. Rev. 65, 228 (1944); 67, 184 (1945). These papers will be called GE II and GE III, 
respectively. 

*'V. Fock, Physik. Zeits. Sowjetunion 6, 425 (1934). : . 

‘P. Dirac, V. Fock, and B. Podolsky, Physik. Zeits. Sowjetunion 2, 468 (1932). This paper will be called DFP. 
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where 
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The last two equations are GE III (4.2) and (4.3), with obvious minor changes in notation. The 
definition of U, is given by GE III (4.1): 


U.=T- (6/824): f [(1/k*) sin (vs— gs) — (1/84) sin (G.— ex) dk, 


where 
¢:=ckt,—K-r,, ¢:=ckt,—k-r,. 
In order to obtain the wave equation for the system, DFP shows that the equations for the 


individual particles are to be added, and the times are to be set equal to the common time 7. Now 


it is readily shown that 
0/8T =8/dt+ >, d/dt,; 


and inasmuch as & is independent of t, the effect of adding the equations and setting the times equal 
is to replace 0/dt, in the individual equations by 0/dT in the combined equation. The resultant 
equation is further simplified by showing that V,U,=0 when the times are set equal, and making use 
of GE III (4.3) and (4.12), namely: 


(1/c)8U,/ate= Da(eu/4x|1.—1u|)[1—exp (— |r—14| /a)] 


The resulting wave equation for a system of particles in the generalized quantum electrodynamics 
is thus® 


LD: {a.:[cp.—e.D(r,, ¢) ]+-m,.c?B,} Q= {ihd/dt— (1/8ra)>, €,” 
_ Li u(eseu/S| r,—r.|)[1—exp (—|r.—ru|/a)]}2. (2.9) 
Representation of Wave Equation in Functional Formalism‘ 


It will be convenient to transform the field variables to their Fourier amplitudes, by means of 
GE II (3.3): 


D(r,, t) =(1/2m)! f { D(k) exp [i(k-r,—kct) ]+D*(k) exp [—i(k-r,—ket) ]}dk 
+(1/2n)! f {Dde) exp [i(k-r,—kcet) ]+D*(k) exp [—i(k-r,—ket) ]}dk. 


From the commutation rules in GE III (3.7),’ it follows that D(k) may be represented by 
. 
D(k) = (ch/2k)* )) 6,(1/k)k Xe,b(k, j), (2.10) 
i=1 


where 87 =1, e; are a set of Cartesian unit base vectors, and the b(k, j) are operators satisfying 


“— Lo(k, 7), O*(k’, 7’) ]=6;75(k—k’); (2.11) 
and D(k) by 


~~ ~ 3 ~ ~ ~ 
D(k) = (ch/2k)! 2 8;(1/ak) (ak Xe;+e;)b(k, j), (2.12) 
j=1 


where 8;?=1, e; are defined previously, and 5(k, j) are operators satisfying 
[o(k, j); b*(k’, j= —6;/5(k—k’). (2.13) 


5 This equation was derived earlier by C. Kikuchi in different manner, but has not been previously published. 
* This is the formalism developed by Fock in the reference of footnote 3. : 
ad Note A typographical error in the second of the equations mentioned; the tilde over the &? in the right-hand side has 
n omitted. 
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For, according to (2.10) and (2.11), 
[D1(k), Dn* (k’) ] = (ch/2k)*(ch/2k’)* Ds, 57 BBs (1/k)(1/k’) (ke Xe;) -ex(k’ Xj) -CmLO(Kk, 7), O*(k’, 7’) ] 
= (ch/2k*) D5, 7 BiBy (kK Xe1-€;)(kK Xen -€5)5575(k—k’) 

= (ch/2k*)(k Xe,) - (kXe,,)5(k—k’) 

= (ch/2k*)(k-ke,-e,,—k-e.k-e,,)5(k—k’) 

= (ch/2k) (Sim —Rikm/k*)5(k—k’), 







which is the same as GE III (3.7). Analogously, using (2.12) and (2.13), we obtain 
[Di(k), Dn* (k’) ] = — (ch/2k) (51m —Rikm/k*)5(k—k’), 


the second of Eqs. GE III (3.7). 
Upon definition of 


G*(k, j) =(1/2m)'(ch/2k)* O°, €.8(1/k)a.-k Xe; exp [i(k-r,—ket) ], (2.14) 
G*(k, j)=(1/21)'(ch/2k)* D. €.8(1/ak)(a,.-ak Xe;+e,-e;)-exp [i(k-r,—ket)], (2.15) 








and 
H=>d., [a.-cp.t+m,c*B, |+(1/8ra) >, e2+Diu (€.¢,/4e|r.—1r.u|)[1—exp (—|r—ru|/a)]J, (2.16) 


the wave equation (2.9) becomes . 





HO—iha0/at= Sf aucerce, j6(6, ) +0, OO, i) 






+G*(k, f)b(k, 7) +G(k, j)b*(k, f)]}2. (2.17) 








The time factors in the exponentials in the G’s may be eliminated through transformations of type 


e~**tb(k, joe, 






where 





w=c> fawceorce, j’)b(k’, 7’) —k’b*(k’, 7’)b(k’, 7’)]. (2.18) 





From the commutation rules for 6 and b, it follows that 
e~twth (Ic, jet = (kk, jeri, e~*'5 (Ik, jet = b(kr, jrerickt, 
e~vth* (kk, jet =b*(k, jeniekt, e~*th* (Kk, jet =5* (kr, joe ie**, 






There is also the relationship 
e~**(thd /dt)e™! = —thd/dt+hw. (2.19) 





Upon definition of 





Go*(k, j) =G*(k, j)e**', and so on, 





the transformed wave equation becomes (where the transformed functional is designated by the 
same symbol as the original functional) 





Ef aulo* ce, j)6(6, j) He, HB Uc, NI} 


j=l 


(H—iha/a1)9+he| 





3 
1 f dk{Go* (i, j)b(K, {)+Go(k, j)b*(ke, #) +Go*(ke, be, 7) +Go(ke, f*de, f)]}9. (2.20) 


jel 
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Explicit Representation for Field Operators and Functional 


A sufficiently general form for th unctional is 


Q2= ia s Qye 
where 


Qn | 2d, i d, J fc, --dk,dl,- 1, -Wre(Ki, 4,°+-K,, tr; hi, jr L, js) 
-5(le;, i:) ++ -5(ke,, 4) (I, jx)+ +B (Ie, je). (2.21) 


Each sum is to be taken from 1 to 3 over the values of ¢ and j, and each integral over the entire 


momentum spaces of k and 1. 
The functional derivatives may be defined by 


b2(b(k, 7) ]/sb(k’, 7’) aesees (1/) {Q[6(k, j) +055;-5(k’ —k) ]—O[b(k, 7)]}, (2.22) 


where k represents the variables of integration and 7 the indices of summation in the functional. 
Definitions (2.21) and (2.22) permit the association 


b(k, i)~5/0b(k, i); b*(k, )~B(k, i); 5, j)~8/0b0, f); O*C, )~—4, A), — (2.23) 


for it may readily be shown that 


. (6/86(i, i))6(e’, i) 9—6 (kW, i’)(6/86(k, é))@=5,6(k—K)9 
= ~ ~ ~ ~ 
(—6/60(1, 7))O(1’, 7’) 2-41, 7’)(— 8/8001, 7)) Q= —4575(1—-1)2. 


These two equations are to be compared with the commutation rules (2.11) and (2.13). 


Application to Wave Equation 


The explicit representation developed in the preceding section may be applied to the wave equation 
(2.20) in order to obtain an ordinary wave equation in k-space. For the immediate purpose of this 
paper, we are interested in the case where the series of functionals is to be broken off after only the 
first three terms. (The technique which is to be used is, however, generalized readily to an arbitrary 


number of terms.) Then 
2 = Qo + Qio+ Qo, 


Qo0=VYo0, 


where 


2Q10= de dk wWi0(Ki, 11)6(Ki, 41), 


Q1= LX fatsdostts, jb(hi, ji). 
a 


The substitution of these expressions into the wave equation gives after some computation the 
following equations: 


(H—thd/dt)Poo= 205 f dk[Go* (k, j)~ro(k, 7) +Go*(k, j)¥or(k, 7), (2.25) 


(H+hek —thd/dt)ro(k, 7) =Go(K, 7) Yoo, (2.26) 
(H+hcek —ihd/dt)~or(k, 7) = —Golk, j) Yoo. (2.27) 





RELATIVISTIC INTERACTION OF ELECTRONS 


3. RELATIVISTIC INTERACTION OF TWO ELECTRONS 


Suppose the system consists only of two electrons and their field. To obtain a first-order approxi- 
imation (i.e., matrix elements proportional to the square of the electronic charge), it is possible to 
treat the last two terms in the definition of H (2.16), and the entire right-hand side of (2.25), as 
perturbations on an unperturbed Hamiltonian consisting of the first two terms of H. In order to 
eliminate Yo: and Yo from the right-hand side of (2.25), the two succeeding equations are solved for 
these two functions, with Yoo approximated by the wave function for the unperturbed system. The 
substitution of these results into (2.25) provides an equation amenable to standard methods of 


perturbation theory. 


Wave Function for the Unperturbed System 


The representative for two free electrons with momenta p,° and p,’, and signs of energy and direc- 


tion of spin designated by s;° and s,°, is given in fr, fe, £1, {2 space by 

. (ri, £1; Te, £2|Pr°, $1°; Po°, S2°) =exp (—iWt/h) go° . (3.1) 
where ; 
go° = (1/2rh)* exp [—i(pi°- r1+-p2®- re) /A Jus,s,(Pr°, 51°; po°, 52°). (3.2) 


Here u is antisymmetric, and has sixteen components corresponding to variables £1, {2, each of which 
has four values. The representative (3.1) is a solution of the wave equation 


(Fit Fr)y= WY, 
where 
F,=a,-cp.+m.c*B,, (3.3) 
and 
W°=W,°+ W.’, (3.4) 


with 
| Fi 90° = Wi go", ~Feeo°= Wego". (3.5) 


Elimination of vo and tn 


Let us define ’ 
, YVi=f=f exp (—1W%t/h), You=g=g° exp (—i1W%t/h). (3.6) 


It is clear from (2.26) and (2.27) that f and g have the same time dependence as (3.1) when yoo has 
been approximated by this wave function. Hence f* and g® are independent of time, and (2.26) and 


(2.27) lead to , 
(Fit Fethck — W°) f° =Go(k, 7) ¢0°, (3.7) 


(Fit F.+hck —W°)g° = —Go(k, 7) ¢0°, (3.8) 


where the static interaction and the self-energy have been neglected, as they are proportional to the 
square of the charge, and would give terms proportional to powers higher than second in the final 
interaction matrix element. 

The solutions for f° and g® are 


f=Og0, g°=Og0", (3.9) 
where 


6=(1/2m)'(hc/2k)!- { (Fithck —Wo')—68;(1/k)(e1-k Xe,) exp (—ik- 11) 
. +(Fot+hck — Wo?)—68 ;(1/k)(a2-kXe;) exp (—ik-re)}, (3.10) 

6=(1/2m)"(hc/2k)*- { (Fithck —Wo')—€:8,(1/ak)(e1-ak Xe;+a1-e;) exp (—ik-r:) 
+ (Fethck — Wo?)-'28j(1/ak)(a2-ak Xe;+e2-e;) exp (—ik-re)}. (3.11) 


Here we have made use of the fact that a; and a2: commute, and that ¢»° satisfies (3.5). 
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The wave equation (2.25) becomes 


Now y differs from yoo only by a perturbation contribution ; and the operator preceding y is itself a 
perturbation operator. Hence we may replace y by yoo, and the equation is in standard form for 
application of perturbation theory. 


Calculation of Interaction Matrix Element 
The perturbing energies are, for two particles of charge €, and é2, 
€1€2[ 1 —exp (-— |ri1—fe| /a) ]/4e|r1—12| =U,, (3.13) 
(62+ €2*) /8ra= Vi, (3.14) 


= fae S (Gut0—Ge") = Ut Vo 3.15) 


where U; represents the part of the left-hand side of (3.15) containing terms in €,€2, and V2 represents 
the part containing terms in ¢€,? and ¢?. Since we are interested in the interaction only, we shall 
calculate only the matrix elements for U; and U:. Further we assume conservation of energy: 
Wi+ W2=W;°+ W.». It is well known that the matrix element 


(pi, 51; D2, S2| €:€2/42 | r1—Tf2! |p.°, $2°; P2", $2°) => (1/22) %e,€26(Pi+P2)(1/hP 1?) (u*, u°), (3.16) 


where P,=p,—p,°, and 


(u*, u°) => us,t,*(Pi, $1; Pa, S2)us,¢,(P1°, $1°; p2’, S2°). 
15s 


From the — ee: exp (—|r1—r2|/a)/4r|r1—r2| term, the contribution to the interaction matrix 
element is calculated to be the negative of (3.16), with P;? replaced by P:?+h?/a?. That is, 


(Pi, 51; Po, S2| — €1€2 Exp (- |r: —re| /a)/4e|r1—Tre| | pi°, 51°; p2°, $2°) 
farses" ee exp (- |t1—fe| /a)/4e|11—Ta| } go° 


= — €1:€2(u*, w(1/2nh)* ff drs: exp (- |r1—fe| /a) ‘exp { —il(pi—p.°) -(f11—T2+12) 
+ (p2—p2°) -t2}/h}/|t1—T2!. 


Define P,=p:—p,"’, P2=p2—p.°, R=r,—r2; the integral immediately above becomes 


Jf faedR exp {iLP-R+ P+) -t2]/A) (O™*/R) 
= (2nh)*5(P,-+Ps)(4h/P,)(P1/h)/C(1/a2) +(P12/h?)] 


Then the matrix element for U; is finally : 
— (1/2m)*e1¢25(Pi+P2)(1/h)(P1?+h?/a*)-'(u*, u°). 
For the part of — {dk + Go*@ which contains €¢2, the contribution is 
—(1/2m)*e¢25(Pi+Pe2) {AL Pi? — (Wi— W1)?/c?}}—*- (u*, (a1 @2—@1*Pyae-P,/Pi*)u°); (3.17) 
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for the part of + /dk >) Go*@ which contains ¢:¢2, the contribution turns out to be the negative of 
(3.17) with P replaced by P;?+h?/a*. The proof follows. 
From the appropriate definitions ((2.14), following (2.19), and (3.10)), we have 


ro Grt=L (1/2m)*(hc/2k*) - { €:2a,-k Xe; exp (—ik-r:)(Fi +hck— W)")'e)-k Xe; exp (—ik-r)) 
+ €,°a2:k Xe; exp (—1k-1,)(F.+hck— W2")-'a2-k Xe; exp (—ik- re) 
+ €,e2a2-k Xe Fi +hck — W,°)-'a,-k Xe; exp [7k- (re—1) ] 
+ €1€20;-k Xe;(F2+hck — W2")-'a2-k Xe; exp [tk- (r2—1) }}. 


We are interested only in the terms in €,¢2:; the matrix element for the first of these is, since F; is . 
Hermitian, 


> 5 7. J ce((Fthck— W,°)—) go*ae-k Xe ja; -kXe; exp [ —ik- (r2—11) ]¢0° 
where the second sum is to be taken over the spin variables, and the integral over the r, fz variables. 


The angular brackets indicate, of course, the Hermitian conjugate. With the help of (3.5), the ex- 
pression becomes 


€1€2 Ef (Withek— W,°)—! go* (a1 Xk) ° (a2 Xk) exp [ik- (r2—11) ]¢o° 
= ase f (Withee _ W,°)"! exp [— i(pi . ri+pe2 ° r2)/h | 


-(u*, (a> a2k? — a: kKee-k) exp [ik-(re—r,) ] exp [4(pi°- ri +pe- re) /h ju) 


The complete expression for the term without tildes is 


—(1/2m)3(hc/2k) eres f dk f dr,dr2(Wi+hck—W,°)--exp [—i(P,—hk) -1,/h] 
“exp [ —i(P2+hk) . ro/h \(u*, (a; *@o— Q@ -kae -k/k?)u°) ° 


Integration with respect to r,; produces a 6(P,—/k) factor, and integration with respect to Ak 
replaces the Ak by P,. The final integration with respect to re gives a 6(P2+P,) factor. The result is 


(1/29) e,€0(c/h)6(Pi1 +P.2)(Wi— Wi°+c!|pi—pi®|)-'- (2! p:—py®| 3) ~1 
X (u*, (a+ @2—@,*Pya2-P,/P,*)u°). 


For the second term containing €,€2, the subscripts 1 and 2 are interchanged. Upon use of the 
expression for conservation of energy, we find the sum of the two to be (3.17). 


The calculations for + (dk >> Go*@ are of the same type. 
The part containing the Dirac matrices can be simplified ‘further ; for by (3.2) and (3.3) 


F,go° = (1/2rh)* exp (ipe®- re/h) -Ler-cpi exp (ipi®-11/h) +myc*B, exp (ipi®-1r,/h) Ju® 
=(1/27h)* exp (ipe®-re/h)-Lexp (ip.°-11/h)a.-cpiu® 
+a;-cp,° exp (ip,?-r:/h)u°+my,c*B, exp (ip.®-1r,/h)u®} 
= (1/2mh)§ exp [i(pi°-r1t+-pe®- rz) /h (Fite: cpi")u® 


_ Wi" v0", 








‘ 
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whence 
(Fite . cpi°)u° = Wu; 







and since, in the expression (3.3) for F,, 


a1: Cpiuv° =a, ° (h/1)du°/dr, =(, 






it follows that e-px°w?=(1/c)(Wi°—my,c*B;)u°. Analogous relations hold for ai-pi, @2-P2°, ae-py, 
Now 






(u*, —@1-Pya2-P,u°) = (u*, [e1-(Pi—pi°) a2: (Pp2—pe”) |v), 





by definition of P, and application of the principle of conservation of momentum in accordance 
with the 6(P:+P:) factor in the matrix element. The expansion of the quantity in brackets gives 







@1*Pi@2*P2— @1°Pi1°a@2-P2— @1-Piae* Po’ + a1-Pi°a@2: po”. 







Then 
c?(u*, @1° Dias: Pov”) = (La2-Ppoay- piu |*, u°)c? 
_ (Lee: po(Wi—myc*B;)u }*, u°)c, 
= ([(W2—mec*B2)(Wi—mic*B,)u }*, 4°), 
=(u*, (We—myec*Be)(Wi—myc*B,)u°) ; 
and 






—c?(u*, a1-pi°@e- pow) = — (Lae-pou ]*, a1-piu®)c?, 
= — (L(We—mec*Be)u |*, (W1°—my\c*B,)u°), 
= —(u*, (We—me2c*B2)(W1°— myc*B;)u°). 
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Analogous expressions are computed in similar fashion for the other two terms in the expansion. 
Upon combining the four terms, the §’s disappear, and the result is 


(u*, (Wi— W,1°)(W2— W2°)u). 












By conservation of energy, Wi1— W1°= —(W2— W,°), whence 
(u*, — a+ Pya2-Pyu°) = —(1/c?)(u*, (Wi— W,°)*u). 


=~. ame 





“ “pe 





Upon combining our various results, we get for the interaction matrix element the following: 
(1/2m)*(1/h)5(pit+p2—pi°—pe") { |pi—Pi®|°— (Wi — W")?/c?} 
*{1+(@?/h*)||pi—pi?|*— (Wi— Wi°)?/c?}}-*(u*, (1 — 1-2) 0°). 


This result is a generalization of Mgller’s formula.® It will be noticed that it is relativistically in- 
variant, and reduces to Mgller’s expression as a—0. 
The author is grateful to Dr. Boris Podolsky for his continued encouragement and assistance. 











(3.18) 










*C. Moller, Zeits. f. Physik 70, 786 (1931). - 
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Combination of Betatron and Synchrotron for 
Electron Acceleration 


H. C. PoLitock 
Research Laboratory, General Electric Company, Schenectady, New York 
November 13, 1945 


N a recent letter E. M. McMillan" discusses a possible 

design for a 300-Mev electron accelerator known as a 
synchrotron. If electrons are injected into the synchrotron 
with 300-kev energy, the radius of their orbit varies from 
78 cm to 100 cm as their velocity gradually approaches that 
of light. Thus the magnetic field must cover a ring of width 
greater than 22 cm in order to shape the field properly. 

It should be possible to reduce the width of this magnetic 
flux ring by using the betatron principle of acceleration 
until the electron energy has reached 2 or 3 Mev. Electrons 
could be injected into a betatron** orbit by a small gun 
pulsed to 50 or 100 kv. They would be further accelerated 
on this orbit of radius slightly less than 100 cm by in- 
creasing the flux in a relatively small central magnetic core, 
designed to saturate when the electron energy increased to 
some appropriate value over 2 Mev. Before the core satu- 
rated, R.F. accelerating voltage would be applied with 
proper timing to the gap in the toroidal vacuum tube. 
Thereafter the energy gain of the electrons would come 
almost entirely from passage across the accelerating gap. 
As the equilibrium energy of the particles increased, the 
proper energy change per cycle would be controlled by the 
phase stability which characterizes the synchrotron. 

The combination of betatron acceleration at low energies 
with synchrotron acceleration at high energies, where 
radiative losses become serious, should allow a saving of 
considerable magnet weight in machines designed for the 
billion-volt range. To test the design, a 70 Mev accelerator 
of this type is being constructed at the General Electric 
Research Laboratory at Schenectady, New York. 

It is understood that Powell and Bohm at Berkeley 
have noticed independently the desirable features of this 
form of electron injection and that it will probably be em- 
ployed in the 300 Mev machine proposed in McMillan's 
letter. 
author bef sy Rev. 68, 143 (1945) (communicated to the 


D. W. Kerst, Rev. Sci. Inst. 13, 387 (1942). 
* W. F. Westendorp and E. E. Charlton, J. App. Phys. 16, 581 (1945). 
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A Theory of the Low Mass Neutral Mesotron 


Curairo KIKUCHI 
Department of Physics, Michigan State College, East Lansing, Michigan 
February 1, 1946 


T has been shown by F. Bopp" and B. Podolsky,? if the 
Lagrangian of the electromagnetic field depends upon 
both the field strengths and their first derivatives, (a) that 
the electrostatic self-energy, which diverges for a Maxwel- 
lian field, is finite; (b) that the field potentials, gq, satisfy 
a fourth-order equation of the form 


(1-a°()) D¢0 (1) 
with the static solution 
e=(e/r)(1—e-"*); (2) 


and (c) that the field is a superposition of fields of quanta 
of rest masses 0 and h/ac. 

The purpose of the present note is to pojnt out a possible 
relationship between this new field theory and the low 
mass neutral mesotrons reported by G. Groetzinger, P. G. 
Kruger, and L. Smith.* 

Several values for the constant a have already been 
proposed. Bopp has given a=r0/2, which leads to 2(137) 
mo for the mass of the heavy particle, where ro is the class- 
ical electron radius and mo, the electron mass. Landé,‘ 
on the other hand, proposed the value #ro for a, giving 
3(137)mo for the neutral mesotron mass. However, it 
occurred to the author some time ago that there is still 
another possible value for a. If the electron-proton inter- 
action potential in the hydrogen atom is assumed to be 
given by Eq. (2), the perturbation energy of the 2S-level 
due to the exponential term is found to be 


AE~a'?é/2ae 
or, in wave numbers, 


Av(2S)~a?é/2ag’ch = (a*a?/2r¢*) (mo c/h), 


where « is the electronic charge, ao, the Bohr radius of the 
hydrogen atom, and Avy the wave number shift. Using 
Pasternack’s value of 0.03 cm™ for the 2.S-level displace- 
ment, we find that a~10ro, giving 


m=h/ac~14mo. 


This agrees remarkably well with the value, 16m», quoted 
by G. Groetzinger and L. Smith.* 

The use of the potential given in Eq. (2) for the elec- 
tron-proton interaction is not free from objections, for, 
according to the well-known argument of Pasternack,’ the 
electron-proton force must become repulsive for small 
distances. In the result published recently by A. Pais,* 
this difficulty is avoided by using 


eg = —(e/r)(1—2e-"*) 


from which he obtains 2.579 for a. This value, however, 
gives 55m for the neutral mesotron mass. 


1F. Bopp, Ann. d. Physik 38, 345 (1940). 

2 B. Podolsky, Phys. Rev. 62, 68 (roe) 
( oan Groetzinger, P. G. Kruger, and L. Smith, Phys. Rev. 67, 52 
1945). 

4A, Landé, Phys. Rev. 60, 121 (1941). 

’C. Kikuchi, Thesis, University of Washington (1944). 

* G. Groetzinger and L. Smith, Phys. Rev. 68, 55 (1945). 

7S. Pasternack, Phys. Rev. 54, 1113 (1938). 

8 A. Pais, Phys. Rev. 68, 227 (1945). 
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Transfer Phenomena and Indeterminacy 


J. G. DAUNT AND K. MENDELSSOHN 
Clarendon Laboratory, Oxford, England 
December 21, 1945 


FEW years ago we advanced the hypothesis of a 
basic analogy between superconductivity and the 
\-phenomenon in liquid helium.' Its salient feature is the 
fact that superconductive electrons and superfluid helium 
atoms remain energetically at absolute zero, even when 
the substance as a whole is raised to a finite temperature. 
Thus, in both cases ‘frictionless transport is provided by 
an assembly of particles with zero entropy (z particles). 
Experimental evidence for this conclusion is furnished by 
our observation of zero Thomson heat in a superconductor* 
and of a mechano-caloric effect in liquid helium II.* The 
latter experiments, which had to be discontinued in 1939, 
have since been confirmed with greater accuracy by Ka- 
pitza.* The analogy is further strengthened. by our dis- 
covery of a critical rate of transfer in helium® which is the 
counterpart of the “current threshold” in superconductors. 
It has been suggested by one of us® that the momentum 
of frictionless transfer is derived from zero-point energy 
and that the rate of transfer of mass or charge (R) gives an 
indication of the number of z particles (m). This view seems 
to be supported by an important paper of F. London.’ 
Taking up our analogy, he has shown that, numerically, 


Anh/2 (1) 


in superconductors as well as in liquid helium. While in 
our opinion the experimental data hardly justify an evalua- 
tion of m to a greater accuracy than an order of magnitude, 
there seems to be little doubt that R/n is of the order of h. 
London’s relation can also be written as 


f mo-dx~h, (2) 


where m and v are the mass and the velocity of the z par- 
ticle and where the integration is carried out over the cross 
section of the specimen along which transfer takes place. 
Taking v as the average velocity throughout the cross 
section, we get 


o~h/md, (3) 


where d is the thickness of the specimen carrying the 
transfer. Thus, v is identical with the zero-point velocity 
conferred by the uncertainty principle on a particle con- 
fined to a space Of linear dimension d. By assuming friction- 
less transport to be zero-point motion, London’s formula 
(Eq. (1)) is then seen to follow as a necessary corollary. 

The peculiar result that the velocity of the z particles is 
inversely proportional to the linear dimensions of the 
specimen is well borne out by Silsbee’s rule in supercon- 
ductors and by the observation of Allen and Misener® that 
frictionless flow of helium is proportional to the diameter 
of the capillary. These effects have usually been regarded 
as surface currents but Eq. (3) indicates a much deeper 
physical significance. Frictionless transport clearly implies 
that, if collisions take place, not only the sum of momenta 
is conserved. but also that each individual momentum is 
conserved, regardless of the identity of the particles carry- 
ing it. Individual momenta will therefore be conserved 
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until a collision with the boundary of the specimen takes 
place. Thus, the observed dependence of v on 1/d is a 
direct consequence of the uncertainty principle and arises 
from the impossibility of locating a z particle in the speci. 
men with greater accuracy than d. 

This fundamental inability to locate z particles within 
the macroscopic dimensions of the specimen causes fric- 
tionless transport to appear in all determinations as a 
surface flow. (In multiply connected bodies the situation 
is more complex, though not fundamentally different.) 
However, what significance has to be ascribed to the 
“‘depth of penetration”’ of this surface flow is by no means 
clear, and recent observations by Casimir? seem to show 
that its physical reality has been interpreted too literally 
up to now. 

+ ¢ . Daunt and K. Mendelssohn, Nature 150, 604 (1942). 

. o Daunt 4 K. Mendelssohn, Nature 141, 116 (1938); Proc. 
Roy. “Soc. A, in p 

. G. Daunt a fc. Mendelssohn, Nature 143, 719 (1939). 

. L. Kapitza, Phys. Rev. 60, 354 (1941). 

J. 6 Daunt and K. Mendelssohn, Proc. Roy. Soc. A170, 423 (1939). 

* K. Mendelssohn, Proc. Phys. Soc. 57, td = ). 

7 F. London, Rev. Mod. Phys. 17, 3c Oe 


on and Misener, Proc. Roy. ‘Alz72, Xr (1939). 
B. G. Casimir, Physica 7, ‘sere (1940). 





A Proposed High Energy Particle 
Accelerator—The Cavitron 


RICHARD F. Post 
Naval Research Laboratory, Washington, D. C. 
January 9, 1946 ’ 


S tools for the investigation of nuclear phenomena, 
particle accelerators are unequaled. The cyclotron, 
the betatron, and most recently, the synchrotron"? are the 
best known examples. In the quest for higher and yet higher 
particle energy outputs, accelerators have, perforce, been 
made larger and larger, until their very size precludes 
construction by any but the largest laboratories. It is the 
purpose of this letter to explore the potentialities of an 
accelerator of much different physical design from any 
of the above, compact in construction, and in principle 
capable of producing much higher peak energies than are 
at present available. 

Little has been published on the advantages of an ac- 
celerator which would incorporate no iron in the magnetic 
circuit. Terletzky' points out that, if radiation effects could 
be eliminated, an iron-less betatron should be capable 
of producing very high energies. Whereas the phenomenon 
of radiative dissipation is a serious limitation upon the 
ultimate attainable output from a betatron,*® it has been 
shown" ® that the mechanism of “resonance” acceleration 
(as employed in the synchrotron) is far less subject to this 
effect. In addition, this means of particle acceleration is 
not limited by relativistic effects. 

The accelerator herein proposed is to be iron-less. 
It is to employ “resonance” acceleration, achieved by slow 
increase of the guiding magnetic field. High intensity 
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(1939), 








electric accelerating fields are to be obtained by replacing 

the usual “dee”’ assembly and acceleration chamber by a 
cavity resonator, similar to types now commonly employed 
in radar. The use of such a cavity becomes practical only 
in an iron-less device, because of the limitation imposed by 
the air-gap in present accelerators. 

The maximum orbit radius is to be small. This requires 

a large guiding magnetic field, but at the same time greatly 
reduces the volume through which this field must be 
maintained. It, therefore, appears practicable to immerse 
the entire accelerator in a liquified gas, such as hydrogen or 
nitrogen. Cooling the assembly materially reduces ohmic 
resistance losses in the field windings, and thus simplifies 
the problem of passing large current pulses. It has also the 
advantage of increasing the “Q”’ of the resonator. 

Although many of the commoner types of cavity resona- 
tors’ exhibit electric field configurations adapted to pro- 
ducing resonance acceleration, the type describable as a 
“sphere (or spheroid) and reentrant cones’’* seems particu- 
larly suitable. In such a cavity, the electric field lines ap- 
pear as arcs of circles, terminating on the surfaces of the 
reentrant cones. Small slots are to be made at points in the 
conical surfaces, through which the accelerated particles 
are to pass in their circular orbits. Injection may con- 
veniently be accomplished from a point external to the 
accelerating fields by placing the injector gun within one 
of the cones. A target may also be conveniently located 
within the cones. It will probably be necessary to partially 
“open-circuit” the cavity walls for low frequency currents 
generated by the guide field. This can readily be done. 

Within a resonator such as described, it should be pos- 
sible to obtain exceptionally high accelerating potential 
drops, of the order of several hundred thousand volts. 
This follows from the fact that a cavity resonator is an 
extremely efficient device, often exhibiting an electrical 
“QO” greater than 100,000 and consequently requiring a 
relatively small energy input to produce large internal 
fields. The upper limit of attainable electric fields within 
an evacuated cavity should be very high, and will probably 
be imposed by field emission effects. In the use of such a 
resonator, it is convenient to employ the cavity itself as 
the frequency controlling element of a power oscillator, 
thus eliminating problems of accurate frequency control. 

Some small advantage may also be gained from enclosing 
the particle orbits within a closed conducting surface, in 
that this will minimize that part of the radiation losses 
contributed by low order harmonics. It has been shown,®* 
however, that these harmonics contribute only slightly to 
the total radiative dissipation. 

The conditions for “resonant” increase of energy’? can 
be satisfied at all times during an acceleration cycle, 
provided the amplitude of the resonance electric field is 
increased simultaneously with the rise of the magnetic 
field. It can be shown that under these conditions the effect 
of the magnetic component of the cavity field will be 
negligible throughout the cycle, provided adequate fo- 
cusing forces are established by the guide field. 

To provide the necessary pulses of current through the 
field windings, three methods could be employed: (a) a 
“pulse transformer’ excited from a d.c. source, (b) a 
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short-circuited d.c. generator (following Kapitza’s tech- 
nique), or (c) a bank of high capacity storage batteries of 
the type used in submarines. A large ignitron and delay- 
line extinguishing circuit would serve as the switching 
element. The limits on the particle yield of the accelerator 
will be set by the available power and by the allowable 
rate of heat dissipation from the field windings. 

Some design calculations have been made for a 500-Mev 
electron accelerator. 

The author wishes to express his appreciation for counsel 
and suggestions given by Mr. H. F. Kaiser, Dr. Ross Gunn, 
and Mr. D. dePackh of the Naval Research Laboratory. 

1 Veksler, J. Phys. Acad. Sci. USSR 9, SS i» Sese8). 
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3 Terletzky, J. P’ Acad. Sci. USSR 9. 159 (1945S). 

*D. Iwanenk FART Pomeranchuk, Phys. Rev. 65, 343 
po. ——— and Pomeranchuk, J. Phys. Acad. Sci. SR 5, 267 
+E. M. McMillan, Phys, Rev. 68, 144 (1945). 

W. Hansen and R. D. Richtmyer, 4 p. Phys. 10, 189 (1939). 
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New York, 1943), p. 288. 
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Nuclear Induction 


F. Biocu, W. W. HANSEN, AND MARTIN PACKARD 
Stanford University, Stanford University, California 
January 29, 1946 


HE nuclear magnetic moments of a substance in a 
constant magnetic field would be expected to give 
rise to a small paramagnetic polarization, provided thermal 
equilibrium be established, or at least approached. By 
superposing on the constant field (z direction) an oscillating 
magnetic field in the x direction, the polarization, originally 
parallel to the constant field, will be forced to precess 
about that field with a latitude which decreases as the 
frequency of the oscillating field approaches the Larmor 
frequency. For frequencies near this magnetic resonance 
frequency one can, therefore, expect an oscillating induced 
voltage in a pick-up coil with axis parallel to the y direction. 
Simple calculation shows that with reasonable apparatus 
dimensions the signal power from the pick-up coil will be 
substantially larger than the thermal noise power in a 
practicable frequency band. 

We have established this new effect using water at room 
temperature and observing the signal induced in a coil by 
the rotation of the proton moments. In some of the experi- 
ments paramagnetic catalysts were used to accelerate the 
establishment of thermal equilibrium. 

By use of conventional radio techniques the induced 
voltage was observed to produce the expected pattern on an 
oscillograph screen. Measurements at two frequencies » 
showed the effect to occur at values H of the z field such 
that the ratio H/» had the same value. Within our experi- 
mental error this ratio agreed with the g value for protons, 
as determined by Kellogg, Rabi, Ramsey, and Zacharias." 

We have thought of various investigations in which this 
effect can be used fruitfully. A detailed account will be 
published in the near future. 


1J. M. B. Kello I. I. Rabi, N. F.{Ramsey, and J. R. Zacharias, 
Phys. Rev. 56, 738 193 9). 
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Concerning the Anomalous Scattering 
of Mesotrons 


R. P. SHutt 
Bartol Research Foundation, Swarthmore, Pennsylvania 
January 14, 1946 


AKING use of a method developed by the writer! 

Sinha? has recently calculated the cross section for 

the anomalous scattering of mesotrons within a band of low 

energies. In this method the energies of the individual 

scattered particles are not determined but the distributions 

in angle of the particles scattered in two lead plates of 
different thicknesses are compared. 

As Sinha has pointed out, the writer’s earlier calculation 
of the number of singly, or anomalously, scattered par- 
ticles was not strictly accurate since in this calculation an 
approximation was required because of the fact that when 
expressed as a function of u=6f~+ the distributions ob- 
served with two thicknesses of lead had different upper and 
lower limits. When this calculation was published it was 
assumed that the error thus introduced was small com- 
pared with the statistical errors of the observations. A more 
rigorous analysis of the problem recently undertaken in 
connection with the reduction of a much larger mass of 
experimental data has revealed that even in these pre- 
liminary results the estimates of the anomalous scattering 
were significantly affected by errors introduced in the 
approximations and the calculated cross sections for the 
anomalous scattering were consequently too small. The 
improved method of analysis will be published soon along 
with a report on the analysis of additional experimental 
data. 

There remain several discrepancies between the pro- 
cedures employed by Sinha and by the writer. The values 
for f(u) given in Table II of Sinha’s article should be multi- 
plied by a factor of 4. This leads to a value for the quantity 
6 of 30.2 percent. Furthermore, it has not been taken into 
account that 6 represents the difference between two 
integrals of f(u) over the same range of u but over different 
ranges of scattered angles. Therefore the writer believes 
that Sinha’s method of obtaining the quantity K from 3 is 
incorrect. Using Sinha’s data the writer has re-calculated a 
scattering cross section, and, because of cancellation of the 
several errors mentioned, accidental agreement is found for 
this and Sinha’s original values except for the amount of 
the standard error which turns out to be almost twice as 
large as that given by Sinha. 

The new results following from Sinha’s and the writer's 
experiments are given in Table I. Results obtained from 
an earlier article by Code* have also been included for 
comparison. Lower and upper limits for the observed 
ranges of the scattering deflections and energies of the 
particles are indicated in columns 3 and 4. Next, the cross 
sections for the anomalous scattering are listed. These 
cross sections represent averages over the ranges given and 
refer to one neutron or proton in lead or tungsten. It is 
seen that satisfactory agreement exists for Code’s and the 
writer’s last values for the scattering cross sections, the 

only two values belonging to nearly equal ranges for 
deflections and energies which may be compared directly. 
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TABLE I. Cross sections for anomalous scattering. 











1 2 3 4 5 6 
Limits of Fraction 
eof Limits of of total, 
defiec- range of Average scattered 
tions energies cross anoma- 

Author Method involved involved section lously 
10-28 cm? 
per 
Degrees 108 ev nucleon % 





Sinha 2 lead 4 +24 0.554155 170475 65430 
plates, 2 
and 4 cm 


thick, resp. 





Code 1 tungsten 9 +90 6 +2 622 80+20 
plate, 3.8 
cm thick, 
energy 
measured 
by magn. 
field 





+o 16+2 55+10 
20+3 88+12 


+o 5+2 81+19 


Shutt 2 lead 4.5 +27 
plates, 1 4.5 +27 
and 5 cm 9 +27 


thick,resp. 9 +27 


POPS 

uiviniv 
I 
8 








Sinha’s value is (10+5) times as large as the writer's 
first value for a similar range of scattering angles. As 
Sinha concludes, the reason for this discrepancy may lie 
in the difference of the energy ranges within which the 
observations were made. 

In the sixth column the fraction of the total number of 
deflected particles is given for the number of mesotrons 
that are scattered anomalously. It appears that approxi- 
mately 50 percent of the mesotrons are scatteted anoma- 
lously if low energies are included while for energies greater 
than ~5X10* ev the scattering seems to be almost ex- 
clusively anomalous. 

1R. P. Shutt, Phys. Rev. 61, 6 (1942). 


2M. S. Sinha, Phys. Rev. 68, 153 (1945). 
3F. L. Code, Phys. Rev. 59, 229 (1941). 





A Theoretical Criterion for the 
Initiation of Slip Bands 


CLARENCE ZENER 
Institute for the Study of Metals, University of Chicago, Chicago, Illinois 
January 31, 1946 


S was first recognized at the end of the last century,' 
initial plastic deformation occurs in most metals 
through the propagation of slip bands. It is therefore 
reasonable to expect that in these metals the over-all 
stress necessary to initiate plastic deformation may be 
computed by the condition that such stress be sufficient to 
propagate slip bands. The purpose of this letter is to 
present the computation for the yield stress of metals 
according to this viewpoint, and to demonstrate that the 
results appear to be in agreement with observations. 
From observations upon the mechanical behavior at 
low stress levels of freshly deformed metals,? it appears 
that the material within a freshly formed slip band behaves 
in a viscous manner,* and that the rate at which shear 
stress is relieved decreases rapidly with the time since its 
formation.‘ It is therefore reasonable to assume that during 
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Fic. 1. Comparison of Eq. (1) with data on silicon brass in reference 7. 
Cross-hatched region observed “‘apparent elastic limit."’ Dotted line, 
Eq. (1) with C =1.8. 


its propagation, which is too rapid to follow visually, all 
shear stress is relieved across a slip band. Such relief of 
shear stress within a slip band necessarily results in a 
large concentration of shear stress just ahead of its ad- 
vancing edge, hereafter called its spearhead. The condition 
that a slip band propagate is that this stress concentration 
be sufficient to induce plastic deformation in front of the 
spearhead. 

Under the above assumptions the stress concentration 
factor is given by® stress concentration factor —(L/p)', 
where L is the lateral diameter of the slip band and p is 
the radius of curvature of the spearhead. The local shear 
stress necessary to induce plastic deformation in a pre- 
viously undisturbed lattice is given by® 


local shear stress<°G, 


where G is the shear modulus. The macroshear stress 
necessary for slip band propagation is therefore 


macro shear stress—(p/L)IG. 


The advancing spearhead will probably be as thin as 
is consistent with the atomic structure of the metal. p will 
therefore be replaced by the spacing d of the atomic planes 
parallel to the slip band. On the other hand, slip bands are 
confined to individual grains, and therefore L is less than 
a grain diameter. The yield stress, which is effectively the 
least tensile stress which induces plastic deformation, is, 
according to the above viewpoint, therefore given by 


yield stress = C(d/D)!-G, (1) 


where D is any convenient measure of the average grain 
diameter, e.g., the standard “grain size,” and C is a numeri- 
cal coefficient of the order of magnitude of unity. A 
comparison of Eq. (1) with the extensive data’ on copper 
alloys indicates the essential correctness of the above 
approach. Figure 1 is presented as an example of the type 
of agreement obtained. 

According to the above viewpoint, a slip band must 
already have attained a critical length before it can propa- 
gate under a given macro stress. The same type of difficulty 
is encountered in a description of a change of phase—a 
particle of a new phase can grow only after it has already 
attained a critical size. The nucleus of a slip band presum- 
ably forms, as does the nucleus of a new phase, by some 
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sort of thermal fluctuation, perhaps of dislocations. The 
comparatively slight effect of changes in temperature and 
in strain rate upon the observed yield stress probably 
enters through their effect upon the rate of formation of 
new nuclei. 

Resistance to deformation by twinning is of the same 
order of magnitude as is the resistance to deformation by 
slip bands. The spearhead of a twin and of a slip band must 
therefore propagate in essentially the same manner. 
Depending upon the plane of the spearhead, the material 
further back deforms either in an orderly manner, thereby 
resulting in a twin band, or in a disorderly manner, thereby 
resulting in a slip band. 

Except in metals containing obvious gross cavities, 
such as cast iron, and in metals of complex structure, 
fracture is always preceded by at least a slight plastic de- 
formation. It appears that the stress concentration at the 
end of a slip band, stopped either by a grain boundary or 
by a precipitate, forms a small micro-crack of sufficient 
size that the conditions of Griffith® allow its further propa- 
gation. Equation (1) therefore gives a lower limit for the 
fracture stress. 


1J. A. Ewing and W. Rosenhain, on nt Soc. 65, 85 (1899). 

2? J. Muir, Phil. Trans. Roy. Soc. (1900). 

3 W. Rosenhain, J. Iron and ay ay 7: (2), 189 (1906). 

*C. Zener, Aneiasticity of Metals, A.1.M.E. in press. 

5 C, E. Inglis, Trans. Inst. Naval Architects 55, Part 1, 219 (1913). 

* E. Schmid and W. Boas, Kristallplastisitat (Berlin, 1935), p. 283. 

7R. A. Wilkins and E. S. Bunn, Copper and new Base Alloys 
(McGraw-Hill Book Company, Inc., New York, 1943). 

8A. A. Griffith, Phil. Trans. Roy. Soc. London 221, 163 (1920). 





The Low Energy 8-Spectrum of Cu®' 


HAROLD LEwis AND Davip BoxuM 
Department of Physics, University of California, Berkely, California 
January 12, 1946 


HE measurements of Backus! on the low energy end 

of the 8-spectra of Cu™, which branches to Zn and 
Ni, show a type of discrimination between electron and 
positron emission that is not contained in the Fermi theory. 
Between 5 kev and 50 kev the electron intensity is very 
closely proportional to the kinetic energy of the emitted 
electron, while the theory predicts a sensibly constant 
intensity over this range (within a few percent). However, 
the positron spectrum is in excellent agreement with the 
theory, so that it is unlikely that the electron behavior is 
caused by an idiosyncrasy of the apparatus. Also, the 
measurements of Tyler? on Cu contain this discrimina- 
tion, though, in Tyler’s case, it is masked by scattering 
from the source and support. 

That this constitutes a crucial test of the theory may be 
seen from the fact that the five relativistically invariant 
interaction hypotheses available in a theory of the Fermi 
type each predict the constant electron behavior. However, 
there are two linear combinations, the difference between 
the scalar and vector interactions, and the difference be- 
tween the tensor and pseudovector interactions, which 
will fit the observed electron spectrum, and not impair the 
agreement with the positron spectrum. These interactions 
essentially multiply the Fermi distribution by 1—mc*/E 
for electrons, and by 1+mc*/E in the case of positron emis- 
sion, where E is the energy of the emitted 8-ray, including 
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the rest energy. This is the only way in which these spectra 
can be fitted into a Fermi type theory. 

However, these linear combinations of interactions 
predict much smaller K-capture probabilities than are ob- 
served, and too strong an energy variation of lifetime for 
ordinary decay, so that neither is likely to be the domi- 
nating term in the majority of nuclei. If they do appear, it 
is probably in combination with another interaction term 
of the ordinary type, which latter must then be screened 
out by a selection rule in Cu“. Further development of the 


. question must await low energy experiments on other ele- 


ments, to determine whether or not the Cu® behavior is 
anomalous. 

In any case, it is clear that the low energy spectra of Cu 
force us to a choice of two severely restricted interaction 
forms, and the incompatibility of these interactions with 
other nuclear data indicates the important role of the low 
energy end of the spectrum in the investigation of the 
nature of @-decay. 

1J. Backus, Phys. Rev. 68, 59 (1945). 


2A. W. Tyler, Phys. Rev. 56, 125 (1939). 
3M. Fierz, Zeits. f. Physik 104, 553 (1937). 





Preliminary Investigation of Spark 
Discharges along Shock Waves 


E. E. HACKMAN 
Rouss Physical Laboratory, University of Virginia, Charlottesville, Virgina 
January 21, 1946 


OLLOWING is an outline of methods used and re- 

sults of preliminary investigations on spark discharges 

in the vicinity of shock waves formed in an air stream 
moving at supersonic velocities. 

Air was discharged through a 2” nozzle into the open at 
Mach numbers in the range of approximately 1.2-1.9. A 
metal rod }’’ diameter was supported in the air stream 
generally longitudinally along the axis of symmetsical 
flow, with its conical point at distances of 4” to 4” from 
the face of the nozzle. In some cases the rod was displaced 
outward from the center (axial) position. The rod was 
negatively charged to a potential of approximately 105 
volts by the usual method (Marx circuit) for recurrent 
sparking, discharging through the air stream to the nozzle 


Fic. 1. View of spark ane along shock cone. A-nozzle; B-rod 
i 


(negatively charged); C-path of discharge from rod to nozzle; D-shock 
f 


waves; E-direction of air stream. 
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Fic. 2. View of spark discharge along shock cone. A-nozzle; B-rod 
(negatively charged); C-path of discharge from rod to nozzle; D-shock 
waves; E-direction of air stream. 


face. In series with the high potential side of the circuit 
a spark gap was inserted to provide a light source for 
shadowgraph photographs. A representative photograph 
is shown in Fig. 1. As will be noted, the positions of the 
principal equipment are recorded together with the main 
spark discharge and the principal shock waves, in profile, 
of the air stream from the nozzle. Figure 2 is a sketch of 
the significant features. 

Of primary interest, as observed both visually and 
photographically, is the path of the main discharge. It will 
be noted that it proceeds generally directly ahead into the 
air stream, meeting the nozzle shock cone, and thence 
along the cone to the nozzle. In this instance the rod was 
displaced vertically approximately 0.15 inch to obtain a 
profile view. 

The foregoing phenomena were observed in numerous 
cases for various nozzle Mach numbers, in the range indi- 
cated above, and for varying distances of the point of the 
rod from the nozzle. The geometrical path varied as the 
position of the rod was shifted with respect to the shock 
cone. 

Further investigations are intended to explore in greater 
detail the observed phenomena with a view in part to 
determining the local conditions along the path of the 
discharge, and the conditions favoring such a path. 





A Relation between Characteristic Bond 
Constants and Electronegativities of 
the Bonded Atoms 


WALTER GorDY 
Department of Physics, Duke University, Durham, North Carolina 


February 1, 1946 
, I ‘HE relation, 


k=1.67N(xaxp/daz*)' +0.30, 


where & is the bond-stretching force constant in dynes/cm 
10-5, N the bond order, dag the bond length in Ang- 
stroms, and x4 and xg are the electronegativities of the 
bonded atoms A and B according to the Pauling scale,’ 
has been found to hold rather accurately for a large number 
of diatomic and simple polyatomic molecules in their 
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Fic. 1. Plot of bond-stretching force constant k ag 


CBr of CH:Br; CC bonds of C2H2, CsH«, CaHs, C 


ground states. Figure 1 shows & plotted as a function of 
(xaxe/dazs*)' for the bonds listed. 

When the same atom is involved in more than one of 
the bonds included, the electronegativity value given for 
it by Pauling was adjusted to give the best over-all fit 
with the relation. This adjustment merely required the 
extension of the values to an additional place as it was not 
found necessary in any instance to alter the last significant 
figure given by Pauling. This is remarkable as Pauling’s 
values were determined in a completely different manner, 
from the extra ionic resonance energies of bonds having 
partial ionic character. 

Identical bonds in different molecules are not repeated 
in the figure. The several CC bonds included all have differ- 
ent bond orders. The orders used for them are those re- 
cently determined by Mulliken, Rieke,and Brown? from the 
molecular orbital theory. The resulting good agreement of 
the CC bonds supports their rather surprising predictions 
that the usual CC “single” bond in hydrocarbons has 
appreciable double bond character, and that the usual CC 
“double” bond has appreciable triple bond character. 

The above relation does not apply to bonds containing 
atoms with unsatisfied valencies, to diatomic molecules in 
which both atoms have +1 valence, nor to polyatomic 


ainst N(xaxp/das*)? for: AsHa; BH of BsHe; BeO; Br:; 
3CCH, CeHe; CCI of CHsCl; CF of CHaF; CHa; CN of 
CH:NH: and of HCN; CO of CH3OH and of HsCO; CS of CHS and of CS2; CSe of C2HsSe; Cle; GeH«; HBr; 
_ A me 1Cl; I:; KBr; KCl; KI; NaBr; NaCl; Nal; N2; NHs; H2O0; OF 2; OCl:; PHs; HaS; HaSe; 
° 23 ™ 
Values for all bond constants and electronegativities and the sources from which they were obtained will be 
listed in a later publication. 


molecules in which there is appreciable interaction between 
the non-bonded atoms. Certain classes of them may, 
however, be fitted to the more general relation, 


k=aN(xaxp/das*)' +0, 


by a suitable choice of the parameters a and 6. 

A more detailed treatment of the relation with some of 
its applications will appear elsewhere. 

1L. Pauling, The Nature of the Chemical Bond (Cornell University 
Press, Ithaca, New York, 1939), p. 64 


?R_S. Mulliken, C. A. Rieke, and W. G. Brown, J. Am. Chem. Soc. 
63, 41 (1941). 





A Calculus of Finite Precision—-A Correction 


BENJAMIN LIEBOWITZ 
350 Fifth Avenue, New York, New York 
January 7, 1946 


N my paper on “A Calculus of Finite Precision,’' I 

introduced the concept of kinematical probability and 
with its aid derived from purely classical considerations 
Schroedinger’s equation for stationary states, or at least 
something very much the same as Schroedinger’s equation 
both in form and meaning. The assumptions underlying 
that derivation are: 
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(a) That the dynamical system is non-conservative 
but nevertheless has a potential energy function (contain- 
ing the time explicitly) so that its motions are describable 
in terms of the Hamilton-Jacoby equation. 

(b) That in a stationary state the system hovers around 
a definite energy E, so that the action function may be 
written as W=—Et+R(x;,t), which differs from the 
usual form in that R still contains the time explicitly; 
hence the H-J equation can be written 


(1/2m)2Z(8R/dx;?+ V-—E=—9dR/at, 


all masses being assumed equal. 
(c) Thus, —dR/dt is the difference between the instan- 


taneous energy (kinetic plus potential) of the system and - 


the mean energy E around which it hovers. It was assumed 
that in a stationary state the space average of this differ- 
ence is a minimum. 

The derivation, as given in the paper, appears to be 
exact, but on making the correction to which this note 
calls attention, it is seen to involve an approximation like 
that in the derivation given earlier in the paper of a 
Schroedinger-like equation for a planetary system in which 
the potential of the mutual interactions among the planets 
is a minimum. 

On making the substitution R= +k log U the integral 
to be minimized in order to minimize the average of 
—9@R/dt becomes that given by Eq. (3.8) of the paper. 
The function U? is then tentatively identified with the 
kinematical probability P; this leads to the Schroedinger 
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Eq. (3.10). In verifying this tentative identification, use 
is again made of the H-J equation, which leads to 


P=constant xU*+ f\(2m/p*)(@R/adt)(d@R/k). (A) 


In the paper itself the last term on the right of Eq. (A) was 
erroneously omitted. 
The integrand in Eq. (A) may be written 


[(E—E;)/REx](dR/ds)ds = ((E—E;)/REx |pds, 


where E as before is the mean energy around which the 
system hovers, E; is the instantaneous sum of its kinetic 
and potential energies, E; its kinetic energy, and » its 
total momentum. Thus the approximation is a good one 
so long as the difference between instantaneous and mean 
energies, or at least the integral of this difference, is small, 

To account, in a classical sense, for the existence of 
stationary states, I introduced in the last part of my paper 
the idea of an “atmosphere” of stationary wave fields 
surrounding an atomic nucleus and the concept of “phase. 
mechanical resonance” between the external electrons and 
said wave-fields. Reasoning on this basis, it appears that, 
for excited states, generally the above approximation is a 
valid one, i.e., that (E—E;)/E;, or at least its integral, is 
sufficiently small, since E—E; may take on positive and 
negative values with equal likelihood. For the normal 
state, however, and perhaps for long-lived metastable 
states also, the situation is more obscure. 


1 Benjamin Liebowitz, Phys. Rev. 66, 343 (1944). 
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Proceedings of the American Physical Society 


MINUTES OF THE MEETING AT THE UNIVERSITY OF CALIFORNIA, 
Los ANGELES, CALIFORNIA, JANUARY 12, 1946 


HE 269th meeting of the American Physical 
Society was held at the University of Cali- 
fornia, Los Angeles, on Saturday, January 12, 
1946. The attendance at the meeting varied 
between 100 and 125, possibly the largest number 
at any similar gathering at Los Angeles. Because 
of the continuing critical condition in regard to 
travel, on the Pacific Coast, almost the entire 
audience came from Southern California. On the 
other hand, there were present, for the first time 
since Pearl Harbor, many persons who had been 
and still were engaged on research projects under 
the control of the Armed Forces. These persons 
came from a variety of locations, extending from 
Los Alamos to San Diego. 

Nearly 90 members and guests of the Society 
had luncheon together at the cafeteria in Kerck- 
hoff Hall. In view of the fact that prior arrange- 
ments had been made for a special luncheon for 
not over 50, the satisfactory service rendered on 


short notice to nearly twice that number repre- 
sented a remarkable achievement. 

The morning session opened at 10:00 a.m. and 
adjourned at noon. It was concluded by a most 
interesting invited paper by Dean Vern O. 
Knudsen of the University of California, Los 
Angeles, entitled ‘Propagation of Sound in Air— 
Absorption and Fluctuations.’’ The afternoon 
session opened at 2:00 p.m. and the formal pro- 
gram concluded about 4:00 p.m. It was followed 
by the presentation of a Navy film, with sound, 
of the use of rockets in the war, especially rockets 
shot from airplanes. This film was furnished 
through the courtesy of Dr. W. V. Houston of 
the California Institute of Technology. Abstracts 
of the 15 contributed papers follow. Only No. 9 
was read by title. 


R. T. BiIrGE 
Local Secretary for the Pacific Coast 


Invited Paper 


Propagation of Sound in Air—Absorption and Fluctuations. 
VERN O. KNUDSEN, University of California, Los Angeles. 


ABSTRACTS OF CONTRIBUTED PAPERS 


1. Forbidden Line Source. JosEpH KAPLAN, University ~ 2. Rotating Ring Sine-Wave Generator, Especially 


of California, Los Angeles—The most intense laboratory 
source of the forbidden lines of atomic nitrogen and oxygen, 
which the author has been able to produce, will be demon- 
strated to the audience. It is an afterglow, which is pro- 
duced in a mixture of nitrogen and oxygen at a pressure of 
approximately 100 mm, and the green auroral line can be 
readily observed with a direct vision spectroscope. The 
forbidden lines are relatively weak during the first few 
seconds following the interruption of the discharge, and 
the allowed radiations are relatively weak during the 
remainder of the approximately two minutes duration of 
the afterglow. During the first part of the afterglow, the 
spectrum closely resembles that of the polar aurora, while 
during the last part, it resembles that of the light of the 
night sky. The importance of forbidden lines to the astro- 
physicist, the physicist of the upper atmosphere, and the 
photochemist, and the possible uses of this source will be 
discussed. The mechanism of the excitation of the green 
auroral line will also be discussed. 


Suitable for Measurement of Magnetic Intensities and 
Their Variations. S. J. BARNETT, University of California, 
Los Angeles, and California Institute of Technology, Pasa- 
dena.—A thick conducting toroidal ring or short-circuited 
coil rotates at f r.p.s. about a symmetrical diameter A 
normal to its axis in a magnetic field with intensity V 
normal to A. The center of the toroid coincides with that 
of a surrounding coil C, preferably a Helmholtz pair, with 
axis B normal to V and A. If the constant G of C and the 
intensity V are uniform, a harmonic e.m.f. with frequency 
2f and amplitude EZ is developed in C such that £ is 
proportional to V and G, and, when eddy currents are 
negligible, to the mean square area of the toroid and to 
f?/Z, where Z is the toroid’s impedance. Largely on 
account of eddy currents an accurate calculation of E is 
impracticable for a thick ring. The e.m.f. for 245 gausses 
and 60 r.p.s. obtained with the ring used was 31 volts r.m.s. 
With short-circuited coils at 268 gausses and 60 r.p.s., the 
calculated and experimental values were 51.7 and 52.0 
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volts r.m.s., respectively. The device has great advantages 
over rotating coils with commutators or slip rings. It is 
simple, robust, and reliable, and should serve as an 
excellent vertical intensity variometer in observatories, 
where the ordinary types of instrument are troublesome. 
Practically any sensitivity desired can be attained by 
suitable design and amplification. 


3. Matrix Elements for Positron Decay. W. A. Fow.er, 
California Institute of Technology.—The original data’ on 
certain positron spectra have been re-examined with the 
purpose of estimating the relative intensities of transitions 
to the excited and ground states of the daughter nuclei. 
These states constitute doublets differing only in the 
relative orientation of the orbital (1 or 2) and spin (4) 
momenta. If Fermi energy distributions are assumed, the 
marked curvature of all the Kurie diagrams ‘indicates 
clearly that in all cases the matrix element of the transition 
to the excited state is the greater. The results are given in 
Table I. The results of N™, O'5, and F!7, but not for C", 


TABLE I, 








Transitions 


1:6. P32 —P 3/2; P32 —Pij2 
as t : Pij2 —P ja; Pij2—P 3/2 
1:2. : 


Ei E: 
(Mev) (Mev) |G|? Intensity 





Pij2—P 12; Pij2—P 3/2 
Daj: —D3;2; D3j2 —Ds;2 








are in agreement with the matrix elements given by 
Wigner? if the square of the Gamow-Teller matrix element 
is twice the square of the Fermi matrix element. 


1 Fowler, Delsasso, and Lauritsen, Phys. Rev. 49, 561 (1936). 
2 E. Wigner, Phys. Rev. 56, 519 (1939). 


4. Some Electrical Characteristics of Reflex Tubes. H. 
V. NEHER, California Institute of Technology.—Following 
the method outlined by J. R. Pierce’ for setting up and 
solving the electronic equations for the reflex oscillator 
tube, expressions are derived for the efficiency and elec- 
tronic tuning when the transit angle across the gap in the 
cavity has its optimum value and the tube is delivering 
maximum power to the load. It is shown that where space 
charge effects are not important this optimum transit 
angle is approximately 135°. It is also shown that when 
the beam current is such that the tube is operating at 
maximum efficiency, maximum power is delivered to the 
load when the Q of the cavity due to the load is decreased 
to half of its unloaded value. 


1J. R. Pierce, Proc. I. R. E. p. 112 (Feb. 1945). 


5. A Variational Method for Linear Dissipative Aniso- 
tropic Elastic Systems. GLEN D. Camp, University of 
California, Division of War Research, San Diego, California. 
—A boundary value problem, governing a quite general 
elastic system, is defined. Although limited to steady-state 
drive (or free vibrations), it is applicable to transient 
drives by the usual method of Fourier superposition. A 
variational principle, rigorously equivalent to this bound- 
ary value problem, is given. It is amenable to the direct or 
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semi-direct method, and thus takes advantage, in th 
choice of trial functions, of any qualitative or semj. 
quantitative knowledge of the system. It furnishes practi. 
cally useful approximations, with reasonable labor, to the 
solution of problems which have not been treated success. 
fully by other methods. An example is worked out jy 
detail and the solution is interpreted physically. This 
work represents one of the results of research carried oyt 
by the University of California, Division of War Research, 
under contracts with the Office of Scientific Research ang 
Development, Section 6.1, National Defense Research 
Committee, and with the Bureau of Ships, Navy Depart. 
ment. 


6. Fundamental Properties of the Vacuum Switch. R. 
KOLLER, Department of Electrical Engineering, University 
of California, Berkeley (Introduced by L. B. Loeb).—An 
exhaustive study of the basic mechanism affecting the 
direct-current operation of the vacuum switch has bees 
undertaken. It is found that the switch operates below 
10-* mm Hg pressure. The evaporated and sputtered 
material with appropriate geometrical design factors acts 
to reduce the pressure due to gases evolved from the 
electrodes by a combination of getter and diffusion-pump 
like action. This acts to maintain a constant working 
pressure. Strict outgassing of the electrodes is not essential 
but an adequate forming treatment is necessary. Losses of 
various electrode materials have been measured under 
different switching conditions and are small. Tungsten 
gave least losses, aluminum most. The mechanism active 
consists of a series of intermittent arc bursts of microsecond 
duration interspersed with extremely short periods of 
extinction. The mechanism is initiated in most cases by 
field émission yielding to magnetically self-focusing beams 
of the Bennett! type, which, leading to an arc, ultimately 
are extinguished, as indicated in Bennett's paper. This 
mechanism is confirmed by the character of losses and 
microscopic pitting of the cathode material, which deposits 
on the anode and on the walls. 


1W. H. Bennett, Phys. Rev. 45, 890 (1934). 


7. Multivibrator Theory. S. C. SNnowpon, California 
Institute of Technology (Introduced by C. D. Anderson).— 
The transient response of a symmetrical multivibrator 
has been analyzed neglecting the effect of shunt capacities 
and grid current. For an arc-tangent approximation to the 
transfer characteristic the equation connecting the two grid 
voltages is found to be du/dv = (u—bv/1+0*)/(v—bu/1+4?). 
This equation has been analyzed for its critical points and 
integrated graphically for suitable values of b. An approxi- 
mate analysis shows that for large times u->—v as is 
observed. The interesting initial condition is « almost 
equal to v. Here the initial transient is very long compared 
with the multivibrator period and the initial sign of « 
may be adjusted by making u slightly greater or slightly 
less than v. 


8. Excitation Energy of Excited Li’. S. Rusin, California 
Institute of Technology.—A measurement of the excitation 
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energy of the low excited state of Li’ has been made by 
comparing the absorption coefficients in lead of the gamma- 
rays from the excited Li’ produced in the decay of Be’ 
with the annihilation radiation from the positrons emitted 
by N¥. The absorption coefficient of the annihilation 
radiation was measured (in 1941) using three samples of 
N'3 produced by deuteron bombardment of carbon. The 
geometry of the absorption measurements was identical 
with that used for the Be’ source, for which a value of 
p=1.635+0.015 cm™ was reported.! It was since found 
that a correction had been omitted, so that the corrected 
value becomes w= 1.570.015 cm. The corrected value 
obtained for the annihilation radiation is ws=1.40+0.03 
cm~!. From the two coefficients, and Read’s* absorption 
coefficient curve in this energy range, the energy of the 
Li’ radiation is 476+10 Kev. This agrees with Zlotowski 
and Williams’* value of 485+5 Kev, determined from the 
relative ranges of the recoil electrons produced by the two 
radiations. 


1S. Rubin, Phys. Rev. 59, 216A (1941). 
2 J. Read, Proc. Roy. Soc. 152, 402 (1935). 
31. Zlotowski and J. H. Williams, Phys. Rev. 62, 29 (1942). 


9. Diabatic Flow of a Compressible Fluid.* Bruce L. 
Hicks, Aircraft Engine Research Laboratory, Cleveland, 
Ohio.—The steady flow of inviscid, compressible fluids 
containing heat sources has been examined using a reduced 
velocity vector W=V/V; similar to Crocco’s vector. In 
such diabatic flow, the variation of V;, the limiting velocity, 
along a streamline is determined by the heat source 
function g which, in burning gases for example, can vary 
throughout the flow. The basic gas dynamical equations 
in the W language reduce to two equations containing 
only W, g, and :, the local stagnation pressure. It is found 
from one (vector) equation that the rate of variation of p; 
along streamlines is proportional to g and, in a direction 
perpendicular to streamlines, to a vorticity, |w| =|VxXW]|. 
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The second (scalar) equation indicates that f(W)g acts 
like a fluid source in the vector field (1—W*)”7"W, and 
that the rate of variation of stream tube area when local 
sonic velocity is reached does not vanish as it does in 
adiabatic flow but is proportional to g. Of the several 
types of irrotational diabatic flow that have been con- 
sidered, one is barytropic, and all lead to linear or quasi- 
linear partial differential equations of the second order for 
the potential function. 


* Read by title. 
1L. Crocco, Z a. M. M. 17, 1 (1937). 


10. A Computer for Solving Linear Simultaneous Equa- 
tions. Ciirrorp E. Berry, Consolidated Engineering 
Corporation, Pasadena (Introduced by Harold W. Wash- 
burn).—The classical iterative method of solving simul- 
taneous equations consists of solving the first equation for 
the first unknown, using arbitrary values of the other 
unknowns, then solving the second equation for the second 
unknown, using the new value for the first unknown and 
the same arbitrary values of the others, etc. This process 
is repeated until further cycles give no appreciable change 
in the unknowns. In a great many problems of practical 
importance, the equations can be readily arranged so that 
the above process converges. An electrical circuit contain- 
ing pairs of potentiometers, one energizing the other, is 
used to represent a series of pairs of products, a;;x;; and 
the sum of these products is obtained by a parallel ‘‘addi- 
tion” network. In this way an analog of one equation is 
set up. In order to set up a number of equations, potenti- 
ometers representing different sets of coefficients, a,j, are 
switched into the same set of unknowns, x;. The accuracy 
of the computer described is such that in using data 
accurate to three significant figures, the errors introduced 
by the computer are negligible. Problems involving twelve 
equations are solved from four to seven times faster than 
by conventional methods. 








11. Rotating Magnetic-Rod Generator Suitable for the 
Measurement of Small Magnetic Intensities and Their 
Variations, and for Other Purposes. S. J. BARNETT, Uni- 
versity of California, Los Angeles, and California Institute 
of Technology, Pasadena.—This generator is formed by 
substituting for the toroid of that described in paper 2 a 
group of thin, straight, parallel rods, preferably of perm- 
alloy or other material of high permeability and small 
retentivity, with their axes normal to the axis of rotation. 
Two Helmholtz pairs were used, C, that used with the 
toroid, and D, with a constant 0.015 as great. In weak 
fields, up to 4 gausses for an 8-rod system and 8 gausses 
for an 18-rod system, E is harmonic with frequency 2f and 
proportional to V, G, and f. In the region of proportionality 
with V, and at f=60, E is 3.8 and 5.2 times as great as for 
the ring. As V increases and saturation approaches, one 
half-cycle extends in breadth and becomes approximately 


Supplementary Programme 








a half-cycle at frequency f, while the other becomes a sharp 
peak on a narrow base—with oscillations present when C 
is used, on account of inductance and capacity in the 
circuit. Related experiments were made with a system of 
2 steel rods. Oscillograms and curves are given showing 
how the results depend on V and f. In strong fields the 
instrument may serve as an impulse generator, or as a 
generator of numerous harmonics; in very weak fields it 
can serve as an excellent variometer. 


12. Impedance Representation of Tangential Boundary 
Surface Forces. GLEN D. Camp, University of Cali- 
fornia, Division of War Research, San Diego, California.— 
The importance of tangential boundary surface forces in 
some systems is noted. The tensor formulation of a general- 
ized impedance boundary condition, involving a tangential 
as well as a normal impedance, is given. The tangential 
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impedance, seen by a plane working into a viscous fluid, 
is calculated. This work represents one of the results of 
research carried out by the University of California, 
Division of War Research, under contracts with the 
Office of Scientific Research and Development, Section 6.1, 
National Defense Research. Committee, and with the 
Bureau of Ships, Navy Department. 


13. A Critical Examination of Sound Velocity Measure- 
ments. R. W. LEonarp, University of California, Los 
Angeles.—It is the purpose of this paper to discuss some 
of the factors limiting the precision of sound velocity 
measurements for freely progressive waves in air. The 
work of Michelson and Hebb, D. C. Miller, and others 
will be reviewed briefly and the limitations of the methods 
discussed. The general conditions which must be satisfied 
for precise measurements will be established. 


14. A Method for the Precision Measurement of the 
Velocity of Sound in Air. R. W. LEonarp, University of 
California, Los Angeles.—It is the purpose.of this paper to 
describe the apparatus used in the measurement of the 
velocity of sound in air. The method involves the use of 
a source generating spherical waves, a movable micro- 
phone, and a means of comparing the phase of the sound 
pressure at the microphone with the motion of the sound 
source. Points of equal phase are located along the axis of 
symmetry of the sound source. These points are separated 
by one wave-length. By measuring over a distance equal 





AMERICAN PHYSICAL SOCIETY 


to several hundred wave-lengths an accurate value May 
be obtained by proper averaging. The results obtained by 
this method will be presented. 


15. Acoustic Reflection from Triplanes, Spheres, ang 
Disks. C. J. BURBANK, University of California, Division 
of War Research, U. S. Navy Radio and Sound Laboratory, 
San Diego.—This work represents one of the results of 
research carried out by the University of California, Divi. 
sion of War Research, under contracts with the Office of 
Scientific Research and Development (Section 6.1 National 
Defense Research Committee) and with the Bureau of 
Ships, Navy Department. The acoustic back reflection 
from a triplane (an object made of three mutually per. 
pendicular square plates with coinciding centers) is com- 
pared with the back reflection from a sphere and a disk, 
The reflection from the triplane is measured for all angles 
in the following three planes: (1) A plane coincident with 
one of the planes, (2) a plane making equal angles with 
two of the planes, (3) a plane making equal angles with all 
three of the planes. The measurements are made by trans- 
mitting a continuous supersonic wave toward the object 
and measuring the level of the signal reflected to a receiver 
mounted just below the source. Cross talk in the system is 
reduced to a negligible level by the use of balanced lines 
and directional units. The validity of the measurements is 
checked by comparing the measured sound level at the 
receiver with the calculated level for a ten-inch disk, a 
square plate, and a three-foot sphere. 














